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STATEMENT  OF  WORK 


1  Use  a  Phase  Doppler  Particle  Analyzer  to  measure  the  spray  size  and  velocity  distribution  from  shear 
coaxial  injectors  under  steady  conditions  as  a  function  of  the  relative  gas/liquid  velocity.  Liquid  oxygen  will  be  used 
for  the  central  liquid  flow  and  either  nitrogen  or  helium  for  the  outer  annular  gas.  The  possibility  of  using  gaseous 
hydrogen  for  the  final  series  of  tests  will  be  examined.  Use  a  laser  diode  LDV  with  a  transparent  injector  to 
examine  the  recirculation  region  at  the  base  of  the  injector.  The  following  parameters  will  be  investigated: 

a)  Chamber  pressures  up  to  1500  psia. 

b)  Injector  design  in  terms  of  flow  area  ratio,  LOX  post  tip  land  width  and  liquid  flow  recess 
distance. 

2.  Use  the  measured  spray  size  and  velocity  distributions  as  input  to  Penn  State  numerical  codes  to  examine 
the  effects  of  atomization  on  predicted  motor  stability  and  to  establish  an  experimental  database  to  provide  guidance 
towards  the  development  of  an  atomization  model.  Compare  results  to  hydrogen  temperature  ranqnng  test  results 
where  available  (NASA  Lewis,  Rocketdyne). 

3.  Examine  the  effect  of  injector  pressure  drop  on  propellant  feed  system  caused  combustitm  instabilities. 
Experimentally  characterize  the  effects  of  operating  conditions  and  injector  design  on  injector  pressure  drop. 

4.  Measure  droplet  sizes  and  velocities  under  combusting  supercritical  conditions  using  a  PDPA.  Investigate 
the  effects  of  relative  gas/liquid  velocities  and  injector  geometry.  Attempt  to  examine  the  recirculation  region  at 
the  base  of  the  LOX  post  by  means  of  a  laser  diode  LDV  system.  Examine  any  effects  of  either  the  combustion 
process  or  the  supercritical  conditions  on  the  injector  pressure  drop. 

5.  Use  the  measured  droplet  sizes  and  velocities  under  combusting  supercritical  conditions  as  input  to  the 
chamber  stability  codes  developed  at  Penn  State  as  well  as  provide  the  guidance  for  the  development  of  injector 
atomization  models. 
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PROGRESS  AND  ACCOMPLISHMENTS 
Experimental  Facility 

Cold  flow  experiments  with  liquid  nitrogen  were  conducted  at  the  Cryogenic  Combustion  Facility  of  The 
Pennsylvania  State  University.  The  facility  was  designed  for  hot  fire  testing  of  subscale  rocket  engines  using  liquid 
oxygen  and  gaseous  hydrogen  as  the  propellant  and  can  accommodate  other  propellant  combmations  under  hot  fire 
or  cold  flow  conditions.  Facility  performance  is  rated  at  0.45  kg/s  of  liquid  oxygen  and  0.11  kg/s  of  gaseous 
hydrogen  for  chamber  pressures  up  to  approximately  100  bars.  The  propellant  feed  system  is  divided  into  two  gas 
supply  lines  for  primary  and  secondary  gas  flows  and  a  cryogenic  liquid  supply  line  which  includes  a  114  liter 
vacuum  jacketed  storage  tank,  tank  pressurization  system,  and  liquid  nitrogen  cooling  jacket.  Flow  metering  was 
accomplished  with  critical  flow  venturies  fitted  with  pressure  transducers  and  thermocouples  upstream  and 
downstream  of  the  venturi. 

The  pressure  vessel  was  designed  for  cold  flow  operation  at  pressures  up  to  100  bars.  The  interior 
dimensions  of  4  inches  square  by  9  inches  long  permit  the  investigation  of  the  conqrlete  spray  field  while  avoiding 
window  contamination  by  impinging  droplets.  There  are  four  gas  ports  around  the  injector  to  provide  additional 
gas  for  chamber  pressurization  and  window  purging.  Pressure  and  temperature  are  measured  in  the  chamber  as  well 
as  in  the  injector  plenum  region  upstream  of  the  fuel  annulus  and  in  the  cryogenic  liquid  supply  line.  To  facilitate 
axial  and  radial  scans  of  the  injector  spray,  the  entire  chamber  was  mounted  on  a  2-D  translation  table.  Positioning 
in  the  radial  direction  was  performed  by  a  stepping  motor  while  axial  positioning  required  manual  operauon. 

Photographic  Studies  of  Shear  Coaxial  Injectors  Using  Liquid/Gaseous  Nitrogen 
Two  flow  visualization  techniques  were  pursued  to  examine  the  overall  structure  and  break-up  processes 
of  the  LNj/GNi  spray  under  varying  conditions.  The  first  approach  to  visualize  the  spray  involved  the  imaging  of 
Mie  scattered  light  from  a  thin  laser  sheet  passing  through  the  axis  of  the  spray.  A  frequency  doubled  Nd:YAG 
laser  (532  nm)  pulsed  at  10  Hz  with  a  pulse  duration  of  10  ns  provided  the  light  source  while  a  lens  pair  consisting 
of  a  spherical  and  cylindrical  lens  formed  the  laser  sheet.  The  scattered  light  was  collected  at  90  degrees  to  the 
sheet  with  a  35  mm  Nikon  camera.  To  eliminate  ambient  light,  a  1  nm  bandpass  filter  centered  at  532  nm  was 
installed  in  front  of  the  camera  lens. 
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Although  the  laser  sheet  photogr^hs  provided  useful  images  of  the  overall  structure  of  the  dense  liquid 
core  region,  minimal  information  was  available  on  the  distribution  of  LNj  droplets  within  the  spray  field.  Since 
this  information  is  of  value  in  determining  optimal  regions  for  PDPA  measurements,  a  modified  q>proach  was  then 
undertaken  wherein  a  strobe  light  replaced  the  pulsed  laser  sheet  for  spray  illumination.  The  strobe  light  was 
directed  toward  the  back  side  of  the  spray,  off  center  from  the  camera  in  the  direction  of  the  injector  face.  With 
this  orientation  both  the  liquid  core  region  and  regions  of  high  number  density  and  large  droplet  concentrations  could 
be  detected.  The  flash  duration  of  the  strobe  was  3  /is  at  a  frequency  of  10  Hz. 

The  objective  of  the  LN^/GNj  tests  was  to  characterize  the  effects  of  chamber  pressure  and  liquid  to  gas 
mixture  ratio  on  the  atomization  and  v^rization  of  liquid  nitrogen  serving  as  a  simulant  to  liquid  oxygen.  The 
baseline  Himpnsinns  of  the  shear  coaxial  injector  were  derived  from  the  dimensions  of  the  prototjqie  injector  element 
of  the  fuel  prebumer  to  the  Space  Shuttle’s  main  engine.  A  second  LX)X  post  was  fabricated  with  a  7*  t^»er  at  the 
tip  to  mafrh  the  shape  of  the  production  injector  element.  Pressure  and  temperature  were  measured  in  the  fuel 
plenum  region  upstream  of  the  fuel  annulus  and  in  the  cryogen  supply  line  at  the  entrance  to  the  LOX  post.  For 
the  elevated  pressures  required  of  LNj  testing,  a  pressure  chamber  was  constructed  to  operate  at  pressures  up  to 
10.0  MPa  and  to  provide  optical  access.  Flow  visualization  tests  were  performed  with  the  nontapeted  LOX  post 
at  chamber  pressures  between  2.2  and  4.1  MPa  and  at  two  liquid  to  gas  mixture  ratios  of  0.80  and  1.58.  To 
determine  if  any  observable  differences  in  liquid  atomization  exists  between  the  two  LOX  post  designs,  several 
stroboscope  tests  were  conducted  with  both  LOX  posts  operating  at  the  same  mixture  ratio  and  at  chamber  pressures 
from  2.2  to  3.29  MPa. 

The  fust  series  of  tests  illustrated  the  effect  of  chamber  pressure  on  the  overall  spray  morphology  of  the 
coaxial  injector  using  the  nontapeted  LOX  post.  Details  are  given  in  Ref.  1  (Appendix  A)  and  are  briefly 
summarized  here.  The  photogrsq)hs  were  imaged  with  the  laser  sheet.  It  is  expected  that  in  the  case  of  large, 
dispersed  sprays  differences  in  the  gas  temperature  may  influence  the  behavior  of  the  spray  downstream  of  the  dense 
liquid  core  region  but  should  not  affect  significantly  the  breakup  behavior  of  the  liquid  jet  or  the  distnbution  of 
droplets  around  the  jet,  since  the  conditions  in  this  region  are  better  represented  by  the  gas  injection  temperature. 
In  general,  the  photographs  depict  an  overall  contraction  of  the  large  scale  liquid  structures  (liquid  core  and  detached 
ligaments)  and  droplet  flow  regions  with  increasing  ambient  pressure.  Regions  of  droplet  flow  are  most  evident 
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downstream  of  the  liquid  structures,  while  very  little  droplet  production  is  observed  along  the  liquid  column  close 
to  the  injector  exit.  With  rising  chamber  pressure  the  gas  density  enhances  the  breakup  of  the  jet,  but  reduces 
secondary  atomization  due  to  a  decrease  in  the  slip  velocity  between  the  gas  and  dropled  The  large  gas  mass  flux 
appears  to  have  rapidly  vaporized  droplets  surrounding  the  core  except  for  those  in  the  shadow  of  the  LN^  jet,  a 
region  characterized  by  lower  gas  temperatures  and  velocities. 

In  a  second  series  of  tests  the  gas  mass  flow  rate  was  decreased  to  approximately  half  of  the  value  of  the 
previous  tests,  while  the  chamber  pressure  was  varied  from  2.34  MPa  to  2.79  MPa.  The  photographic  results  of 
these  tests  were  obtained  as  well  with  the  laser  sheet  visualization  technique.  Overall,  the  reduction  in  gas  mass 
flow  rate  resulted  in  a  larger,  more  dispersed  spray  indicating  a  strong  effect  of  mixture  ratio  on  the  atomization 
and  vaporization  of  liquid  nitrogen.  The  behavior  of  the  spray  with  respect  to  increasing  diambcr  pressure  foUowed 
the  same  trend  as  before,  in  that  the  LN^  spray  contracts  in  size  as  the  chamber  pressure  increases. 

A  third  series  of  tests  were  performed  with  the  tapered  LOX  post  using  the  stroboscope  to  visualize  the 
spray.  The  stroboscopic  image  provides  more  detail  on  the  structure  of  the  liquid  core  and  detached  ligaments  than 
the  laser  sh«yt  images  but  less  information  on  the  droplet  flow  region.  Contradictory  to  the  previously  observed 
contraction  of  the  spray  with  increasing  pressure,  the  liquid  core  breakup  length  significantly  increased  for  a 
chamber  pressure  greater  than  the  critical  pressure  (3.4  MPa).  This  behavior  may  be  attributed  to  the  drop  in 
relative  velocity  between  the  gas  and  liquid  at  very  high  chamber  pressures,  which  may  surpass  any  enhancements 
to  liquid  breakup  due  to  increased  gas  density,  and  may  be  related  as  well  to  the  increase  in  droplet  evaporation 
times  observed  in  the  supercritical  regiorf.  Stroboscope  tests  with  the  nontapered  LOX  post  under  similar  operating 
conditions  revealed  the  same  behavior  for  pressures  above  and  below  the  critical  point  but  indicated  little  differences 
in  spray  morphology  between  the  two  LOX  post  designs. 

PDPA  Studies  of  Shear  Coaxial  Injectors  Using  Water/ Air 

The  injector  was  operated  at  6.0  g/s  of  water  and  3.7  g/s  of  air  using  both  the  tapered  and  non-tapered 
LOX  post  recessed  inside  the  annulus.  Data  on  droplet  mean  velocity,  mean  size,  number  density  and  volume  flux 
were  obtained  as  a  function  of  radial  position  at  5  cm  downstream  of  the  injector  face.  The  spray  data  depicted 
noticeable  differences  between  the  two  sprays.  The  velocities  on  the  injector  centerline  are  approximately  equal 
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for  both  configurations  but  depart  outward  in  the  spray.  A  more  r^id  drop  in  axial  velocity  for  the  tapered  post 
suggests  that  the  droplets  are  traveling  faster  in  the  radial  direction  than  in  the  axial  direction.  With  the  tapered 
post  smaller  Sauter  mean  diameters  are  found  on  the  centerline  and  increase  radially  outward  in  the  spray .  Overall, 
the  tapered  LOX  post  produced  smaller  droplets  within  the  region  of  the  spray  probed  with  the  PDPA.  In  the  case 
of  the  tapered  post  the  distribution  of  volume  flux  is  rather  uniform  throughout  the  spray  as  compared  to  the  steep 
distribution  obtained  for  the  nontapered  post.  The  occurrence  of  this  distribution  is  due  to  the  opposite  trends 
between  the  mean  velocity  and  the  number  density  and  drop  size.  That  is,  the  steep  velocity  profile  compensates 
for  the  hollow  cone  distribution  in  drop  size  and  number  density  to  create  the  uniform  volume  flux  distribution. 
Based  on  this  distribution  the  mixing  efficiency  of  the  tapered  post  is  higher  than  that  of  the  nontapered  post. 

Recirculation  Region 

The  original  injector  design  with  the  straight  LOX  post  was  modified  in  order  to  allow  LDV  access  into 
the  recess  region  by  inserting  a  Pyrex  tube  concentrically  around  the  LOX  post,  such  that  the  fud  annulus  was 
defined  by  a  transparent  boundary  1.02  mm  thick.  The  inner  diameter  of  the  Pyrex  fuel  annulus  was  5.00  mm. 
All  measurements  were  taken  using  Aerometrics’  Phase  Doppler  Particle  Analyzer  (PDPA)  with  the  droplet  diameter 
m>.a<airpTnent  function  disabled.  Measurements  of  mean  axial  velocity  were  taken  inside  the  injector  recess  region 
as  well  as  downstream  of  the  injector.  Water  was  used  as  the  LOX  simulant  and  either  air  or  gaseous  nitrogen  for 
the  gaseous  fuel  simulant.  The  air  flow  rate  was  1.8  g/s  and  the  water  flow  rates  were  either  1.9  or  4.5  g/s. 
Downstream  of  the  injector  face  only  positive  velocities  were  measured.  It  was  found  that  the  velocity  drops  to  a 
negative  value  (reverse  flow)  approximately  0.5  mm  downstream  of  the  LOX  post  base  before  becoming  positive. 
Larger  negative  values  of  the  velocity  were  observed  for  lower  liquid  flow  rates.  Thus  it  aj^jears  that  gas 
recirculation  does  exist  in  the  recess  region  of  a  shear  coaxial  injector  element  and  that  the  strength  of  the 
recirculation  is  a  function  of  the  operating  conditions. 


Injector  Coupling 

Some  previous  experimental  evidence  has  indicated  that  the  spontaneous  instability  condition  for  a 
hydrogen/oxygen  rocket  was  determined  by  the  hydrogen  pressure  drop  through  the  injector*-*.  Bdow  a  critical 
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minimum  value  for  the  pressure  drop,  chamber  oscillations  evidently  can  couple  to  the  propellant  feed  system, 
raining  a  combustion  instability  in  the  chamber.  However,  the  value  of  this  critical  pressure  drop  was  a  function 
of  the  injector  and  chamber  design. 

A  simple  model  based  on  the  Rayleigh  criterion  was  constructed  for  the  flow  dynamics  of  several  processes 
in  the  system  to  determine  the  effect  of  various  design  and  operating  parameters  on  combustion  stability.  The 
mechanisms  included  are  the  fuel  flow  through  the  injector,  the  vaporization  of  oxidizer  in  the  chamber,  and  the 
flow  of  exhaust  gases  through  the  nozzle. 

The  Rayleigh  criterion  states  that  a  wave  will  grow  if  heat  or  mass  is  added  in  phase  with  pressure. 
Conversely,  the  wave  will  damp  if  the  addition  is  out  of  phase.  If  several  mechanisms  are  contributing  to  the  heat 
or  mass  release,  then  the  wave  growth  is  determined  by  the  net  in-phase  or  out-of-phase  heat  or  mass  addition.  By 
assuming  that  both  the  heat  or  mass  addition  and  the  pressure  oscillations  can  be  represented  by  a  perturbation  added 
to  their  steady  swte  values,  linear  models  of  the  dynamic  response  of  the  processes  may  be  obtained.  A  response 
factor  is  defined  as  the  real,  or  in-phase,  part  of  the  mass  flow  peituibation  with  tespea  to  the  diamber  pressure 
perturbation.  The  response  factor  can  be  either  positive  or  negative  dqiending  on  the  phase  angle  between  the 
pressure  and  flow  perturbations. 

A  lumped-element  model  was  developed  similar  to  Reference  6  to  determine  the  dynamic  response  of  the 
fuel  flow  through  the  injector.  The  flow  through  the  dome,  the  orifice  flow  into  the  entrance  of  the  annulus  and 
the  flow  through  the  annulus  are  represented  by  perturbation  forms  of  the  mass  and  momentum  balances.  The 
pressure  at  the  exit  of  the  annulus  was  assumed  to  equal  the  mean  chamber  pressure.  These  equations  may  be 
combined  to  give  the  transfer  function  and  inverting  this  transfer  function  gives  the  amplitude  and  phase  shift. 

The  injector  response  factor  was  calculated  as  a  function  of  frequency  and  gaseous  fuel  temperature  for 
the  experimental  conditions  experienced  during  the  cold-flow  LNj/GNj  experiments.  During  these  experiments,  high 
frequency  pressure  oscillations  were  measured  in  the  chamber  and  the  injector  fuel  plenum.  The  recorded  high 
frequency  pressure  oscillations  were  spectrally  analyzed  to  determine  the  frequencies  and  magnitudes  of  any  acoustic 
oscillations  within  either  the  chamber  or  the  injector  fuel  plenum.  Oscillations  at  the  same  frequency  in  both  the 
chamber  and  the  plenum  were  found  at  two  frequencies.  It  was  found  that  at  the  frequency  where  the  injector  was 
more  strongly  damping,  the  experimentally  measured  ratio  of  the  pressure  oscillation  amplitude  in  the  chamber  to 
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the  amplitude  in  the  plenum  was  less. 


Subscale  Hot-Fire  Rocket  Combustion  Chamber 

Injector  faceplates  of  several  existing  and  proposed  liquid  propellant  rocket  engines  were  reviewed  to 
ascertain  the  local  grouping  of  injector  elements  around  a  single  element.  Although  combustion  chambers  are 
generally  cylindrical  in  design,  the  immediate  boundary  of  a  single  injector  element  is  more  accurately  represented 
by  a  hexagon  formed  by  six  adjacent  injector  elements.  A  hexagon  enclosing  a  2.82  cm  circle  was  selected  for  the 
chamber  interior  cross-section  allowing  flat  windows  for  optical  access  for  laser  based  diagnostic  devices.  With 
the  six  sided  polygon  flat  windows  can  be  mounted  onto  the  chamber  at  opposite  sides  or  at  an  inclination  of  60* 
to  the  principal  axes  of  the  windows. 

nie  chamber,  as  illustrated  in  Fig.  1,  is  modular  in  design  aUowing  for  simple  variation  of  the  chamber 
length,  location  of  optical  access,  injector  design  and  throat  diameter.  The  individual  chamber  segments  include 
the  injector  segment,  window  segment,  igniter  segment,  blank  segment  and  nozzle  segment.  With  aU  segments  m 
place  the  interior  length  of  the  combustion  chamber  measures  25.7  cm.  The  chamber  is  mounted  together  with  a 
hydraulic  press  and  sealed  between  segments  with  standard  polymer  O-rings. 

The  injector  segment  is  separated  into  two  concentric  brass  sections.  The  forward  section  contains  a 
threaded  insen  which  forms  the  outer  diameter  of  the  fuel  annulus  and  may  be  replaced  between  tests  without 
complete  chamber  disassembly.  The  rear  section  provides  the  coupling  of  the  injector  to  the  propdlant  feed  system. 
In  addition,  fuel  plenum  pressure,  both  the  steady  and  oscillating  components,  as  well  as  fuel  injection  temperature 
are  measured  in  the  rear  section.  The  pressure  and  temperature  of  the  liquid  oxygen  are  measured  in  the  liquid 
supply  line  at  the  entrance  to  the  LOX  post.  The  baseline  dimensions  of  the  injector  were  modeled  after  the  mjeaor 
element  of  the  SSME  prebumer.  Additional  fuel  annuli  were  fabricated  for  the  first  series  of  tests  in  which  the 
effect  of  gas-to-liquid  velocity  ratio,  density  ratio  and  mixture  ratio  on  injector  performance  were  to  be  examined. 

Detail  of  the  window  segment  is  provided  in  Figure  2,  illustrating  the  60®  orientation  of  the  window 
principal  axes.  The  segment  is  intended  for  droplet  size  and  velocity  measurements  within  the  liquid  oxygen  spray 
using  a  phase  Doppler  interferometric  device.  To  reduce  or  eliminate  contaminants  from  collecting  on  the  windows 
during  engine  testing,  a  purge  gas  section  was  incorporated  into  the  design  of  the  window  segment.  Eight  small 
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holes  inject  nitrogen  gas  across  the  inner  surface  of  the  windows  within  the  window  recess  region.  Prelinunary  tests 
with  the  window  segment  located  at  the  chamber  midsection  verified  adequate  window  purging  with  nitrogen  flow 
rates  less  than  10%  of  the  total  propellant  flow  rate.  Measurements  of  the  chamber  inner  surface  temperature  and 
the  steady  and  oscillating  components  of  chamber  pressure  are  made  in  this  segment. 

The  nozzle  segment  houses  a  graphite  nozzle  throat  insen  that  is  sized  for  a  specific  propellant  mass  flow 
rate  and  chamber  pressure.  A  stainless  steel  flange  retains  the  insert  within  the  nozzle  segment  while  serving  a 
second  function  as  a  pressure  relief  mechanism.  During  instances  of  excessive  chamber  pressures,  the  nozzle  flange 
is  HwgtipH  to  break  away  from  its  fasteners  allowing  for  nozzle  removal  and  nqtid  chamber  depressurization. 
Propellant  ignition  is  accomplished  with  a  gaseous  oxygen/gaseous  hydrogen  torch  igniter.  Both  gases  are 
introduced  into  a  small  combustion  chamber  attached  to  the  side  of  the  igniter  segment  and  ignited  by  a  high  voltage, 
electrical  discharge.  The  resulting  flame  is  drawn  into  the  main  chamber  to  light  die  main  propellanc. 


Fig.  1  Single  element  rocket  combusuon  chamber. 


Fig.  2  Detail  of  window  cross-secuon. 
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ABSTRACT  NOMENCLATURE 


There  is  no  definitive  knowledge  of  which  of 
several  concurrent  processes  ultimately  results  in 
unstable  combustion  within  liquid  propellant  rocket 
motors  employing  shear  coaxM  injectors.  Possible 
explanations  are  a  detrimraital  change  in  the  atomization 
characteristics  due  to  a  decrease  in  the  velocity  ratio,  the 
disqtpearance  of  a  stabilizing  recirculation  region  at  the 
base  of  the  LOX  post  or  a  change  in  the  gas  side  injector 
pressure  drop  allowing  acoustic  coupling  to  the 
propellant  feed  systan.  The  spray  charactaistics  of  a 
shear  coaxial  injector  were  examined  under 
noncombusting  conditions  with  water  and  air  at 
atmospheric  pressure  and  with  liquid  and  gaseous 
nitrogen  at  elevated  pressures.  Two  injectors  were 
evaluated  with  different  LOX  post  designs,  nontapered 
and  tapered.  Initial  spray  visualizadon  tests  with  water 
and  air  revealed  two  modes  of  atomizadon,  namely  a 
resonant  and  nonresonant  mode,  and  an  inherent 
periodicity  in  the  clustering  of  droplets  within  the  spray. 
Subsequent  measurements  with  a  phase  Doppler 
interferometric  device  quantified  the  differences  in  spray 
characteiisdcs  between  the  two  modes.  Tests  at  elevated 
pressure  using  LN2/GN2  did  not  provide  any  indicadon 
of  the  resonant  mode.  Visualizadon  tests  did  reveal  a 
contracdon  of  the  spray  with  increasing  chamber 
pressures  below  the  cridcal  pressure.  At  pressures  above 
the  cridcal  point,  stroboscope  images  revealed  an 
elongadon  of  the  liquid  breakup  length.  Droplet  size  and 
velocity  measurements  were  made  for  three  oxidizer  to 
fuel  mixture  rados  and  a  range  of  chamber  pressures. 
The  results  indicated  improved  atomizadon  at  lower 
mixture  rados,  but  increased  droplet  sampling  difficulty. 
Velocity  measurements  in  the  recess  region  at  the  base 
of  the  LOX  post  were  obtained  with  water  and  air  at 
atmospheric  pressure  and  revealed  reverse  flow 
indicadve  of  a  recirculadon  region.  Finally,  the  results 
for  the  injector  response  from  a  linearized  lumped- 
element  model  are  presented  as  a  funedon  of  temperature 
and  frequency  and  compared  to  experimental  results. 


A  area  of  annulus 

df  fuel  annulus  inner  diameter 

d,  post  inner  diameter 

dg  post  outer  diameter 

C  capacitance  term  (Eqm.  8) 

g  gravitadonal  constant 

1  inertia  term  (Eqtn.  9) 

L  annulus  length 

m  mass  flow  rate 

N  response  factor 

Pc  chmber  pressure 

Pi  plenum  pressure 

Pi  annulus  entrance  pressure 

s  Laplace  operator 

T  temperature 

U  velocity 

V  plenum  volume 

W  mass  flow  rate  through  injector,  model  notadon 
Wj  mass  flow  rate:jth  mechanism,  model  notadon 
Wt  total  mass  flow  rate,  model  notadon 

p  density 

<7  liquid  surface  tension 

^  specific  heat  rado 

d  phase  angle  between  P'  and  W' 

m  frequency 

H  viscosity 

AP  Pd- Pi 

APg  Pd  -  Pc.  fuel  pressure  drop 

Subscripts 

g  gas  or  fuel 

H2  hydrogen 

I  hquid  or  oxidizer 

lox  liquid  oxygen 

max  peak  amplimde  of  sine  wave 

SiiperscriPLs 


_ _  -  average  value 
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INTRODUCTION 

Shear  coaxial  injectors  are  currently  found  in 
liquid  propellant  rocket  engines  which  use  hquid  oxygen 
and  hydrogen  as  the  propellant,  such  as  the  Space  Shuttle 
Main  Engine  (SSME)  and  the  first  stage  Ariane  5 
Vulcain  engine.  Their  use  in  such  bipropellant  systems 
enables  an  effective  means  of  liquid  oxidizer  atomization 
by  the  gaseous  fuel  and  subsequent  mixing  of  the 
propellant.  Liquid  oxygen  is  injected  into  the 
combustion  chamber  through  a  central  lube  while 
hydrogen  in  a  gaseous  fonn  is  fed  through  an  outside 
annular  passage  at  a  high  velocity  relative  to  the  liquid 
jet.  Impelled  by  liquid  turbulence  and  gas  to  liquid 
interfacial  interactions,  liquid  is  shipped  from  tlie  jet  imd 
entrained  into  the  surrounding  shear  How  where  the 
liquid  ligaments  and  droplets  experience  further  breakup 
imd  atomization.  The  resulting  spray  then  vaporizes  and 
mixes  with  the  gaseous  fuel  to  produce  a  volatile  mixture 
for  combustion.  Since  the  spatial  distribution  of  the 
droplets  downstrejun  of  the  injector  faceplate  establishes 
the  conditions  for  ilie  release  of  mass  and  energy  witliin 
the  chtunber,  the  emergence  of  unstable  pressure 
oscillations  within  the  chamber  is  strongly  coupled  to  the 
atomization,  vaporization  and  burning  characteristics  of 
the  injector  spray.  Such  combustion  instabilities  cire  an 
inherent  phenomenon  of  liquid  propellant  rocket 
engines,  in  which  strong  exothermic  reactions  occur 
within  a  small  volume,  high  pressure  environment. 
Small  fluctuations  at  the  onset  of  or  during  steady  state 
operation  may  initiate  resonant  pressure  oscillations  that 
grow  rapidly  in  amplitude  and  eventually  lead  to  the 
fiiilurc  of  the  engine. 

Extensive  engine  testing  during  the  I960's 
discovered  a  strong  dependency  of  engine  stability  on 
the  fuel  temperature  at  injection.  This  observation  lead 
to  a  stability  rating  technique,  in  which  the  fuel 
temperature  was  ramped  to  decreasing  values  during 
steady  stale  engine  operation.  The  temperature  at  which 
the  engine  transitioned  into  an  unstable  mode  of 
operation  became  a  measure  of  the  stability  of  the 
engine,  lowest  temperatures  meaning  highest  engine 
stability.  The  parametric  tests  on  a  full  scale,  20,000  lb. 
engine  at  NASA  Lewis  isolated  the  effects  of  several 
injector  design  and  operating  parameters  on  the 
hydrogen  temperature  stability  limit.  The  most 
significant  appeared  to  be  the  relative  velocity  between 
the  fuel  and  oxidizer,  a  presumption  supported  by  the 
argument  that  a  large  momentum  exchange  between  the 
two  fluids  is  necessary  for  the  disruption  of  the  liquid  jet, 
Wanhainen  et  al.^  showed  that  unstable  combustion 
iu^ose  for  a  variety  of  injector  mass  flow  rates  and  exit 
areas  when  the  ratio  of  fuel  to  oxidizer  velocities 
dropped  below  6.5.  In  terms  of  engine  design 
piiraineters,  the  LOX  post  recess  depth  and  fuel  injection 
angle  relative  to  the  LOX  post  both  enhanced  the 
stability  of  tiie  engine  by  lowering  liic  fuel  temperature 


stability  limit.  Ultimately,  these  parameters  were 
incorporated  into  a  single  correlation,  which  accounted 
for  ail  variations  in  injector  design  and  operation 
investigated  in  this  work^.  Based  on  the  hydrogen 
threshold  temperature,  the  correlation  is  stated  as 
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which  illustrates  that  for  a  given  engine  chamber 
pressure  ,  and  mixture  ratio,  ,  the  fuel  side 

pressure  drop,  ,  governs  the  stability  characteristics 

of  the  engine.  This  result  would  suggest  that  instability 
is  driven  by  an  acoustic  coupling  between  the  injector 
hydrodynamics  and  spray  formation  and  combustion 
processes  within  the  chamber  since  higher  pressure 
drops  can  be  effected  by  modifications  to  the  fuel 
annulus  inlet  configuration  in  addition  to  higher  gas  exit 
velocities.  The  contribution  of  liquid  jet  diameter  to  the 
hydrogen  threshold  temperature  is  less  apparent  due  to 
the  uncertainty  of  the  effect  of  larger  jet  diameter  on 
liquid  jet  atomization.  By  increasing  the  jet  diameter  for 
a  fixed  liquid  mass  flow  rate,  tlie  liquid  will  exit  the 
injector  at  a  lower  velocity  and  likewise  produce  a  larger 
relative  velocity  that  may  improve  atomization.  On  the 
other  hand,  the  resultant  increase  in  the  interfacial 
surface  area  may  reduce  the  momentum  exchange  per 
unit  area  and  consequently  worsen  the  atomization  of  the 
jet 

Another  proposed  explanation  for  the  emergence 
of  unstable  operation  in  liquid  propellant  rocket  engines 
is  that  a  recirculation  region  acting  as  a  flameholder 
exists  downstream  of  the  LOX  post  tip.  Below  a 
minimum  relative  velocity  between  the  liquid  oxygen 
and  gaseous  hydrogen,  the  recirculation  region  becomes 
too  weak  to  act  as  a  flameholder  and  the  combustion 
zone  moves  away  from  the  injector  face  to  a  location 
where  it  can  interact  more  strongly  with  the  chamber 
acoustics  modes.  Liang  and  Schumann  have  examined 
this  idea  with  an  experimental  and  computational 
investigation  of  gaseous  oxygen  and  hydrogen  coaxial 
injectors^.  They  examined  several  injectors  designed  to 
produce  recirculation  regions  of  dif  ferent  sizes  but  found 
that  all  injectors  tested  showed  the  combustion  region 
anchored  to  the  base  of  the  injector.  Cox  made  LDV 
measurements  of  the  gas  velocity  in  the  fuel  annulus  in 
an  SSME  injector  element  and  found  flow  separation  at 
the  start  of  the  LOX  post  taper^. 

In  general,  the  underlying  physical  processes 
influencing  combustion  instability  ine  not  completely 
understood  due  to  a  deficiency  in  experimental  data  on 
the  characteristics  of  the  injector  spray  and  in  computer 
codes  capable  of  analyzing  the  complex,  unsteady,  three 
dimensional  processes  within  the  engine.  With  the 


present  developmental  status  of  high  performance 
computer  architectures  and  advanced  numerical 
algorithms,  sophisticated  computer  models  are  being 
generated  to  provide  performance  and  subility 
predictions  in  the  design  of  future  liquid  propellant 
rocket  engines^.  In  order  to  validate  these  models, 
accurate  experimental  data  is  required  on  the 
characteristics  of  the  injector  flowfield^.  This 
information  includes  the  spatial  and  time  varying 
distribution  of  droplet  size  and  velocity,  turbulence  scale 
and  intensity,  structure  of  the  chemical  reaction  zones, 
and  the  distribution  of  temperature  and  concentration  of 
the  reactants  and  products  throughout  the  flow  field. 

Advanced  laser  based  diagnostic  techniques 
have  been  devised  concurrent  to  the  improved 
development  of  computer  codes  and  are  capable  of 
obtaining  the  required  information  on  the  spray.  As  a 
consequence,  the  impetus  in  academia  and  industry  is  to 
exploit  the  wide  variety  of  available  laser  based 
diagnostic  techniques  for  the  characterization  of  the 
complete  injector  atomization,  vaporization  and 
combustion  processes  under  conditions  similar  to  actual 
engine  operation.  As  a  precursor  to  combusting 
conditions,  many  cold  flow  experiments  have  been 
performed  with  a  single  shear  coaxial  injector  using 
water  as  the  liquid  oxygen  simulant  in  order  to 
investigate  the  mechanisms  of  liquid  jet  atomization 
without  the  influence  of  combustion^*^°.  Common  fuel 
simulants  have  been  air,  nitrogen  and  helium.  Since 
liquid  nitrogen  offers  better  similitude  with  liquid 
oxygen  in  terms  of  its  low  boiling  temperature  and  the 
important  physical  properties  which  influence 
atomization,  namely  density,  viscosity  and  surface 
tension,  injectors  have  also  been  tested  with  LN2  and 
either  gaseous  nitrogen  or  gaseous  helium  as  the  tuel 
simulant* 

At  the  Propulsion  Engineering  Research  Center 
of  The  Pennsylvania  State  University,  multiple  laser 
based  diagnostic  techniques  are  being  used  in  the 
characterization  of  a  unielement  shear  coaxial  injector 
under  combusting  and  noncombusting  conditions  and 
elevated  pressures*^  *^.  The  focus  of  the  present 
research  effort  was  to  characterize  the  evolutionary 
behavior  of  an  atomizing  and  evaporating  spray 
produced  by  a  single  shear  coaxial  injector  operating 
under  noncombusting  conditions.  To  this  end  a 
laboratory  injector  was  designed  and  fabricated  for 
operation  at  atmospheric  and  elevated  pressures.  Water 
and  liquid  nitrogen  were  selected  to  serve  as  simulants  to 
liquid  oxygen  while  air  and  gaseous  nitrogen  supplanted 
gaseous  hydrogen  as  the  fuel.  Preliminary  investigations 
began  with  the  combination  of  water  and  air  to  examine 
injector  atomization  characteristics  without  the  influence 
of  liquid  vaporization.  Through  these  tests  general 
trends  in  atomization,  both  qualitative  and  quantitative, 
were  to  be  observed  for  comparison  to  the  subsequent 


high  pressure  investigations  with  liquid  and  gaseous 
nitrogen.  Furthermore,  investigations  of  the 
recirculation  region  with  water  and  air  provide  an 
evaluation  of  the  feasibility  of  this  technique  for  future 
investigations  under  reacting  conditions.  Nonintrusive 
diagnostic  techniques  were  implemented  to  obtain 
information  on  the  global  behavior  of  the  spray  and  to 
quantify  the  important  spray  parameters,  namely  droplet 
mean  size  and  mean  axial  velocity,  at  several  locations 
within  the  spray.  The  techniques  included  planar  laser 
imaging  and  stroboscopic  photography  to  visualize  the 
morphology  of  the  spray  and  phase  Doppler 
interferometry  to  measure  the  droplet  parameters.  In 
order  to  detect  potential  coupling  between  the  chamber 
and  injector  acoustics,  a  simplifled  stability  model  was 
constructed  for  the  flow  dynamics  of  several  processes  in 
a  cold  flow  pressure  chamber.  In  support  of  this  effort, 
the  chamber  and  injector  pressure  oscillations  were 
measured  to  evaluate  the  analytical  results  of  the  model. 

RFSEARPH  APPROACH 

On  the  Simulation  of  LOX/GH2 

To  observe  the  injector  spray  characteristics 
without  the  influence  of  combustion,  cold  flow  tests  are 
often  performed  with  alternative  liquids  and  gases  that 
offer  a  high  degree  of  similitude  with  the  actual 
propellants  (LOX/GH2)  with  respect  to  the  important 
physical  properties  which  influence  the  spray  formation 
processes,  namely  liquid  viscosity  and  surface  tension 
and  gas  density.  In  the  present  research  effort  water  and 
liquid  nitrogen  were  used  as  simulants  to  liquid  oxygen 
while  air  and  gaseous  nitrogen  served  as  the  simulants  to 
hydrogen.  The  relevant  physical  properties  of  these 
simulants  are  compared  to  those  of  LOX/GH2  in  Table 
1.  Pressure  and  temperature,  at  which  the  properties 
were  evaluated,  are  included  in  the  table  and  correspond 
to  the  operating  conditions  of  the  SSME  prebumer  in  the 
case  of  LOX/GH2  and  to  the  expected  baseline  run 
conditions  in  the  case  of  the  simulants.  These 
conditions  are  listed  in  Table  2. 

Based  on  the  large  differences  in  liquid 
viscosity  and  surface  tension  between  liquid  oxygen  and 
water,  water  is  not  expected  to  yield  spray  characteristics 
similar  to  liquid  oxygen.  Furthermore,  water  is  not 
subject  to  evaporation  processes  at  the  injection 
temperature  of  the  air.  Nevertheless,  its  availability  and 
ease  of  use  make  water  a  preferred  liquid  for  the  study  of 
the  fundamental  processes  of  liquid  jet  breakup  and 
atomization.  For  this  reason  water  has  been  selected  for 
use  with  air  in  the  initial  spray  characterization  tests.  At 
the  baseline  pressure  and  temperature  specified  in  Table 
1  liquid  nitrogen  offers  physical  properties  closer  to 
those  of  liquid  oxygen,  especially  in  the  case  of  liquid 
viscosity  and  surface  tension.  Since  the  boiling 
temperature  of  LN2  2.5  MPa  is  significantly  less  than 
the  coflowing,  nitrogen  gas  temperature,  liquid  nitrogen 
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will  undergo  evaporation,  thereby  providing  more 
insight  into  evaporation  etfects  on  spray  evolution.  The 
lower  density  of  LN2  may  not  introduce  spray 
behavioral  differences  when  considering  the  ratio  of 
liquid  to  gas  density  to  be  a  more  influential  parameter  in 
liquid  atomization  than  liquid  density  alone.  For 
LN2/GN2  the  density  ratio  is  20.1  while  that  for 
LOX/GH2  is  26.5. 

In  Table  2  the  conditions  for  LN2/GN2  were 
selected  to  approximate  the  injection  velocities  of  LOX 
and  H2.  Although  the  sound  speed  for  gaseous  nitrogen 
limited  the  maximum  exit  velocity  to  approximately  300 
m/s,  the  resulting  mass  flow  rate  of  gaseous  nitrogen 
provided  a  mixture  ratio  almost  equal  to  that  of  the 
prebunier.  The  conditions  for  water/air  were  chosen  to 
produce  a  similar  mixture  ratio  and  relative  velocity 
while  insuring  proper  atomization  of  the  water  jet. 
Overall,  a  large  degree  of  similitude  is  achieved  between 
LOX/GH2  and  LN2/GN2.  Most  nondimensional 
parameters,  including  mixture  ratio  and  velocity  ratio, 
iu-e  of  the  same  order  of  magnitude  between  tliese  two 
fluid  combinations  with  the  only  exception  being  the  gas 
Reynolds  number.  According  to  Farago  and  Chigier^^, 
die  brcjikup  of  the  liquid  jet  can  be  correlated  by  the 
liquid  Reynolds  number  and  aerodyniunic  Weber 
number.  The  magnitude  of  these  parameters  are  very 
close  for  LN2/GN2  and  LOX/GH2  suggesting  that  both 
combinations  would  possess  similar  breakup 
characteristics.  The  water/air  combination  differs 
significantly  from  LOX/GH2  in  terms  of  the 
nondimensional  parameters  and  consequently  is 
expected  to  yield  different  characteristics  in  atomization. 

Experimental  Set-Uo. 

The  baseline  dimensions  of  the  shear  coaxial 
injector  are  shown  in  Fig.  1  cUid  were  derived  from  the 
dimensions  of  tlie  prototype  injector  element  of  ihe  fuel 
prebumer  to  the  Space  Shuttle's  main  engine.  A  second 
LOX  post  was  fabricated  with  a  T  taper  at  the  tip  to 
match  tlie  shape  of  the  production  injector  element. 
Pressure,  both  ^e  mean  and  oscillating  components,  and 
temperature  were  measured  in  the  fuel  plenum  region 
upstream  of  the  fuel  annulus.  The  mean  pressure  and 
suitic  temperature  of  the  liquid  simulant  were  measured 
in  the  liquid  supply  line  at  the  entrance  to  the  LOX  post. 

For  the  elevated  pressures  required  of  LN2 
testing,  a  pressure  chamber  was  constructed  to  operate  at 
pressures  up  to  lO.O  MPa  and  to  provide  optical  access 
for  laser  ba.sed  diagnostic  techniques.  The  chamber, 
which  is  depicted  in  Fig.  2,  is  22.9  cm  in  length  and  has 
a  square  cross-section  with  sides  measuring  10.2  cm. 
Optical  access  is  through  2.54  cm  thick  plexiglass  plates 
mounted  to  both  sides  of  the  chamber  mid  reinforced 
will)  high  strengtii,  carbon  steel  flanges,  1.27  cm  thick. 
Slits  measuring  3.2  cm  wide  by  15.2  cm  long  were  cut 


into  the  flanges  to  maintain  optical  accessibility  to  the 
chamber.  The  position  of  these  slits  depended  on  the 
diagnostic  technique  in  use,  as  is  illustrate  in  Fig.  2  for 
the  case  of  phase  Doppler  interferometry.  Four  ports 
drilled  into  the  injector  end  of  the  chamber  serve  to 
provide  additional  gas  for  chamber  pressurization  and  to 
suppress  recirculating  flow.  By  the  replacement  of  the 
nozzle  plate  a  range  of  chamber  pressures  can  be 
investigated  for  specific  propellant  mass  flow  rates. 
Mean  chamber  pressure  and  internal  gas  temperature 
were  measured  at  the  center  of  the  chamber  as  illustrated 
in  Fig.  2.  High  frequency  pressure  oscillations  were 
measured  at  the  exit  plane  of  the  injector.  To  facilitate 
particle  sizing  and  velocity  measurement  at  several  axial 
and  radial  locations  of  the  spray,  the  entire  chamber  was 
mounted  on  a  2-D  translation  platform  and  connected  by 
flexible  lines  to  the  propellant  feed  system. 

Test  Facility 

The  initial  spray  characterization  tests  were 
performed  at  atmospheric  pressure  using  water  and  air. 
Water  was  lapped  from  the  local  area  network  and 
metered  into  a  graduated  cylinder  to  determine  the  liquid 
mass  flow  rate.  Air  was  supplied  by  compressors  within 
the  Center  and  regulated  by  means  of  a  standard  pressure 
gauge  rated  at  100  psig.  Knowing  the  supply  pressure 
and  temperature  of  air  within  the  plenum  region  of  the 
injector,  the  air  mass  flow  rale  was  calculated  from  the 
isentropic  relations  for  a  truncated  nozzle,  corrected  for 
two-dimensional  and  turbulence  effects. 

The  cold  flow  experiments  with  liquid  nitrogen 
were  performed  at  the  Cryogenic  Combustion 
Laboratory  of  The  Pennsylvania  State  University.  The 
laboratory  was  designed  for  the  experimentation  of 
subscale  rocket  engines  and  can  accommodate  liquid 
oxygen  and  gaseous  hydrogen  and  other  propellant 
combinations  under  reacting  or  cold  flow  conditions. 
Facility  performance  is  rated  at  0.45  kg/s  of  liquid 
oxygen  and  0.11  kg/s  of  gaseous  hydrogen  for  chamber 
pressures  up  to  approximately  8.0  MPa.  The  supply  of 
gaseous  nitrogen  to  the  injector  and  the  four  purge  ports 
was  from  the  primary  and  secondiiry  supply  lines  of  tlie 
gaseous  hydrogen  feed  subsystem  of  the  laboratory. 
Two  subsystems  comprise  the  overall  gas  feed  system 
and  are  supplied  by  clusters  of  pressurized  gas  cylinders. 
Liquid  nitrogen  was  supplied  to  the  test  section  from  a 
114  liter  vacuum  jacketed  storage  tank  by  a  gaseous 
nitrogen  dome  loading  system.  To  reduce  the  time 
required  for  cooling  down  the  cryogen  supply  lines,  a 
liquid  nitrogen  cooling  jacket  was  installed  onto  the  lines 
between  the  storage  tank  and  the  test  section.  Flow 
metering  for  both  gaseous  and  liquid  propellants  was 
accomplished  with  critical  flow  venturies  fitted  with 
pressure  transducers  and  thermocouples  upstream  and 
downstream  of  the  venturi.  Since  cavitation  is  a 
significant  concern  in  LN2  atomization  experiments  at 


Baseline  Dimensions 


Rgure  1.  Shear  coaxial  injector  schematic  with  baseline  dimensions. 


Table  1  Comparison  of  propellant/simulant  physical  properties. 


Table  2  Comparison  of  injector  operating  conditions. 


Operating  Conditions 
(units) 

Symbol/ 

Formula 

SSME 

Prebumer 

LN2/N2 

Simulants 

H20/Air 

Simulants 

Chamber  Pressure  (MPa) 

Pc 

34.1 

2.5 

0.1 

Fuel  Pressure  Drop  (MPa) 

2.95 

U6 

0.062 

Fuel  mass  flow  (kg/s) 

0.138 

0.077 

0.0033 

Ox.  mass  flow  (kg/s) 

rhi 

0.135 

0.071 

0.0031 

Fuel  exit  velocity  (m/s) 

Us 

360 

300 

260 

Ox.  exit  velocity  (m/s) 

Ui 

30 

30 

0.8 

Relative  velocity  (m/s) 

U,-U, 

330 

270 

259 

Velocity  ratio 

UsIUi 

12 

10 

325 

Density  ratio 

PilPs 

26.5 

20.2 

726 

Mixture  ratio 

0.98 

0.92 

0.94 

Momentum  ratio 

m.UJrhiUi 

12.2 

10.8 

347 

Gas  Re  No  (xlO^) 

P.uX‘if-do)IP^ 

3.37 

0.62 

0.028 

Liquid  Re  No  (xlO^) 

PiU,dilHi 

6.60 

8.56 

0.013 

Weber  No.  (xlO^) 

Ps(Us-Ui)d-^h- 

14.12 

34.61 

0.026 

Ohnesorge  No.  (xlO"^^ 

8.42 

10.84 

33.62 
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subcritical  conditions,  a  criterion  was  established  based 
on  the  work  by  Nurick*^  to  determine  whether  a  test  was 
influenced  by  LN2  cavitation.  The  criterion  required 
that  the  difference  between  the  LOX  post  upstrciun 
pressure  :md  vapor  pressure  must  exceed  die  pressure 
drop  through  the  LOX  post  by  a  factor  determined  from 
the  entrance  geometry  and  discharge  coefficient  of  the 
LOX  post.  Although  actual  LOX/GH2  rocket  engines 
operate  in  the  supercritical  regime,  the  dynamics  of 
cavitation  at  subcritical  pressures  could  yield  spray 
structures,  drop  size  distributions,  and  overall  propellant 
mixing  quality  not  concurrent  with  ideal  injector 
operation. 

Spray  Visualizadon  Technique. 

The  global  characteristics  of  the  shear  coaxial 
injector  .spray  were  examined  by  imaging  the  scattered 
light  from  a  laser  sheet  passing  through  the  axis  of  the 
spray  tind.  in  a  separate  test  series,  from  a  stroboscope 
illuminating  the  backside  of  the  spray.  With  this 
approach  the  momentary  structure  of  the  spray,  including 
the  liquid  core  and  large,  detached  ligaments,  could  be 
visualized  for  evaluation  of  the  effects  of  injector  design, 
propellant  flow  rates  and  ambient  conditions. 
Furthermore,  dilute  droplet  flow  regions  could  be 
identified  within  die  spray  for  subsequent  droplet  sizing 
and  velocity  measurement  by  means  of  phase  Doppler 
interferometry.  For  the  laser  sheet  imaging  tests  a 
frequency  doubled  Nd:YAG  laser  (532  nm),  pulsed  at  10 
Hz  with  a  pulse  duration  of  10  ns.  was  used  with  a 
spherical  and  cylindrical  lens  combination  to  lorm  a 
laser  sheet  with  a  thickness  of  approximately  0.15  imn 
;uid  negligible  divergence  at  the  injector  centerline.  The 
scattered  light  from  the  spray  was  vi.sualized 
perpendicuhtf  to  die  sheet  with  a  35  mm  Nikon  camera. 
A  1  nm  bandpass  tiller  centered  at  532  nm  was  insudled 


in  front  of  the  ciunera  lens  to  eliminate  ambient  light. 
This  technique  as  well  as  the  stroboscope  technique  were 
employed  during  the  liquid  nihogen  tests  at  elevated 
pressures.  For  the  water/air  tests  only  the  stroboscope 
was  used  to  illuminate  the  spray.  The  orientation  of  the 
stroboscope  for  both  tests  was  toward  the  backside  of  the 
spray  with  respect  to  the  camera  and  in  the  direction  of 
the  injector  face.  Preliminary  flow  visualization 
experiments  with  water  and  air  at  atmospheric  pressure 
revealed  that  both  the  dispersed  droplet  region  and  the 
structure  of  the  liquid  core  may  be  imaged  at  this 
orientadon  with  the  35  mm  camera.  The  flash  rate  of  the 
strobe  was  set  at  10  Hz  with  a  flash  duration  of  3  jls. 

Phase  Doppler  Interferometry. 

Quantitative  information  on  the  spray  and 
recirculation  region  was  pursued  with  a  phase  Doppler 
interferometric  (PDI)  device,  developed  by  Aerometrics, 
which  enabled  simultaneous  measurement  of  the 
temporally  averaged  droplet  size  and  velocity 
distribution  at  a  point  within  the  spray  1^.  The  optical 
configuration  for  these  tests  included  a  CW  argon  ion 
laser  emitting  at  514.5  nm,  a  Bragg  cell  based 
transmitter,  ;uid  a  receiver  oriented  at  30°  from  the 
forward  direction  of  the  transmitter  optical  axis.  The 
transmitter  was  fitted  with  a  500  mm  f.l.  lens  to  torm  a 
focused  beam  waist  of  234  ^tm.  Using  a  1(X)  |lm  x  1 
mm  .spatial  filter  and  a  lens  pair  consisting  of  a  495  mm 
f.l.  collection  lens  and  a  250  mm  f.l.  focusing  lens  within 
the  receiver,  the  effective  probe  area  became 
approximately  212  fim  x  234  jim.  Doppler  signal 
filtering,  acquisition  and  analysis  were  performed  by 
dedicated  processing  electronics  and  software  developed 
by  Aerometrics.  During  the  LN2/GN2  tests  data 
sampling  was  over  a  2.5  second  interval  triggered  by  a 


programmable  laboratory  controller  at  2.5  seconds  into 
each  test,  at  which  point  the  pressure  within  the  chamber 
bad  reached  a  steady  value. 

Analytical  Modeling  of  Injector  Response 

A  simple  model  based  on  the  Rayleigh  criterion 
was  constructed  for  the  flow  dynamics  of  two  processes 
in  the  cold  flow  liquid  rocket  chamber  to  determine  the 
effect  of  various  design  and  operating  parameters  on 
combustion  stability.  The  mechanisms  included  are  the 
fuel  flow  through  the  injector,  and  the  flow  of  exhaust 
gases  through  the  nozzle. 

The  Rayleigh  criterion  states  that  a  wave  will 
grow  if  heat  or  mass  is  added  in  phase  with  pressure. 
Conversely,  the  wave  will  damp  if  the  addition  is  out  of 
phase.  If  several  mechanisms  are  contributing  to  the 
heat  or  mass  release,  then  the  wave  growth  is  determined 
by  the  net  in-phase  or  out-of-phase  heat  or  mass 
addition.  By  assuming  that  both  the  heat  or  mass 
addition  and  the  pressure  oscillations  can  be  represented 
by  a  perturbation  added  to  their  steady  state  values, 
linear  models  of  the  dynamic  response  of  the  processes 
may  be  obtained. 

A  response  factor  is  defined  as  the  real,  or  in- 
phase,  part  of  the  mass  flow  perturbation  with  response 
to  the  chamber  pressure  perturbation.  The  response 
factor  for  a  transverse  mode  was  derived  in  Reference 
18: 

Al= - cose  (2) 

c,  max 


The  response  factor  can  be  either  positive  or  negative 
(tfpfnrting  on  the  phase  angle  between  the  pressure  and 
flow  perturbations.  The  analytical  derivation  of  the 
response  factor  for  the  various  flow  processes  is  detailed 
in  the  Appendix. 

The  stability  limit  for  a  transverse  mode  of 
oscillation  within  the  chamber  was  derived  in  Reference 
19  from  a  mass  balance  of  an  axial  column  of  gas: 
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This  relation  thus  represents  the  total  response  as  a  linear 
combination  of  each  response  factor  weighted  by  the 
fraction  of  mass  flow  contributed  by  that  process.  A 
positive  response  corresponds  to  an  unstable  system, 
while  a  negative  response  indicates  stability. 

water; AIR  RESULTS 


Spray  Visualization 

The  majority  of  the  spray  visualization  tests 
were  performed  with  the  shear  coaxial  injector  using  the 


nontapered  LOX  post  recessed  to  2.54  mm  inside  the 
fuel  annulus.  The  water  flow  rate  was  varied  between 
3.0  g/s  and  54.0  g/s  with  air  supply  pressures  up  to  184 
kPa.  Higher  supply  pressures  were  not  possible  due  to 
flow  choking  considerations  within  the  fuel  annulus.  To 
examine  the  effect  of  injector  design  on  spray 
morphology,  additional  tests  were  conducted  without  the 
post  recessed  and  with  the  tapered  LOX  post  recessed  to 
2.54  mm. 

Aside  from  the  expected  improvement  in 
atomization  with  decreasing  water  flow  rate  or 
increasing  air  flow  rate,  several  interesting  phenomena 
were  observed  during  these  tests.  At  all  water  mass  flow 
rates  the  spray  produced  regions  of  high  and  low  droplet 
number  density  occurring  at  wave-like  frequency.  These 
droplet  clusters  appeared  most  distinctly  at  the  edge  of 
the  spray  and  at  the  highest  mass  flow  rate  of  air.  The 
wavelength,  or  distance  between  adjacent  clusters,  was 
shortest  at  this  flow  rate  and  increased  with  decreasing 
air  mass  flow  rate.  Furthermore,  it  appeared  that  the 
bandwidth  of  these  waves  widened  at  lowering  air  mass 
flow  rates.  This  behavior  is  depicted  in  Fig.  3  for  6.0  g/s 
of  water  and  3.7  g/s  of  air  using  the  nontapered  LOX 
post  without  a  recess.  The  figure  illustrates  the  inherent 
periodicity  of  liquid  jet  disruption  and  droplet  flow 
within  the  spray.  Faragd  and  Chigia:  have  identified  this 
behavior  as  a  normal  mode  of  jet  disintegration  from 
their  experiments  with  a  twin  fluid  atomizer^  The 
occurrence  of  periodic,  coherent  structures  within  a  shear 
coaxial  injector  spray  supports  their  findings  while 
revealing  an  unsteadiness  of  the  spray  that  can  affect  the 
evaporation  and  combustion  characteristics  of  the  spray 
and  likewise  the  overall  stability  characteristics  of  the 
engine. 


During  these  tests  an  unusual  whistling  noise, 
characteristic  of  a  resonating  pipe,  emanated  from  the 
injector  when  the  post  was  recessed  inside  the  fuel 
anniiliis-  This  behavior  arose  with  both  post  shapes  but 
was  less  pronounced  with  the  nontapered  post.  No 
noticeable  change  in  the  structure  of  the  spray  was 
observed  with  the  nontapered  post  upon  the  emergence 
of  the  whisding  phenomenon.  Furthermore,  the 
amplitude  of  the  noise  was  relatively  low,  which  in 
conjunction  with  the  steady  spray  behavior  would 
suggest  that  the  nont^iered  post  operated  close  to  a  point 
of  resonance.  Tests  performed  without  a  recessed  post 
produced  no  whistling  noise  over  the  entire  range  of  air 
supply  pressures  and  water  mass  flow  rates. 

Operation  with  the  tapered  LOX  post  at  the 
nominal  recess  depth  of  2.54  mm  revealed  the  most 
noticeable  change  in  the  behavior  of  the  spray.  With 
increasing  air  supply  pressure,  the  injector  transitioned 
from  a  nonresonant  mode  of  operation  to  a  resonant 
mode  at  a  specific  air  pressure.  At  the  point  of 
resonance  the  injector  spray  shifted  from  a  solid  cone  to 
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a  hollow  cone  with  a  significant  increase  in  the  spray's 
angle  of  divergence.  The  liquid  jet  experienced 
enhanced  atomization  as  the  breakup  length  approached 
the  injector  exit.  Furthermore,  the  amplitude  of  the 
whistling  noise  became  more  audible  with  this  post 
shape  and  increased  with  increasing  air  mass  flow  rate. 
Fig.  3b  depicts  the  tapered  post  operating  at  resonance 
with  6.0  g/s  of  water  and  3.7  g/s  of  air.  Resonance  was 
detected  at  126  kPa  of  air  for  this  water  flow  rate.  As 
seen  from  Figs.  3a  and  3b,  the  onset  of  resonant 
atomization  has  resulted  in  a  marked  change  in  the 
overall  behavior  of  the  spray.  The  consequences  of  this 
behavior  are  profound  and  encompass  a  possible  means 
of  improved  injector  performance  or  a  mechanism  for 
injector  destabilization. 

Droplet  Size  and  Velocity  Measurements 

Since  a  complete  characterization  of  the  injector 
spray  field  has  been  conducted  by  Zaller  and  Klem^  at 
several  water  and  air  mass  flow  rates,  the  emphasis  of 
this  effort  was  to  obtain  a  database  of  spray 
characteristics  for  comparison  to  the  subsequent  elevated 
pressure  tests  with  liquid  nitrogen  and  to  quantify  the 
observed  differences  between  the  two  modes  of 
atomization.  The  injector  was  operated  at  6.0  g/s  of 
water  and  3.7  g/s  of  air  using  botli  the  tapered  and  non- 
tapered  LOX  post  reces.scd  inside  the  minulus.  Based  on 
tile  visualization  tests,  the  nontapered  post  is  not 
expected  to  transition  into  a  resonating  mode  of 
atomization  as  the  tapered  LOX  post  under  similar 


conditions.  Data  on  droplet  mean  velocity  and  mean  size 
were  obtained  as  a  function  of  radial  position  at  5  cm 
downstream  of  the  injector  face.  The  spray  data  are 
plotted  in  the  following  four  figures  and  depict 
noticeable  differences  between  the  two  sprays.  In  Fig.  4 
the  data  on  droplet  mean  axial  velocity  illustrate  a 
steepening  of  the  velocity  profile  when  transitioning  into 
the  resonating  condition  with  the  tapered  post  from  the 
non  resonating  condition  with  the  nontapered  post. 
Interestingly,  the  velocities  on  the  injector  centerline  are 
approximately  equal  for  both  configurations  but  depart 
outward  in  the  spray.  The  more  rapid  drop  in  axial 
velocity  for  the  tapered  post  would  suggest  that  the 
droplets  are  U*aveling  faster  in  the  radial  direction  than  in 
the  axial  direction. 

Examination  of  the  resultant  mean  drop  size 
distribution  in  Fig.  5  reveals  a  more  pronounced 
departure  from  the  normal  behavior  of  the  spray. 
According  to  the  results  of  Ref.  9,  the  Sauter  mean 
diameter  is  expected  to  decrease  at  increasing  radial 
distances  from  the  injector  centerline,  a  trend  which  is 
observed  in  Fig.  5  for  the  nontapered  post.  With  the 
tapered  post  smaller  Sauter  mean  diameters  are  found  on 
the  centerline  and  increase  radially  outward  in  the  spray. 
Overall,  the  tapered  LOX  post  produced  smaller  droplets 
within  the  region  of  the  spray  probed  with  the  PDPA.  At 
the  centerline  the  difference  in  sizes  is  most  substantial 
and  decreases  when  moving  out  in  the  spray. 


(a)  Nontapered  post  w/o  recess 


(b)  Tapered  post  w/recess 


Fig.  3  Photographic  results  of  aUnosphcric  vvatcr/air  spray  at  6.0  g/s  of  water  :uid  3.7  g/s  of  air. 


Hg.  4  Radial  variation  in  droplet  mean  axial  velocity 
for  two  LOX  post  shapes. 


Fig.  5  Radial  variation  in  droplet  Sauter  mean 
diameter  for  two  LOX  post  shapes. 


Fig.  6  Radial  variation  in  droplet  number  density  for 
two  LOX  post  shapes. 


Fig.  7  Radial  variation  in  droplet  volume  flux  for  two 
LOX  post  shapes. 


The  distribution  of  droplet  number  density  is 
presented  in  Fig.  6  for  both  injector  configurations.  The 
distribution  for  the  tapered  post  corroborates  the  hollow 
cone  behavior  observed  in  the  spray  visualization  tests. 
That  is,  the  majority  of  drops  reside  at  the  edge  of  the 
spray  with  a  decreasing  number  toward  the  center.  The 
volume  flux  of  droplets  at  each  point  of  measurement  is 
plotted  in  Fig.  7  for  both  injector  configurations.  In  the 
case  of  the  tapered  post  the  distribution  of  volume  flux  is 
rather  uniform  throughout  the  spray  as  compared 
to  the  steep  distribution  for  the  nontapered  post.  The 
occurrence  of  this  distribution  is  due  to  the  opposite 
trends  between  the  mean  velocity  and  the  number 
density  and  drop  size.  That  is,  the  steep  velocity  profile 
compensates  for  the  hollow  cone  distribution  in  drop  size 
and  number  density  to  create  the  uniform  volume  flux 
distribution.  Based  on  this  disfiibution  the  burning 


efficiency  of  the  tapered  post  under  resonating 
conditions  may  be  expected  to  be  higher  than  that  of  the 
nontapered  post.  But  if  the  range  over  which  the 
resonance  condition  occurs  is  narrow,  this  condition  may 
result  in  a  mechanism  of  combustion  instability  when  the 
injector  transitions  intermittently  between  the  two  modes 
as  a  result  of  pressure  fluctuations  within  the  injector  and 
combustion  chamber. 

Recirculation  Region 

The  (siginai  injector  design  with  the  straight  LOX 
post  was  modified  in  order  to  allow  LDV  access  into  the 
recess  region.  This  was  accomplished  by  inserting  a 
Pyrex  tube  concentrically  around  the  LOX  post,  such 
that  the  fuel  annulus  was  defined  by  a  transparent 
boundary  1.02  mm  thick.  The  inner  diameter  of  the 
Pyrex  fuel  annulus  was  5.00  mm.  All  measurements 
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were  Uiken  using  Aerometrics’  Phase  Doppler  Particle 
Analyzer  (PDPA)  with  the  droplet  diameter 
measurement  function  disabled.  The  system 
incorporates  a  40  MHz  frequency  shift  of  one  beam 
through  the  use  for  a  Bragg  Cell  so  tliat  tlie  direction  of 
die  pitrticlc  may  be  obtained  along  with  the  velocity 
measurement.  Both  the  transmitting  and  receiving  optics 
were  aligned  at  angles  of  15  degrees  from  the  horizontal. 

With  the  optics  in  this  configuration  and  the 
presence  of  the  curved  Pyrex  surtace,  the  refraction  of 
the  beams  as  well  as  any  probe  volume  effects  must  be 
considered.  In  this  case,  the  cylinder  is  larger  than  the 
diameter  of  the  laser  beam.  Therefore,  the  probe  volume 
will  not  undergo  any  appreciable  distortion.  Refraction 
of  the  beam,  however,  does  occur  and  is  a  function  of  the 
radius  of  curvature,  ihe  refractive  indices  of  the  media 
and  die  height  above  the  plane  which  crosses  the 
centerline  of  the  injector  at  an  angle  of  15  degrees 
(defined  by  the  angle  of  the  incoming  beams). 
Translation  of  the  probe  volume  widiin  die  fuel  annulus 
addressed  diis  refractive  effect  to  determine  accurately 
the  location  of  diita  sampling  within  the  recess  region. 

Measurements  of  mean  axial  velocity  were  taken 
inside  die  injector  recess  region  cts  well  as  downsiremn 
of  the  injector.  Downstream  of  the  injector  face  only 
positive  velocities  were  measured.  Water  was  used  as 
the  LOX  simulant  and  either  iiir  or  gaseous  nitrogen  for 
the  gaseous  fuel  simulant. 

Fig.  8  shows  the  line  iilong  which  die  mean  axial 
velocity  measurements  were  taken.  These  measurements 
are  depicted  in  Fig.  9  as  a  function  of  the  axial  distance 
from  the  LOX  post  base.  The  air  How  rate  was  1.8  g/s 
luid  die  water  flow  rates  were  eidicr  1 .9  or  4.5  g/s.  It  can 
be  seen  that  the  velocity  drops  to  a  negative  value 
approximately  0.5  mm  downstream  of  the  LOX  post 
base  before  becoming  positive.  Larger  negative  values 
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Fig.  8  Injector  clement  orifice  showing  location  of 
velocity  measurement  scans. 


Fig.  9  Droplet  mean  axial  velocity  as  a  function  of 
axial  distance  within  the  injector  recess  region. 

of  the  velocity  were  observed  for  lower  liquid  flow  rates. 
Thus  it  appears  that  gas  recirculation  does  exist  in  the 
recess  region  of  a  shear  coaxial  injector  element  and  that 
the  strength  of  the  recirculation  is  a  function  of  the 
operating  conditions. 


Spray  Visualization. 

Visualization  tests  using  both  techniques  were 
performed  with  the  nontapered  LOX  post  at  chamber 
pressures  between  2.2  and  4.1  MPa  and  at  two  liquid  to 
gas  mixture  ratios  of  0.80  and  1.58.  To  determine  if  any 
observable  differences  in  liquid  atomization  existed 
between  the  two  LOX  post  designs,  specifically  with 
respect  to  the  emergence  of  a  resonating  mode  of 
operation,  several  stroboscope  tests  were  conducted  with 
both  LOX  posts  operating  at  the  same  mixture  ratio  and 
at  a  range  of  chamber  pressures. 

The  first  series  of  photographs  in  Fig.  10 
correspond  to  the  run  conditions  of  tests  A1  to  A3  of 
Table  1  and  illustrate  the  effect  of  chamber  pressure  on 
the  overall  spray  morphology  of  the  coaxial  injector 
using  the  nontapered  LOX  post  recessed  within  the 
annulus.  The  photographs  were  imaged  with  the  laser 
sheet  and  have  been  digitally  reproduced  from  color 
photographs  to  improve  image  quality  for 
photoduplication  and  to  accentuate  noticeable  structures 
within  tlie  spray.  The  pressure  drop  from  2.93  to  2.32 
MPa  Wits  a  result  of  lowering  average  gas  temperatures 
within  the  chamber  from  test  to  test.  It  is  expected  tliat 
in  tlic  case  of  large,  dispersed  sprays  differences  in  the 
gas  temperature  may  influence  (he  behavior  of  the  spray 
downstream  of  the  dense  liquid  core  region  but  should 
not  affect  significandy  the  breakup  behavior  of  the  liquid 


10 


(a)  2.93  MPa  (test  Al) 


(b)  2.60  MPa  (test  A2) 


(c)  2.32  MPa  (test  A3) 


Fig.  10  Effect  of  chamber  pressure  on  LN2/GN2  spray  (mean  mixture  ratio :  0.78). 


jet  or  the  distribution  of  droplets  around  the  jet,  since  the 
conriitinns  in  this  region  are  better  represented  by  the  gas 
injection  temperature. 

In  general.  Fig.  10  depicts  an  overall 
contraction  of  a  dense  liquid  core  region  with  increasing 
chamber  pressure.  Regions  of  droplet  flow  are  most 
evident  downstream  of  this  region  and  decrease  in 
number  density  and  size  with  increasing  pressure.  Based 
on  the  comparatively  small  size  of  the  spray  depicted  in 


Fig.  10,  evaporation  is  considered  an  important 
mechanism  in  the  evolutionary  behavior  of  the  spray.  In 
the  classical  description  of  shear  coaxial  injector 
atomization,  liquid  ligaments  and  droplets  ate  produced 
from  the  effluent  jet  as  a  result  of  unstable  surface  wave 
growths  induced  by  liquid  turbulence  and  a  shearing 
force  between  the  sunounding  gas  and  liquid.  These 
ligaments  and  droplets  in  turn  undergo  secondary 
atomization  through  surface  deformation  and  breakup  by 
aerodynamic  forces  resulting  in  a  spray  of  polydisperse 


(a)  2.79  MPa  (test  A4) 


(b)  2.63  MPa  (test  A5) 


(c)2.34  MPa  (test  A6) 


Fig.  1 1  Effect  of  chamber  pressure  on  LN2/GN2  spray  (mean  mixture  ratio :  1.57). 


Table  3.  Operating  conditions  for  LN2/GN2  spray  visualization  tests. 


Test  No. 

(MPa) 

m,  (kg/s) 

rhg  (kg/s) 

Ug-U,(,mJs) 

PjPi 

A1 

2.93 

WSSM 

0.062 

MBM 

8.99 

164.8 

0.061 

A2 

2.60 

■EH 

0.062 

10.06 

184.3 

0.056 

A3 

2.32 

0.063 

msSM 

11.58 

203.9 

0.052 

A4 

2.79 

0.052 

0.033 

1.58 

4.82 

82.3 

0.060 

A5 

2.63 

0.052 

0.033 

1.58 

5.06 

87.7 

0.057 

A6 

2.34 

0.050 

0.033 

1.54 

5.87 

0.051 

A7 

2.94 

0.076 

0.094 

0.82 

8.68 

0.065 

A8 

3.97 

0.076 

0.094 

0.81 

6.99 

183.6 

0.081 

droplets  around  the  liquid  jet.  With  rising  chamber 
pressure  the  gas  density  enlumccs  die  brejikup  of  tlie 
jet^^,  which  is  observed  in  Fig.  10,  but  reduces  secondary 
atomization  due  to  a  decrease  in  the  slip  velocity 
between  the  gas  and  droplet.  As  a  consequence,  larger 
droplets  pervade  the  spray  of  a  nonevaporating  liquid  at 
higher  chamber  pressures^ ^  although  shorter  liquid 
breakup  lengths  may  exist.  As  seen  in  Fig  10,  tlie  hage 
gas  mass  flux  appears  to  have  rapidly  vaporized  the 
droplets  surrounding  the  core  except  for  those  in  the 
shadow  of  the  LN2  jet,  a  region  characterized  by  lower 
gas  temperatures  and  velocities.  The  rate  of  liquid 
vaporization  is  further  enlnuiced  by  the  increasing 


(a)  2.94  MPa  (test  A7) 

Fig.  1 1  Stroboscope  results  for  tapered  LOX  post 


chamber  pressure  through  a  combined  reduction  in  the 
latent  heat  of  vaporization  and  in  the  liquid  surface 
tension  resulting  in  the  observed  contraction  of  the 
droplet  flow  region. 

In  tests  A4  to  A6  of  Table  1  the  gas  mass  flow 
rate  was  decreased  to  approximately  half  of  the  value  of 
tests  A1  to  A3,  while  the  chamber  pressure  was  varied 
from  2,34  MPa  to  2.79  MPa.  The  photographic  results 
of  these  tests  have  been  digitally  reproduced  in  Fig.  1 1 
and  were  obtained  as  well  with  the  laser  sheet 
visualization  technique.  Overall,  the  reduction  in  gas 
mass  How  rate  has  resulted  in  a  larger,  more  dispersed 

0,0 


2.0 


4.0 


6.0 
(cm) 

(h)  3.97  MPa  (test  A8) 

at  sub-  and  supercritical  pressures  (mixture  ratio  :  0.81), 


Table  4.  Operating  Conditions  and  Droplet  Size  and  Velocity  Data  for  LN2/GN2  PDI  Tests. 


1  Test  No. 

P,  (MPa) 

m,  (kg/s) 

ihg  (kg/s) 

m,lm^ 

f/,-[/,(m/s) 

PjPi 

2.63 

0.073 

0.033 

2.24 

3.87 

86.0 

0.053 

H:9 

2.40 

0.072 

0.033 

2.19 

4.20 

95.5 

0.050 

mm 

231 

0.070 

0.033 

2.15 

4.30 

96.5 

0.049 

B4 

2.41 

0.071 

0.032 

2.20 

4.23 

95.3 

0.050 

B5 

2.45 

0.075 

0.047 

1.60 

5.74 

141.6 

0.050 

B6 

2.37 

0.071 

0.063 

1.14 

7.53 

193.5 

0.053 

Test 

No. 

Pos.  (cm) 

Axial  Velocity  (m/s) 

Droplet  Size  ( (im) 

No.  of 
Droplets 

Percent 

Validation 

Z 

R 

Mean 

Deviation 

Arithmetic 

Sauter 

B1 

12 

0 

15.9 

3.5 

10.1 

96.7 

512 

28 

B2 

12 

0 

14.4 

3.6 

11.2 

97.3 

2514 

25 

B3 

12 

0.5 

9.3 

4.5 

12.9 

91.6 

375 

4 

B4 

10 

0 

14.8 

4.3 

14.1 

105.5 

562 

4 

B5 

10 

0 

20.8 

7.1 

7.5 

110.3 

151 

19 

B6 

6 

0 

37.1 

11.1 

8.4 

118.3 

121 

8 

spray  indicating  a  strong  effect  of  mixture  ratio  on  the 
3^nmi7atinn  and  vaporization  of  liquid  nitrogen.  The 
behavior  of  the  spray  with  respect  to  increasing  chamber 
pressure  follows  the  same  trend  as  tests  A1  to  A3,  in  that 
the  LN2  spray  contracts  in  size  as  the  chamber  pressure 
inaeases.  In  terms  of  droplet  sampling  with  the  phase 
Doppler  interferometric  device,  these  reduced  gas  flow 
conditions  are  advantageous  in  that  the  observable  dilute 
spray  region  is  spread  over  a  larger  volume.  The  droplet 
flow  regions  of  tests  A1  to  A3  appear  confined  close  to 
the  liquid  core  amounting  to  potential  difficulties  in 
droplet  sampling. 

Tests  A7  and  A8  were  performed  with  the 
tapered  LOX  post  using  the  stroboscope  to  visualize  the 
spray.  Test  A7  depicted  in  Fig.  12a  matches  closely  the 
chamber  pressure  and  mixture  ratio  of  test  Al,  but  was 
conducted  with  substantially  higher  liquid  and  gas  mass 
flow  rates.  Regardless  of  this  difference,  the  lengths  of 
the  intact  liquid  core  appear  very  similar  between  the 
two  tests.  It  is  evident  though  that  the  stroboscopic 
image  provides  more  detail  on  the  structure  of  the  Uquid 
core  and  detached  ligaments  than  the  laser  sheet  images 
but  less  information  on  the  droplet  flow  region.  Fig.  12b 
shows  the  structure  of  the  jet  issuing  from  the  injector  at 
a  supercritical  pressure  of  3.97  MPa.  Contradictory  to 
the  observed  contraction  of  the  spray  with  increasing 
pressure,  the  Uquid  core  breakup  length  has  significandy 
increased  for  a  chamber  pressure  greater  than  the  critical 
pressure  (3.4  MPa).  This  behavior  may  be  attributed  to 
the  drop  in  relative  velocity  between  the  gas  and  Uquid  at 
increasing  chamber  pressures,  which  may  surpass  any 
enhancements  to  Uquid  breakup  due  to  increased  gas 
density,  and  may  be  related  as  well  to  the  increase  in 
droplet  evaporation  times  observed  in  the  supercritical 


region^*.  Stroboscope  tests  with  the  nontapered  LOX 
post  under  similar  operating  conditions  revealed  the 
samfi  behavior  for  pressures  above  and  below  the  critical 
point  but  indicated  Uttle  differences  in  spray  morphology 
between  the  two  LOX  post  designs. 

Droplet  Size  and  Velocity  Measurements. 

For  droplet  sizing  and  velocity  measurement 
only  the  nontapered  LOX  post  injector  design  was 
examined  at  liquid  to  gas  mixture  ratios  of 
approximately  2.2,  1.6  and  1.1  for  a  nominal  Uquid 
nitrogen  flow  rate  of  0.072  kg/s.  Chamber  pressures 
ranged  from  2.20  to  3.29  MPa.  Based  on  past  difficulty 
in  obtaining  a  large  statisUcal  sampUng  of  droplets  with 
the  injector  operating  with  a  mixture  ratio  close  to  unity, 
the  present  effort  began  with  the  gas  mass  flow  rate  at 
almost  half  of  the  Uquid  mass  flow  rate.  According  to 
the  photographic  results  of  the  flow  visualization  tests 
(Fig.  11),  these  conations  result  in  a  dilute  droplet  flow 
region  dispersed  over  a  larger  volume,  which  effectively 
reduces  both  the  droplet  number  density  and  signal 
attenuation  due  to  multiple  particle  scattering.  The 
average  gas  and  Uquid  mass  flow  rates  were  0.032  and 
0.072  kg/s,  respectively,  yielding  a  Uquid  to  gas  mixture 
ratio  of  approximately  2.20.  Fig.  13  iUustrates  the  effect 
of  chamber  pressure  on  the  number  of  droplets  collected 
over  a  2.5  second  sampling  interval  at  12  cm 
downstream  of  the  injector  on  the  injector  centerUne. 
With  an  increase  in  chamber  pressure  the  number  of 
droplets  sampled  over  this  interval  decreased 
dramatically.  This  trend  would  suggest  that  the  droplet 
number  density  has  declined  due  to  the  contraction  of  the 
'  spray  at  higher  pressures  as  observed  in  Fig.  11  and/or 
that  the  optical  alignment  has  been  affected  by  the 
increased  gas  density  (i.e.  changing  refractive  index). 
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Fig.  13, Effect  of  chamber  pressure  on  number  of 

droplets  smnpled  over  2.5  sec  interval. 

The  effect  of  chamber  pressure  on  droplet  size  and  mean 
iixiid  velocity  was  found  to  be  negligible  as  is  seen  in  tlie 
results  of  tests  B 1  ;uid  B2  of  Table  2. 

Additional  tests  were  performed  to  examine  the 
effect  of  position  mid  gas  mass  flow  rate  on  tlic  droplet 
siimpling  rate  and  characteristic  parameters.  Tests  B2, 
B3,  and  B4  provide  a  comparison  of  the  spray 
measurements  taken  at  three  different  locations  in  the 
spray  for  approximately  equal  operating  conditions. 
Best  data  was  obtained  at  12  cm  downstream  from  the 
injector  on  the  centerline  as  seen  by  the  large  number  of 
droplets  and  good  ratio  of  valid  droplet  samples  to 
attempted  samples  (percent  validation)  of  test  B2.  For 
the  same  axial  distance  away  trom  the  injector  but  5  mm 
off  axis,  the  number  of  droplets  acquired  over  the 
siunpling  interval  decreased.  The  droplet  data  show 
tliat  tlie  droplet  mean  axiid  velocity  decreased  in  the 
radial  direction  while  the  droplet  size  represented  by  the 
arithmetic  mean  diameter  changed  negligibly.  The 
results  of  test  B4  for  an  axial  position  of  10  cm  exhibit 
the  expected  increase  in  droplet  mean  velocity  and 
droplet  size  when  probing  closer  to  the  injector  face. 


The  effect  of  gas  mass  flow  rate  is  seen  from 
the  results  of  tests  B4,  B5  and  B6.  Test  B5  was 
conducted  at  a  higher  gas  flow  rate  of  0.047  kg/s  at  the 
same  location  as  test  B4.  The  flow  rate  was  further 
increased  to  0.063  kg/s  in  test  B6.  But,  due  to  the 
marked  contraction  of  the  spray  at  this  gas  flow  rate, 
measurements  were  made  at  6  cm  downstream  from  the 
injector  for  this  test  The  results  of  test  B5  and  B6  reveal 
the  expected  increase  in  droplet  mean  velocity  with 
increasing  gas  mass  flow  rate  due  to  the  higher  gas 
injection  velocities.  Although  the  low  number  of 
droplets  collected  in  the  last  two  tests  does  not  permit  an 
accurate  evaluation  of  the  mean  droplet  size  for  all  three 
tests,  the  droplet  size  distributions  have  been 
superimposed  in  Fig  14  to  illustrate  the  general  behavior 
of  tlie  spray  with  varying  liquid  to  gas  mixture  ratios  At 
the  lower  mixture  ratio  of  test  B6,  the  majority  of 
droplets  fall  in  the  smallest  size  range  of  Fig.  14.  As  the 
gas  mass  flow  rate  decreases  to  approximately  3/4  and 
1/2  of  that  in  test  B6,  the  droplet  distribution  shifts  to 
hirger  droplet  sizes.  In  the  case  of  tests  B4  and  B5, 
wliich  were  both  conducted  at  the  same  position  in  the 
spray,  tiiis  trend  confirms  die  reduction  in  droplet  size 
and  number  due  to  enhanced  vaporization  and 
atomization  with  increasing  gas  mass  flow  rate.  The 
large  difference  between  the  inithmetic  mean  diameter 
and  Sauter  mean  diameter  in  the  tests  arose  from  the 
bimodal  distribution  in  the  sampled  droplets  as  seen  in 
Fig.  14.  The  small  number  of  large  droplets  at  the  upper 
end  of  the  size  range,  observed  during  all  of  the  tests, 
may  be  due  to  remnants  of  the  liquid  core  and/or 
coalesced  small  droplets.  As  seen  in  Table  4,  their  effect 
on  the  Sauter  mean  diameter  is  most  pronounced  with 
decreasing  mixture  ratio,  although  dieir  number  is  lowest 
at  these  mixture  ratios. 


The  injector  response  factor  was  calculated  as  a 
function  of  frequency  and  gaseous  fuel  temperature  for 
the  experimental  conditions  experienced  during  the  cold 
flow  LN2/GN2  experiments.  During  these  experiments, 
high  frequency  pressure  oscillations  were  measured  in 
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Fig.  14,  Droplet  size  distribution  tor  three  oxidizer  to  fuel  mixture  ratios  (Pc  =  2,4  MPa). 
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the  chamber  and  the  injector  fuel  plenum.  Figure  15 
shows  the  injector  response  (real  part)  as  a  lunction  of 
fuel  temperature  for  several  frequencies.  It  can  be  seen 
that  at  high  frequencies  the  response  is  near  zero,  at 
intermediate  frequencies  the  response  is  negative  and 
almost  constant  and  at  low  frequencies  the  response 
becomes  strongly  positive  as  the  fuel  temperature 
decreases.  Fig.  16  shows  the  injector  response  as  a 
function  of  frequency  for  the  gas  temperature  used 
during  the  experiments.  The  response  is  positive  tor  low 
frequencies  but  then  switches  to  negative  before 
decaying  to  zero  as  the  frequency  increases. 

The  recorded  high  frequency  pressure  oscillations 
were  spectrally  analyzed  to  determine  the  frequencies 
and  magnitudes  of  any  acoustic  oscillations  within  either 
the  chamber  or  the  injector  fuel  plenum.  Oscillations  at 
the  same  frequency  in  both  the  chamber  and  the  plenum 
were  found  at  two  frequencies.  Table  5  compares  the 
calculated  injector  response  at  those  two  frequencies  to 
the  measured  ratios  of  the  pressure  oscillation  amplitude 
in  the  chamber  to  the  amplitude  in  the  plenum.  It  can  be 
seen  that  at  the  frequency  where  the  injector  was  more 
strongly  damping,  the  measured  ratio  of  amplitudes  was 
less.  Future  work  for  both  the  analysis  and  the 
experiments  will  include  the  effects  of  combustion. 

Table  5  Comparison  of  calculated  injector  response  to 
the  ratio  of  measured  pressure  oscillation  in  the 
chamber  to  the  amplitude  in  the  plenum. 

Frequency  Injector  ^  Measmed 

(Hz) _  Response  P'cham^'olenum 

- 2000  ^o!25  ol6 

5800  -0.04  0.62 


CONCLUSIONS 

An  experimental  investigation  of  the  spray 
produced  by  a  single  shear  coaxial  injector  has  been 
performed  under  noncombusting  conditions.  Two  fluid 
combinations  were  selected  for  spray  characterization 
using  stfoboscopic  photography  and  planar  laser  imaging 
to  visualize  the  global  spray  characteristics  and  phase 
Doppler  interferometry  to  measure  droplet  size  and 
velocity.  Initial  tests  were  performed  with  water  and  air 
at  atmospheric  pressure.  Using  the  tapered  LOX  post 
recessed  into  the  fuel  annulus,  the  injector  transitioned 
into  a  condition  of  resonance  characterized  by  a 
whistling  noise  and  a  significant  modification  to  the 
overall  structure  of  the  spray.  Injector  operation  with  the 
nontapered  post  produced  a  near  resonance  condition 
that  exhibited  a  lower  amplitude  whistling  noise  but  did 
not  effect  a  change  in  the  spray  morphology.  No 
evidence  of  injector  resonance  was  encountered  upon 
elimination  of  the  LOX  post  recess.  Stroboscopic 
images  emphasized  the  magnitude  of  this  phenomenon 
on  the  atomization  characteristics  of  the  liquid  jet  with 
two  post/annulus  configurations.  Measurements  of 
droplet  size  and  velocity  at  a  single  cross-section  of  the 
spray  illustrated  a  significant  quantitative  difference 
between  nonresonant  atomization  and  resonant 
atomization. 

spray  visualization  tests  with  LN2/GN2  at 
elevated  pressures  illustrated  the  strong  influence  of 
liquid  to  gas  mixture  ratio  and  chamber  pressure  on  the 
overall  morphology  of  the  spray.  Higher  mature  ratios 
yielded  large,  poorly  atomized  and  vaporized  s^ys 
while  increasing  chamber  pressures  below  the  critical 
point  resulted  in  a  contraction  of  the  spray.  For  a 
mixture  ratio  near  unity  and  pressures  above  the  critical 
point  the  spray  ceased  to  contract  and  extended  deeper 
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into  the  chamber.  Measurements  ot  droplet  size  and 
mean  velocity  as  a  function  of  mixture  ratio  and  chamber 
pressure  have  been  performed.  The  results  confirm  the 
improved  atomization  of  the  spray  with  decreasing 
mixture  ratio,  but  reveal  potential  concerns  of  droplet 
sampling  in  dense  cryogenic  sprays  at  elevated  pressures 
and  reduced  temperatures. 

Results  from  a  linearized  lumped-element 
model  of  the  response  for  a  cryogenic  shear  coaxial 
injector  show  that  the  injector  response  may  shift  from 
stabilizing  to  destabilizing  as  the  gaseous  fuel 
temperature  decreases  for  certain  frequencies.  The 
analytical  results  compared  well  to  chamber  high 
frequency  pressure  measurements. 

The  first  velocity  measurements  ever  taken  in 
die  recess  region  of  a  recirculation  region  of  m  injector 
element  showed  reverse  flow  indicative  of  a 
recirculation  region.  Measurements  will  be  expanded  to 
quantify  the  size  and  strength  of  the  recirculation  region 
as  a  function  of  operating  conditions  iind  injector  design. 


The  pressure  at  the  exit  of  the  annulus  was  assumed  to 
equal  the  mean  chamber  pressure. 

These  equations  may  be  combined  to  give  the 
transfer  function: 
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Inverting  tliis  transfer  function  gives  the  following 
amplitude  iuid  phase  shift: 
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APPENDIX 


Rgt>i?.Qn^e,Eat;M:a 

A  lumped-element  model  was  developed  simihir 
to  Reference  22  to  determine  the  dynamic  response  of 
die  fuel  flow  through  the  injector.  The  dynamics  of  the 
flow  dirough  the  dome  is  represented  by  a  mass  balance. 
In  perturbation  form,  this  relation  is: 

^sFj=-\V  (4) 


The  orifice  flow  relation  represents  the  flow  from  the 
dome  into  die  entrance  of  the  annulus: 
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The  flow  through  die  annulus  is  represented  by  a 
momentum  biiltUicc  assuming  plug  flow: 
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Assuming  adiabatic  flow,  die  response  factor  for 
the  nozzle  can  be  expressed  as: 
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(12) 
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ABSTRACT 


Vaporization  and  combustion  of  liquid  droplets  in  both  subcritical  and  su¬ 
percritical  environments  have  been  studied  systematically.  A  variety  of  liquid 
propellants  and  propellant  simulants,  including  hydrocarbon  and  cryogenic  fluids, 
in  both  steady  and  oscillatory  environments  are  treated.  The  formulation  is  based 
on  the  full  conservation  equations  for  both  gas  and  liquid  phases,  and  accommo¬ 
dates  variable  properties  and  finite-rate  chemical  kinetics.  Full  account  is  taken  of 
thermodynamic  non-idealities  and  transport  anomalies  at  high  pressures,  as  well  as 
liquid/vapor  phase  equilibrium  for  multi- component  mixtures.  Because  the  model 
allows  solutions  from  first  principles,  a  systematic  examination  of  the  droplet  be¬ 
havior  over  wide  ranges  of  pressure,  temperature,  and  ambient  flow  velocity  is 
made  possible.  Results  can  not  only  enhance  the  basic  understanding  of  the  prob¬ 
lem,  but  also  serve  as  a  basis  for  establishing  droplet  vaporization  and  combustion 
correlations  for  the  study  of  liquid  rocket  engine  combustion,  performance,  and 
stability. 

A  series  of  calculations  have  been  performed  to  understand  the  effects  of 
ambient  flow  conditions  on  droplet  gasification  behavior.  Results  indicate  that  the 
velocity  and  thermodynamic  state  of  the  ambient  flow  have  strong  influence  on  the 
mass,  momentum,  and  energy  transport  in  the  droplet  gasification  and  burning 
processes.  The  droplet  gasification  rate  increases  progressively  with  pressure  and 
ambient  Reynolds  number.  In  addition,  the  amplitude  of  pressure-coupled  droplet 
vaporization  response  increases  with  increasing  pressure  owing  to  the  susceptibility 
of  enthalpy  of  vaporization  to  ambient  flow  oscillations  at  high  pressures.  However, 
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the  effect  of  mean  pressure  on  the  phase  angle  of  the  droplet  vaporization  response 
appears  quite  limited. 

Detailed  flow  structures  and  thermodynamic  property  variations  are  exam¬ 
ined  to  reveal  underlying  mechanisms  for  droplet  gasification  and  burning,  as  well 
as  deformation  and  breakup  dynamics  at  supercritical  pressures.  Correlations  of 
droplet  lifetime  and  aerodynamic  drag  coefficient  are  developed  as  functions  of 
fluid  thermodynamic  state,  Reynolds  number,  and  vaporization  transfer  number. 
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CHAPTER  1 


INTRODUCTION 


Many  practical  liquid-propellant  rocket  propulsion  systems  involve  droplet 
vaporization  and  spray  combustion  in  high-pressure  environments [1-3].  Liquid 
propellants  are  usually  delivered  to  liquid-rocket  combustion  chambers  as  a  spray 
of  droplets,  which  then  undergo  a  sequence  of  vaporization,  ignition,  and  com¬ 
bustion  processes  at  pressure  levels  well  above  the  thermodynamic  critical  points 
of  the  liquids.  Under  these  conditions,  droplets  initially  injected  at  subcritical 
temperatures  may  heat  up  and  experience  a  thermodynamic  state  transition  into 
the  supercritical  regime  during  their  lifetimes,  as  shown  in  Fig.  1.1.  Consequently, 
the  sharp  distinction  between  gas  and  liquid  disappears,  and  the  entire  system  ex¬ 
hibits  many  distinct  characteristics  which  conventional  droplet  theories  developed 
for  low-pressure  cases  can  not  deal  with. 

To  clarify  this  problem,  it  is  useful  to  examine  the  behavior  of  an  isolated 
liquid  droplet  when  suddenly  confronted  with  a  stream  having  a  thermodynamic 
state  in  the  supercritical  regime  of  the  liquid  (see  Figure  1.2  ).  As  the  droplet  is 
heated  by  the  ambient  gas,  its  temperature  increases  and  finally  exceeds  the  critical 
point.  Several  important  aspects  must  be  noted  during  this  process.  First,  when 
the  droplet  surface  approaches  its  thermodynamic  critical  state,  the  difference  in 
densities  of  gas  and  liquid  phases  become  smaller.  The  characteristic  times  of  the 
transport  processes  near  the  droplet  surface  in  both  phases  have  the  same  order  of 
magnitude.  Therefore,  the  transient  effects  in  the  gas  phase  are  as  important  as 
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Figure  1.1.  Operating  Conditions  of  Advanced  Air-breathing  and  Rocket  Engines. 


those  in  the  liquid  phase;  the  quasi-steady  condition  may  never  be  reached  during 
the  lifetime  of  the  droplet [4].  Second,  the  latent  heat  of  vaporization  decreases 
to  zero  at  the  critical  point.  Conventional  low-pressure  models  may  erroneously 
predict  the  vaporization  rate  if  the  variation  of  latent  heat  with  pressure  is  not 
properly  taken  into  account.  In  addition,  if  the  droplet  is  moving,  the  behavior 
of  liquid  deformation  and  breakup  may  be  altered  considerably  due  to  the  dimin¬ 
ished  value  of  surface  tension.  Third,  at  high  pressure,  effects  of  thermodynamic 
non-idealities  and  property  variations  play  decisive  roles  in  determining  transport 
properties  and  interface  thermodynamic  relationships.  The  solubility  of  the  am¬ 
bient  gas  in  the  liquid  phase  increases  with  pressure;  the  classical  Raoult’s  law 
for  ideal  mixtures  is  not  applicable  in  the  phase-equilibrium  analysis.  One  must 
develop  a  more  comprehensive  model  for  vapor-liquid  interface  conditions  in  terms 
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Figure  1.2.  Liquid  Spray  from  an  Co-axial  Injector. 


of  fugacity.  Finally,  when  the  droplet  exceeds  its  critical  state,  it  essentially  be¬ 
comes  a  puff  of  dense  fluid.  The  entire  field  becomes  a  continuous  medium,  and 
no  distinct  liquid/vapor  interface  can  be  identified. 

1.1  Literature  Survey 

Attempts  to  study  droplet  vaporization  and  combustion  in  the  supercritical 
regime  have  been  made  for  three  decades.  The  earliest  development  of  a  predictive 
model  was  initiated  by  Spalding[5].  He  considered  the  problem  of  supercritical 
combustion  of  an  isolated  fuel  droplet  in  an  stagnant  environment  by  approximat¬ 
ing  the  droplet  as  an  instantaneous  point  source  of  dense  gas.  The  formulation 
was  based  on  a  simplified  one-dimensional  model  with  constant  physical  properties; 
real-gas  effects  and  density  variations  were  ignored.  The  analysis  was  then  refined 
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by  Rosner[6]  to  include  finite  distributed  sources.  Results  of  both  studies,  substan¬ 
tially  in  agreement,  demonstrated  the  transient  characteristics  of  droplet  burning 
under  supercritical  conditions.  The  same  problem  was  examined  by  a  number  of 
researchers  (notably  Brzustowski[7],  Chervinsky[8],  and  Polymeropoulos  and  Pe- 
skin[9])  in  order  to  investigate  the  influences  of  convection,  density  variation,  and 
finite-rate  chemical  kinetics.  Recently,  Lee  et  a/.[10]  investigated  the  stagnant- 
point  evaporation  of  a  liquid  fuel  at  near  and  super-critical  conditions.  Because  of 
the  complexity  of  this  problem,  many  important  high-pressure  phenomena  (includ¬ 
ing  solubility  of  ambient  gases  in  the  liquid  phase,  thermodynamic  non-idealities, 
and  property  variations)  were  ignored. 

In  1963,  Wieber[ll]  proposed,  in  accordance  with  a  low-pressure  quasi- 
steady  model,  that  a  droplet  may  undergo  a  certain  amount  of  evaporation  prior 
to  reaching  the  critical  condition.  This  observation  was  later  confirmed  experimen¬ 
tally  by  Faeth  et  a/. [12]  in  a  study  of  supercritical  bipropellant  droplet  combustion. 
Systematic  treatment  of  droplet  vaporization  at  near  critical  conditions  was  initi¬ 
ated  by  Manrique  and  Borman[13]  and  Savery  and  Borman[14],  based  on  a  quasi¬ 
steady  model.  They  concluded  that  the  effects  of  thermodynamic  non-idealities, 
property  variations,  and  high-pressure  corrections  for  phase  equilibrium  modify  the 
vaporization  mechanism  significantly.  In  particular,  an  appreciable  amount  of  am¬ 
bient  gas  may  dissolve  in  the  liquid  phase,  causing  an  intrinsic  change  in  the  entire 
process.  In  light  of  these  findings,  Lazar  and  Faeth[15]  and  Canada  and  Faeth[16] 
conducted  a  series  of  experimental  and  theoretical  studies  on  droplet  combustion  in 
both  stagnant  and  forced-convective  environments,  with  special  attention  focused 
on  the  high-pressure  phenomena  of  phase  equilibrium.  The  effects  of  forced  con¬ 
vection  in  the  gas  phase  were  treated  by  conventional  multiplicative  corrections. 
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The  assumptions  of  quasi-steadiness  and  uniform  property  distributions  adopted  in 
the  analyses  of  Borman  and  Faeth  were  later  relaxed  to  include  several  important 
aspects  of  droplet  vaporization  and  combustion.  In  works  by  Matlosz  et  a/.  [17], 
Rosner  and  Chang[18],  Kadata  and  Hiroyasu[19],  and  Curtis  and  Farrel[20],  the 
effects  of  transient  processes,  natural  convection,  and  the  conditions  under  which 
a  droplet  may  be  driven  to  its  critical  point  were  examined. 

Many  practical  liquid  propellants  consist  of  multicomponent  mixtures.  Even 
for  pure  cryogenic  propellants  such  as  liquid  oxygen,  the  solubility  of  ambient 
gases  in  the  liquid  phase  causes  the  droplet  to  become  a  multicomponent  system. 
As  a  result  of  different  volatilities  and  diffusivities  of  the  liquid  constituents,  the 
vaporization  and  combustion  processes  present  several  features  distinct  from  those 
of  a  simple  component  system.  The  major  problem  in  treating  multicomponent 
systems  is  to  specify  the  state  of  mixedness  of  the  liquid  phase.  Mass  diffusion 
in  the  liquid  phase  is  usually  slower  than  heat  diffusion;  therefore,  convection  and 
transient  effects  are  potentially  more  important  for  mass  transfer  than  for  heat 
transfer.  At  high  pressures,  the  problem  becomes  even  more  difficult,  due  to  the 
dissolution  of  ambient  gases  in  the  liquid  phase  and  rapid  variations  of  physical 
properties  near  the  critical  point.  Jin  and  Borman[21]  carried  out  a  simplified 
analysis  of  multicomponent  fuel  droplet  vaporization  at  elevated  pressures  and 
temperatures,  noting  the  significance  of  liquid-phase  transport.  Several  calculated 
results  were  given  for  pentane/octane  mixtures.  The  potential  for  this  particular 
fuel  to  undergo  microexplosions  was  also  discussed  briefly. 

Hsieh  et  al.[4]  recently  developed  a  comprehensive  analysis  to  investigate 
hydrocarbon  fuel-droplet  vaporization  at  near-critical  conditions.  The  model  is 
based  on  full  time-dependent  conservation  equations,  with  thorough  treatment  of 
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property  variations  and  high-pressure  vapor-liquid  phase  equilibrium.  Because  of 
its  completeness,  the  model  enables  systematic  examination  of  droplet  vaporiza¬ 
tion  in  a  high-pressure  environment.  Results  of  Hsieh’s  work  indicate  that  ambient 
pressure  strongly  influences  the  vaporization  process  and  that  the  droplet  evap¬ 
oration  rate  increases  progressively  with  pressure.  In  the  extreme  case,  droplets 
may  reach  their  critical  states  before  the  end  of  their  lifetimes,  thereby  causing  a 
significant  change  in  the  vaporization  process. 

Although  the  preceding  studies  have  provided  significant  information  to 
understanding  the  physics  and  chemistry  of  droplet  vaporization  and  combustion 
at  high  pressures,  a  number  of  fundamental  problems  remain  unresolved.  First, 
much  of  the  previous  work  employed  certain  assumptions  and  empirical  correla¬ 
tions  which  were  extrapolated  from  low-pressure  causes  and  involved  a  considerable 
number  of  uncertainties.  Second,  no  unified  treatment  of  the  entire  droplet  his¬ 
tory,  in  particular  the  transition  from  subcritical  to  supercritical  state,  was  made. 
Almost  all  of  the  models  for  supercritical  droplet  vaporization  and  combustion 
assumed  that  the  droplet  reaches  its  critical  state  instantaneously  upon  intro¬ 
duction  into  a  supercritical  environment.  Neither  initial  heat-up  transients  nor 
non-uniformities  of  property  distributions  in  the  droplet  were  taken  into  account 
in  these  analyses.  Third,  very  little  is  known  about  multicomponent  systems, 
ignition  transients,  and  effects  of  forced  convection.  It  is  difficult,  therefore,  to 
apply  the  results  for  a  single  droplet  to  a  spray  environment  without  introducing 
errors.  Fourth,  almost  no  attempt  was  made  to  investigate  the  dynamic  responses 
of  droplet  transport  to  ambient  flow  oscillations.  In  the  study  of  combustion  insta¬ 
bilities,  these  responses  are  particularly  important  since  the  interactions  between 
droplet  vaporization/combustion  and  acoustic  wave  oscillations  represent  one  of 
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the  major  (if  not  the  most  significant)  mechanisms  for  driving  unsteady  motions 
in  a  rocket  motor. 

In  order  to  correct  the  deficiencies  of  existing  models  for  high-pressure 
droplet  vaporization  and  combustion,  the  following  research  topics  are  required: 

1.  comprehensive  analysis  of  multicomponent  droplet  vaporization  and 
combustion  in  both  quiescent  and  forced-convective  environments,  in¬ 
cluding  calculations  of  detailed  physical  and  chemical  processes; 

2.  examination  of  the  effects  of  droplet  interactions  on  interface  mass,  mo¬ 
mentum,  and  energy  transfer  at  high  pressure; 

3.  establishment  of  correlations  for  droplet  vaporization  and  combustion 
in  terms  of  ambient  flow  conditions  and  droplet  states; 

4.  investigation  of  droplet  ignition  transients;  and 

5.  study  of  the  dynamic  responses  of  droplet  transport  rates  to  ambient 
flow  oscillations. 

1.2  Research  Objectives 

The  physical  and  chemical  mechanisms  associated  with  droplet  vaporiza¬ 
tion  and  combustion  at  high  pressure  are  extremely  complicated.  Many  intricate 
flow,  transport,  and  combustion  processes  occur  in  various  parts  of  the  liquid-  and 
gas-phase  regions.  In  view  of  the  inadequacy  of  existing  theories  and  experimen¬ 
tal  techniques,  a  theoretical  investigation  of  the  problem  using  direct  numerical 
simulation  methods  is  conducted  here. 
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The  overall  objective  of  this  work  is  to  answer  some  of  the  fundamental 
questions  regarding  liquid-propellant  droplet  vaporization  and  combustion  in  a 
high-pressure  environment.  Specific  program  objectives  are: 

1.  to  acquire  basic  understanding  of  physical  and  chemical  mechanisms 
involved  in  the  vaporization  and  combustion  of  isolated  multicompo¬ 
nent  liquid-propellant  droplets  in  both  stagnant  and  forced- convective 
environments; 

2.  to  establish  droplet  vaporization  and  combustion  correlations  for  the 
study  of  liquid-propellant  spray  combustion  and  two-phase  flow  fields 
in  rocket  engines; 

3.  to  investigate  dynamic  responses  of  droplet  vaporization  to  ambient  flow 
oscillations;  and 

4.  to  incorporate  the  droplet  vaporization  response  functions  developed  in 
this  project  into  existing  stability  analyses  for  examining  the  unsteady 
behavior  of  a  rocket  engine. 

The  present  theoretical  research  represents  a  series  of  investigations  into 
the  behavior  of  liquid- propellant  droplets  at  elevated  pressures  and  temperatures. 
In  addition  to  contributions  to  the  knowledge  base  of  the  fundamental  mecha¬ 
nisms  of  droplet  vaporization  and  combustion,  results  can  be  applied  eff^ectively  to 
existing  spray  models  for  engine  performance  and  stability  predictions.  In  the  fol¬ 
lowing  chapters,  droplet  behavior  for  both  hydrocarbon  and  cryogenic  propellants 
are  investigated.  The  theoretical  models,  including  conservation  laws  of  physics 
and  a  unified  evaluation  scheme  for  thermophysical  properties,  are  first  described 
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in  Chapter  2.  Since  all  the  governing  equations  are  strongly  coupled  together,  a 
robust  numerical  scheme  capable  of  treating  transient,  low  Mach-number  compress¬ 
ible  flows  is  adopted.  Chapter  3  explains  this  numerical  method  and  associated 
discretization  procedure.  Chapter  4  focuses  on  droplet  vaporization  and  com¬ 
bustion  in  quiescent  environments  where  the  pressure  and  temperature  are  either 
above  or  below  the  critical  points  of  the  propellants.  Effects  of  ambient  pressure 
on  droplet  behavior  are  studied  systematically.  Many  intricate  flow,  transport, 
and  combustion  processes  occurring  in  various  parts  of  the  liquid  and  gas-phase 
regions  are  discussed  in  detail.  Dynamic  responses  of  liquid-fuel  droplet  vapor¬ 
ization  to  ambient  pressure  oscillations  are  examined  in  Chapter  5.  The  analysis 
extends  the  model  of  droplet  combustion  described  in  Chapter  4,  and  includes  a 
periodic  pressure  oscillation  in  the  ambient  flow.  The  purpose  is  to  provide  a  de¬ 
tailed  understanding  of  pressure-coupled  droplet  vaporization  response,  including 
the  amplification  gain  and  time  delay  due  to  the  transient  processes  involved  in 
the  vaporization  process.  Chapter  6  investigates  the  effects  of  flow  convection  on 
supercritical  droplet  gasification  mechanisms.  The  work  treats  multi-dimensional 
flowfields  in  both  the  droplet  interior  and  surrounding  gases.  Because  the  model  is 
based  on  the  complete  conservation  laws  of  mass,  momentum,  energy,  and  species, 
several  important  transient  phenomena,  including  droplet  dynamics  and  interphase 
transport,  can  be  thoroughly  studied.  Finally,  conclusions  of  the  present  work  are 
summarized  in  Chapter  7. 


CHAPTER  2 


THEORETICAL  FORMULATION 

2.1  Droplet  Behavior  in  Quiescent  Environments 

Vaporization  and  combustion  of  an  isolated  liquid  droplet  in  a  quiescent 
environment  are  first  considered.  This  configuration  allows  researchers  to  focus 
on  the  effects  of  pressure  on  the  thermodynamic  and  transport  processes  involved. 
The  initial  temperature  of  the  droplet  is  subcritical,  but  the  ambient  tempera¬ 
ture  and  pressure  are  in  the  supercritical  regime  of  the  liquid  species,  as  shown  in 
Fig.  2.1(a).  As  a  result  of  heat  transfer  from  the  surrounding  gases,  the  surface 
temperature  increases,  and  may  reach  its  critical  mixing  point  during  the  droplet 
lifetime.  Prior  to  occurrence  of  the  critical  mixing  condition,  the  droplet  sur¬ 
face  provides  a  well-defined  interfacial  boundary  which  separates  the  liquid  phase 
from  the  ambient  gases.  To  facilitate  formulation,  these  two  regions  are  treated 
separately,  then  linked  together  at  the  droplet  surface  by  requiring  liquid-vapor 
phase  equilibrium  and  continuities  of  mass  and  energy.  The  procedure  eventually 
determines  the  droplet  vaporization  rate  and  interface  conditions,  e.g.,  surface  tem¬ 
perature,  species  concentrations,  and  enthalpy  of  vaporization.  When  the  surface 
becomes  critical,  the  sharp  distinction  between  gas  and  liquid  vanishes.  The  latent 
heat  of  vaporization  reduces  to  zero,  with  no  abrupt  phase  change  involved  in  the 
vaporization  process.  The  density  and  temperature  of  the  entire  field  (including 
both  liquid  droplet  and  ambient  gases)  as  well  as  their  gradients  vary  continu¬ 
ously  across  the  droplet  surface.  The  droplet  interior,  however,  remains  at  the 


liquid  state  with  a  subcritical  temperature  distribution.  Under  this  situation,  the 
droplet  regression  is  best  characterized  by  the  motion  of  the  surface  which  attains 
the  critical  mixing  temperature  of  the  system. 


2.1.1  Governing  Equations 


The  analysis  is  based  on  the  time-dependent  conservation  laws  of  physics  for 
a  multicomponent  system.  For  problems  involving  moving  boundaries  and  com¬ 
plicated  two-phase  flows,  the  finite-volume  formulation  offers  definite  advantages 
over  the  differential  formulation  due  to  the  physical  clarity  of  the  method.  If  body 
forces,  viscous  dissipation  and  radiation  are  ignored,  the  dynamic  behavior  of  the 
flowfield  can  be  described  by  the  following  formulation. 

Mass: 


Momentum: 


d 


dt  IaV  pUr{Ur  —  Ug)dA  -t-  TrrdA 


(2.2) 


Energy: 


^  petdV  +  pet{ur  —  Ug)dA  -|-  pUrdA  =  q^dA  (2.3) 


Species  Concentration: 

I  pYidV  +  /  pYi{ur  -  Ug)dA  =  [  qm,idA+  f  widV 

Jav  Ja  Ja  Jav 


dt 


(2.4) 


where  the  surface  and  volume  integrations  are  carried  out  over  a  finite  spherical 
control  volume  concentric  with  the  droplet  (see  Fig.  2.2).  The  physical  variables 


Hydrocarbon 

Droplet 


(b) 


Hydrogen 


=  10-140  atm 
=  1500K 
=  300K 
=  100, 300  fim 


=  100-400  atm 
=  1000K 
=  100K 
=  100  ^im 
=  0.1-15  m/s 


Figure  2.1.  Schematic  Diagrams  of  Droplet  Vaporization  and  Combustion:  (a) 
Quiescent  Environments,  (b)  in  Forced  Convective  Environments. 
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Size  of  Control 


Figure  2.2.  Schematic  of  a  Spherical  Droplet. 


p,  Ur,  p,  Yi,  and  Wi  are  density,  radial  velocity,  pressure,  species  concentration  and 
rate  of  production  of  species  i,  respectively,  and  Ug  the  grid  moving  velocity.  The 
specific  total  internal  energy  is  written  as 

et  =  h-^^\ul  (2.5) 

P  2 

with  h  the  specific  enthalpy.  In  a  spherically  symmetric  coordinate,  the  normal 
stress  Trr  takes  the  form 


d  (  2  dur\  2p,Ur 

ssr  '‘arj ” 


(2.6) 


where  p  is  the  viscosity  and  r  the  radial  coordinate  relative  to  the  droplet  cen¬ 
ter.  Thermal  and  species  diffusions  are  evaluated  using  Fourier’s  and  Fick’s  laws, 
respectively. 
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9771,2  — 


pT^im 


dr 


(2.8) 


where  the  index  N  is  the  number  of  species  considered,  and  T,  A,  hi,  and  Vim  the 
temperature,  thermal  conductivity,  specific  enthalpy  and  effective  mass  diffusion 
coefficient  for  species  i,  respectively. 

A  moving  grid  system  is  employed  to  simplify  the  problem  associated  with 
droplet  surface  motion.  In  the  gas  phase,  all  the  finite- volume  cells  move  uniformly 
toward  the  center  at  a  speed  of  Ug,  taken  to  be  the  droplet  surface  regression 
rate,  Ug  =  dR/dt.  In  the  liquid  phase,  owing  to  the  finite  radius  of  the  droplet, 
a  time-varying  grid  system  is  used,  with  the  grid  moving  velocity  equal  to  the 
instantaneous  local  flow  velocity. 

For  convenience,  the  conservation  equations  are  collected  into  a  vector  form 

f  QdV+  f  (E-E„)dA-Ug  (  QdA=  (  SdV  (2.9) 

dtJAV  Ja  ^Ja  Jav 

where 
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All  the  variables  of  Q  in  the  vector  equation  are  treated  simultaneously  in  a  fully 


coupled  manner. 


2.1.2  Boundary  Conditions 


Liquid  Phase 

A  symmetric  boundary  condition  is  implemented  at  the  droplet  center(r  = 
0).  That  is,  the  radial  velocity  and  gradients  of  pressure,  temperature  and  species 
concentration  must  be  zero  at  the  droplet  center. 


Ut 

=  0 

r=0 

dp 

dT 

_  dYi 

dr 

r=0  ~ 

dv 

r=0 

r=0 

Gas  Phase 


(2.10) 

(2.11) 


Since  the  convective  velocity  is  negligible  small  in  the  far  field,  the  pressure  is 
set  to  be  the  ambient  pressure.  In  addition,  the  gradients  of  temperature,  velocity. 
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and  species  concentration  are  zero.  In  mathematical  form,  they  can  be  expressed 
as  follows. 


dur 


dr 


(2.12) 

(2.13) 


2.1.3  Matching  Conditions  at  Liquid/ Vapor  Interface 


The  physical  processes  in  each  phase  must  be  linked  at  the  droplet  surface 
to  provide  the  interfacial  boundary  conditions  (including  temperature,  density,  and 
species  concentration)  and  the  vaporization  rate  at  subcritical  conditions.  Match¬ 
ing  is  accomplished  by  requiring  the  prevalence  of  thermodynamic  phase  equilib¬ 
rium  and  continuities  of  mass  and  energy  fluxes,  as  shown  in  Figs.  2.3(a)-(d). 
Mass  Balance  for  the  Mixture 
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where  rh"  is  the  droplet  mciss  gasification  rate  per  unit  surface  area,  and  Ahv,i  the 
specific  enthalpy  of  vaporization  of  species  i.  The  subscripts  and  R-  represent 
conditions  at  the  droplet  surface  on  the  gas  and  liquid  sides,  respectively. 

When  the  droplet  surface  reaches  the  critical  mixing  point,  the  sharp  dis¬ 
tinction  between  gas  and  liquid  phases  disappears.  The  entire  flowfield  (including 
both  liquid  droplet  and  ambient  gases)  essentially  becomes  a  continuous  medium, 
with  no  abrupt  phase  transition.  Under  this  situation,  interfacial  matching  condi¬ 
tions  are  no  longer  needed,  and  the  entire  flow  field  can  be  treated  as  a  single-phase 
fluid. 

2.2  Droplet  Behavior  in  Forced-Convective  Environments 

The  second  case  deals  with  cryogenic  fluid  droplets  in  supercritical  streams. 
A  liquid  oxygen  (LOX)  droplet  initially  with  a  uniform  distribution  of  subcritical 
temperature  is  introduced  into  hydrogen  flows,  where  the  pressure  and  tempera¬ 
ture  are  both  in  the  supercritical  regime  of  oxygen,  as  shown  in  Fig.  2.1(b).  As 
a  result  of  heat  transfer  from  the  ambient  stream,  the  droplet  surface  reaches  the 
thermodynamic  critical  mixing  state  almost  instantaneously  upon  introduction  to 
the  hydrogen  gas  due  to  the  cryogenic  fluid  property  of  oxygen.  Once  this  oc¬ 
curs,  the  enthalpy  of  vaporization  and  surface  tension  vanish.  The  entire  flowfield 
(including  both  droplet  and  surrounding  gases)  becomes  essentially  a  continuous 
medium,  and  no  well-defined  liquid/vapor  interfacial  boundary  exists  as  for  a  sub- 
critical  droplet.  The  temperature  and  density,  as  well  as  their  spatial  gradients, 
vary  smoothly  throughout  the  domain  from  the  liquid  core  to  the  far  field.  In  the 
following,  a  single-phase  multidimensional  fluid  model  is  adopted  to  facilitate  the 
analysis. 


LIQUID  r,  VAPOR  w  LIQUID  .LT  VAPOR 
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Figure  2.3.  Matching  Conditions  at  Liquid/Vapor  Interface:  (a)  Mass  Fluxes  across  Droplet  Surface,  (b)  Mass  Fluxes  of 
Species  i  at  Liquid/Vapor  Interface,  (c)  Force  Balance  at  Droplet  Surface,  (d)  Energy  Transport  across  a  Liquid/Vapor 
Interface. 
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2.2.1  Governing  Equations 


The  flow  field  is  assumed  to  be  laminar  and  axisymmetric.  If  body  forces 
and  radiation  are  ignored,  the  governing  equations  can  be  written  in  the  following 
conservative  form. 


Mass: 


dp  dpu  dpv  ^ 

dt  dx  dr 


(2.18) 


Axial-momentum: 


dpu  d  ,  0  ,  dpuv 


^'^XX  ^'^TX  2d.  .  /O  1  0^ 

-fc  +  -  3 


Radial- momentum: 


dpu  dpuv  d  ,  2  \  dxxr  dxrr 

ir+-5r  +  +''>  =  -&+ir 

4  u  2  du  2  dv  2d, 


(2.20) 


Energy: 


2  3,  ,  23,  2^ 

-  33;<'‘“'’>  “  sa;*"'’ > 


(2.21) 


Species  Concentration: 


dpYi  dpuYi  dpvYi 

dt  dx  dr 


%  m,i  )x  .  %  m,t) 

dx  dr 


(2.22) 


The  physical  variables  t,  u,  v,  and  Yi  are  the  time,  velocity  components  in  the  axial 
and  radial  («-  and  r-)  directions,  respectively,  and  the  mass  fraction  of  component 
i.  The  specific  total  internal  energy  is  written  as 


20 


et  =  h-^  +  liu^  +  v^)  (2.23) 

P  2 

The  normal  stresses,  txx  and  T^r,  and  shear  stresses,  txt  and  Tn  are 


/4  5u  2  5u\ 

(2.24) 

Txx  “ 

^\3dx  3  dr) 

/du  dv\ 

(2.25) 

Txr  ~ 

/idv  2du\ 

(2.26) 

Trr  ~ 

^\3dr  3dx) 

where  fi  is  the  viscosity. 

Fourier’s  and  Pick’s  laws  are  applied  to  evaluate  the  thermal  and 

species 

diffusion  fluxes. 

{(ie)x  = 

dT  ^  dY 

1=1 

(2.27) 

{Qe)r  = 

,3r  *  ay; 

~^dr 

t=l 

(2.28) 

T, 

(2.29) 

•n  ^ 

pT^im 

(2.30) 

where  T,  A,  and  Vim  are  the  temperature,  thermal  conductivity,  and  mass  diffu- 
sivity  for  species  i  in  mixture,  respectively. 
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2.2.2  Boundary  Conditions 


Centerline  Boundary  Condition 


The  centerline  boundary  conditions  are  specified  based  on  flow  symmetry, 
which  are  given  by  setting  the  normal  velocity  and  gradients  of  the  pressure,  tan¬ 
gential  velocity,  temperature,  and  species  concentration  to  zero. 


V 


=  0 


r=0 


dp 

du 

_  dT 

dr 

dv 

dv 

r=0 

dv 

r=0 

(2.31) 

(2.32) 


Inflow  Boundary  Condition 


For  a  subsonic  inflow,  the  temperature,  species  concentration,  and  axial  and 
tangential  velocities  are  specified.  The  pressure  is  extrapolated  from  the  interior 
cells. 


Outflow  Boundary  Condition 

Since  the  outflow  is  subsonic,  a  back  pressure  is  specified,  with  other  flow 
variables  extrapolated  from  the  interior. 

2.3  Thermophysical  Properties 

One  major  challenge  in  the  present  analysis  is  the  establishment  of  a  unified 
property  evaluation  scheme  capable  of  treating  thermophysical  properties  over  the 
entire  mixture  thermodynamic  state  from  compressed  liquid  to  dilute  gas.  These 
properties  can  be  generally  divided  into  two  classes;  thermodynamic  (i.e.,  density, 
internal  energy,  and  specific  heat)  and  transport  (i.e.,  viscosity,  thermal  conduc¬ 
tivity,  and  mass  diffusivity)  properties.  At  high  pressure,  models  normally  used 
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to  represent  ideal-gas  behavior  may  encounter  significant  errors.  From  the  micro¬ 
scopic  point  of  view,  the  intermolecular  mean  free  paths  tend  to  decrease  with 
increasing  pressure;  the  molecular  volume  and  intermolecular  forces  are  no  longer 
negligible  as  for  the  ideal  fluids.  For  convenience,  each  property  is  expressed  as 
the  sum  of  an  ideal-gas  property  at  the  same  temperature  and  a  thermodynamic 
departure  function  which  takes  into  account  the  dense-fluid  correction.  This  de¬ 
parture  function  quantifies  thermodynamic  nonideality  and  can  be  derived  directly 
from  a  given  equation  of  state. 

For  a  multicomponent  system,  accounting  for  changes  in  properties  due  to 
mixing  is  much  more  complicated.  A  pseudo  pure  substance  model  is  adopted 
to  evaluate  the  properties  of  a  mixture,  treating  the  mixture  as  a  single-phase 
fluid  with  its  own  set  of  properties  evaluated  via  the  extended  corresponding-state 
(ECS)  principle.  This  method  improves  the  prediction  accuracy  and  only  limited 
data  (i.e.,  critical  properties  and  Pitzer’s  acentric  factors)  for  each  constituent 
component  are  needed.  Successful  application  of  the  corresponding- state  argument 
for  the  evaluation  of  fluid  p-V-T  properties  also  encourages  similar  improvement  in 
the  prediction  of  thermophysical  data.  Density  serves  as  an  important  parameter 
in  the  evaluation  process,  mainly  due  to  the  fact  that  the  temperature  dependence 
of  properties  at  constant  density  is  much  smaller  when  compared  with  that  at 
constant  pressure.  In  the  following,  the  corresponding-state  method  in  conjunction 
with  the  mixture  combining  rule  is  first  discussed,  followed  by  a  brief  summary  of 
the  equations  of  state  used  in  the  present  study.  Finally,  the  evaluation  schemes 
for  both  thermodynamic  and  transport  properties  are  addressed  thoroughly. 
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2.3.1  Extended  Corresponding-State  Principle 

Two  scaling  factors  are  used  to  characterize  the  conformal  mappings  of  tem¬ 
perature  and  density  between  the  mixture  of  interest  and  the  reference  fluid[21,22]. 

/m  =  |^  ;  am  =  -  (2.33) 

lo  P 

The  former  represents  the  conformation  of  temperature,  while  the  latter  character¬ 
izes  the  effect  of  mixture  molecular  size.  The  subscripts  m  and  0  refer  respectively 
to  the  fluid  of  interest  and  the  reference  fluid.  Assuming  all  components  in  a  mix¬ 
ture  obey  the  extended  corresponding-state  principle,  the  following  mixing  rules 
for  a  multicomponent  system  are  applied  to  the  scaling  factors 

u  =  E  E  Xi  X,  fi, 

j  =  l 

N  N 

(2-34) 

1=1  j=l 

where  N  and  Xi  are  the  number  of  species  and  the  mole  fraction  of  species  i.  The 
two  binary  scaling  factors  fij  and  hij  are  defined  as 

fij  — 

Su  =  +  (2-35) 

where  kij  and  £ij  are  binary  interaction  parameters  accounting  for  the  effects  of 
energy  and  molecular  size,  respectively.  In  order  to  model  the  quantum  behavior  of 
hydrogen,  interaction  parameters  for  a  hydrogen-containing  mixture  are  specially 
treated  by  a  generalized  correlation  proposed  by  Valderrama  and  Reyes  [23]. 

kij  =  A- BIT* 


24 


lij  =  l  +  Q.OnnliDMg/Dsmall)]  (2.36) 

where 


[Dbig/ D small]  —  [Vc.feij/Vc.sma//]  ^ 

A  =  0.1805  +  3.2  Ituj +2.437w| 

B  =  0.1323  +  0.5507a;j +3.5994u;|  (2.37) 

where  [Dbig/ D small]  is  always  greater  than  1  and  Vc  is  the  critical  molar  volume. 
The  variable  T*  is  the  reduced  temperature  and  u;  Pitzer’s  acentric  factor.  The 
subscripts  i  and  j  denote  the  hydrogen  and  non-hydrogen  components,  respectively. 

The  scaling  factors  (/*  and  hi)  for  each  individual  component  can  be  ob¬ 
tained  by  a  two-parameter  corresponding  state  principle[21,22]. 

/.  =  {Tc,^/Tc,o)0i{T*,vlu^) 

hi  =  (vc,t/vc,o)  (f>i{Ti,v*i,ui)  (2.38) 

where  Tc  is  the  critical  temperature  and  Oi  and  </)i  are  the  so-called  shape  factors 
which  are  functions  of  Pitzer’s  acentric  factor  (wj)  and  reduced  temperature  and 
molar  volume  {T*  and  v*).  The  analytical  expressions  for  Oi  and  <f>i  proposed  by 
Leach  and  Leland[24]  take  the  forms 

ei{T’,viui)  =  i  +  (ui,-wo)f(7;*,v?) 

=  [l  +  (wi-«o)G(i;*,vn]Z«,o/Ze,i  (2.39) 

and 


F{T*,V^)  =  ai  +  hlnTt{ci+di/Tt){yt-0.5) 
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G{T*,v:)  =  a2(v+  +  62)  +  C2(vt  +  rf2)lni;+  (2.40) 

Tt  =  min{2,max{T*,0.5}} 

v+  =  mm{2,  maa;{v*,  0.5}}  (2-41) 

where  Z  is  the  compressibility  factor,  and  the  asterisk  *  represents  reduced  prop¬ 
erties  defined  as  the  values  which  are  divided  by  the  corresponding  quantity  at 
the  critical  point.  Constants  a,',  bi,  ci,  and  d,  are  coefficients  for  shape-factor 
correlations  and  can  be  found  in  Ref.  [24]. 

2.3.2  Equations  of  State 

An  adequate  equation  of  state  is  required  to  complete  the  formulation  for 
the  property  evaluation.  In  the  following,  two  commonly  used  equations  of  state 
are  described:  the  Soave-Redlich-Kwong  (SRK)  and  Benedict- Webb- Rubin  (BWR) 
equations  of  state. 

2.3. 2.1  Soave-Redlich-Kwong  Equation  of  State 

The  Soave-Redlich-Kwong  equation  of  state  is  an  extension  of  the  original 
Redlich-Kwong  equation  of  state  by  introducing  an  acentric  factor  and  a  tem¬ 
perature  dependence  into  the  cohesive  energy  term  to  account  for  the  effect  of 
non-sphericity  on  fluid  p-V-T  properties.  It  takes  the  form 

_  ‘RuT  oca 
^  Y  —  b  v(v  4-  6) 

where  Ru  and  v  are  the  universal  gas  constant  and  molar  volume,  respectively. 

The  mass-averaged  parameters  aa  and  b  for  a  mixture  are  determined  using  the 

van  der  Waals  one-fluid  combining  rule[25]. 
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N  N  N 

aa  =  X!  E  WjCirja^,,  b  =  X  (2.43) 

2=1  7  =  1  2=1 

where  Xi  is  the  mole  fraction  of  species  e,  and  the  cross  mixture  parameter  aijaij 
for  the  ij  pair  is 


OLijdij  =  (1  -  kij)y/aiajaiaj  (2.44) 

where  kij  is  the  interaction  coefficient  which  accounts  for  the  deviation  from  the 
geometric-mean  combining  rule  for  aa,  and  a,-  is  an  empirical  function  of  reduced 
temperature  and  Pitzer’s  acentric  factor.  For  most  of  hydrocarbon  fuels  and  certain 
non- hydrocarbons,  ai  can  be  expressed  as 


oi  =  [1  +  Siil  -  ^^)]2  (2.45) 

with  T*  the  reduced  temperature.  The  expression  for  Si  suggested  by  Graboski 
and  Daubert[26]  is  given  below. 


Si  =  0.48508  -f  1.55171u»,-  -  0.15613wf  (2.46) 

where  Ui  is  Pitzer’s  acentric  factor  of  component  i.  In  Eqs.  (2.43)  and  (2.44),  the 
energy  and  volume  constants  ai  and  hi  depend  on  the  critical  properties  of  each 
species  and  take  the  form 


(Xi 


hi 


0.42747  ’ 

Pc,i 

'R'uTcA 


0.08664 


Pc,i 


(2.47) 


with  Tc,i  and  pc,i  the  critical  temperature  and  pressure  of  component  i,  respectively. 
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If  the  compressibility  factor  is  introduced  to  characterize  the  thermody¬ 
namic  non-ideality,  the  equation  of  state  is  rewritten  as 

pv  =  ZTluT  (2.48) 

Substitution  of  (2.48)  into  (2.42)  and  rearrangement  of  the  result  lead  to  a  cubic 
polynomial  equation  for  the  compressibility  factor  Z. 

Z^-Z^  +  {A-B-  B^)Z  -AB  =  0  (2.49) 


and 


c  =  A-B-B^-\ 

O 

d  =  ^  p{A-B-B^)-  21 AB  -  2  (2.53) 


For  a  system  involving  two  phases,  at  a  given  temperature  and  pressure,  the  small¬ 
est  real  root  is  designated  to  the  liquid  phase,  while  the  largest  to  the  gas  phase. 
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2.3. 2. 2  Benedict- Webb-Rubin  Equation  of  State 

Although  the  Benedict- Webb-Rubin  (BWR)  equation  of  state  is  more  com¬ 
plicated  than  classical  cubic  equations  of  state  (e.g.,  the  SRK  equation  of  state), 
it  has  demonstrated  superior  performance  over  wide  ranges  of  temperature  and 
pressure.  However,  the  constants  for  the  BWR  equation  of  state  are  valid  only 
for  a  limited  number  of  pure  compounds [27-29].  To  overcome  this  constraint,  an 
extended  corresponding-state  (ECS)  principle  is  used.  The  basic  idea  is  to  assume 
the  properties  of  a  single  phase  fluid  can  be  evaluated  via  conformal  mappings  of 
temperature  and  density  to  those  of  a  given  reference  fluid.  As  a  result,  only  the 
BWR  constants  for  the  reference  substance  are  needed. 

Under  the  assumption  of  the  ECS  principle,  the  density  of  a  mixture  can 
be  evaluated  by 


p{T,p,Xi)^ 


po{To,po) 

hm{T,p,Xi) 


(2.54) 


where  hm  is  the  scaling  factor  and  po,  To,  and  po  are  the  density,  temperature  and 
pressure,  respectively,  of  the  reference  fluid  at  the  corresponding  state  defined  by 
Eq.  (2.33).  The  related  temperature  at  the  conformal  state  takes  the  form 


To  = - - - r  (2.55) 

"  fm{T,p,Xi)  ^  ’ 

Based  on  the  general  compressibility  theory [30],  the  corresponding  pressure  for 
reference  fluid  is  evaluated  by 


TABLE  2.1 

BWR  Parameters  as  Functions  of  Temperature. 
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a,(T) 

TIT 

ciT) 

= 

NiT  +  +  iVa  +  N4/T  +  N^/T^ 

a,(T) 

= 

NeT  +  Nr  +  Ns/T  +  N^/T^ 

MT) 

= 

NioT  +  iVn  +  Nn/T 

as(T) 

iVl3 

«6(r) 

= 

iVl4/T  +  iVi5/T2 

C7(T) 

N,efT 

a,(T) 

N^r/T  +  N,s/T^ 

as(T) 

ciio(T) 

iVao/r^  +  iVai/T" 

«ii(T) 

== 

N22/T^  +  N2z/T^ 

ai2(r) 

=: 

N24fT^  +  N25fT^ 

ai3(T) 

N2e/T^  +  N2r/T* 

GuiT) 

=z 

N2s/T^  +  N29/T^ 

G15(T) 

=r 

Nzo/T^  +  NsilT^  +  Nz2lT^ 

MTo,  po)  =  E  +  E  <*»(ro)/>S”“‘"e-"o  (2.57) 

n=l  n=10 

where  7  is  0.04  and  the  units  for  pressure,  density,  and  temperature  are 

pressure  (p);  bar 

density  (p):  liter/mole 

temperature  (T):  Kelvin 


The  functional  forms  of  ai{T)  are  given  in  Table  2.1,  where  the  coefficients 
Ni  (i  =  1...32)  are  dependent  on  the  reference  fluid.  In  the  current  study,  methane 
is  chosen  as  the  reference  fluid  because  of  the  availability  of  experimental  data 
for  both  volumetric  and  transport  properties  over  a  broad  range  of  fluid  state. 
Table  2.2  provides  the  temperature  constants  of  methane  for  the  BWR  equation 


of  state. 


TABLE  2.2 

Temperature  Constants  of  Methane  for  the  BWR  Equation  of  State. 
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Ni  =  -1.184347314485x10-2 

N2  =  7.540377272657x10-^ 

N3  =  -1.225769717554x10+^ 

^4  =  6.260681393432x10+2 

N5  =  -3.490654409121xl0+'‘ 

Ne  =  5.301046385532x10-'* 

N7  =  -2.875764479978x10-^ 

Ns  =  5.011947936427x10+^ 

Ns  =  -2.821562800903x10+^ 

iVio  =  -2.064957753744x10+5 

Nn  =  1.285951844828x10-2 

Ni2  =  -1.106266656726x10+° 
Ni3  =  3.060813353408x10-'* 

Nu  =  -3.174982181302x10-5 
iVi5  =  5.191608004779x10+° 

iVi6  =  -3.074944210271x10-'* 


iVi7  =  1.071143181503x10-5 

Nis  =  -9.290851745353x10-5 
Nis  =  1.610140169312x10-'* 

N20  =  3.469830970789x10+'* 

N21  =  -1.370878559048x10+° 
N22  =  1.790105676252x10+2 

N23  =  1.615880743238x10+° 

N24  =  6.265306650288x10-1 

N25  =  1.820173769533x10+1 

N26  =  1.449888505811x10-5 

N27  =  -3.159999123798x10+1 
N2S  =  -5.290335668451x10-° 
N29  =  1.694350244152x10-5 

^30  =  8.612049038886x10-° 

N31  =  -2.598235689063x10-° 
N32  =  3.153374374912x10-5 


Although  this  equation  of  state  must  be  solved  iteratively  for  density  at 
a  given  pressure  and  temperature,  the  prediction  covers  a  wide  range  of  thermo¬ 
dynamic  state,  and  as  such  promotes  the  establishment  of  a  unified  evaluation 
scheme  of  thermophysical  properties.  Figure  2.4  shows  the  error  of  density  predic¬ 
tion  for  oxygen  at  supercritical  pressures  of  interest  (100  and  400  atm)  by  three 
commonly  used  equations  of  state,  i.e.,  the  Benedict- Webb-Rubin  (BWR)  equation 
of  state  with  extended  corresponding-state  principle  (ECS),  Peng-Robinson  (PR), 
and  Soave-Redlich-Kwong  (SRK)  equations  of  state.  The  BWR  equation  of  state 
together  with  the  ECS  principle  shows  the  deviation  as  low  as  1.5%  for  both  cases 
(100  and  400  atm).  However,  the  SRK  and  PR  equations  of  state  may  miscalculate 
the  density  by  13%  and  17%,  respectively.  Figure  2.5  shows  the  comparison  of  oxy¬ 
gen  density  between  experimental  data  and  the  prediction  by  the  BWR  equation 
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Figure  2.4.  Relative  Errors  of  Density  Predictions  by  Three  Different  Equations 
of  State. 


of  state  in  conjunction  with  the  ECS  principle.  Excellent  agreement  is  obtained 
over  a  wide  range  of  fluid  states,  from  compressed  liquid  to  dilute  gas. 

2.3.3  Thermodynamic  Properties 

Thermodynamic  properties  are  very  important  quantities  in  the  analyses  of 
supercritical  droplet  behavior.  Variations  in  these  properties  can  often  be  related 
to  local  fluid  properties  (i.e.,  pressure,  temperature,  and  species  concentration). 
Properties  like  enthalpy,  internal  energy,  and  specific  heat  of  a  mixture  are  ex¬ 
pressible  as  the  sum  of  the  ideal-gas  property  and  a  correction  term  accounting 
for  the  thermodynamic  non-idealities.  A  detailed  derivation  can  be  found  in  Ap¬ 
pendix  A. 
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Temperature,  K 


Figure  2.5.  Comparison  of  Density  Data  Predicted  by  the  BWR  Equation  of 
State  and  Experimental  Measurement. 
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N 


EK 

Z  =  1 


^Tef,i 


(2.61) 


N 

- 

i=l 


°ref,i 


+ 


(2.62) 


1  =  1 


0 

p,i 


(2.63) 


where  Yi  is  the  mass  fraction  of  component  and  h^gf  i,  ^ref,v  ^p,i 
ideal-gas  enthalpy,  internal  energy,  and  specific  heat  at  a  reference  temperature  of 
Tref,  respectively. 


2.3.4  Itansport  Properties 

Ely  and  Hanley  proposed  a  corresponding-state  model  which  is  capable 
of  predicting  the  transport  properties  of  fluid  mixtures  over  the  entire  p-V-T 
states[21,32].  This  method  requires  only  the  common  characterization  parame¬ 
ters  (i.e.,  critical  properties  and  Pitzer’s  acentric  factor)  of  each  component  as 
input  parameters. 


Viscosity 

The  mixture  viscosity  can  be  evaluated  via  the  corresponding-state  argu¬ 
ment  as 

n{p,T,Xi,Mi)  =  iJ,o{pQ,To)J^fi{p,T,  Xi,  Mi)  (2.64) 

where  po  is  the  value  of  the  reference  fluid  at  the  corresponding  thermodynamic 
state.  The  scaling  factor  Tfi  is  defined  herein  as 
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T  —  f— 


(2.65) 


where  Mo  is  the  molecular  weight  of  the  reference  fluid,  and  the  value  of  the 


mixture  denoted  by  M  can  be  evaluated  by 


z=l  j=l 


(2.66) 


The  scaling  factors  fm,  fij,  and  kij  follow  the  same  definitions  in  Eqs.  (2.33) 
and  (2.35).  A  binary  mixing  rule  is  used  for  the  molecular  weight. 


Mij  =  2MiMj/{Mi  +  Mj) 


(2.67) 


where  Mi  is  the  molecular  weight  of  component  i.  The  viscosity  for  the  reference 
fluid  can  be  written  as [21] 


//o(po,2o)  =  Ho  (^o)  +  Mo  (^o)mo  +  ^ho{po,To)Xh 


(2.68) 


where  Ho  Mo  represent  the  viscosity  at  the  dilute-gas  limit  and  the  first  den¬ 

sity  correction,  respectively.  The  third  term  is  the  high  density  correction,  with 
Amo  being  a  dense-fluid  remainder.  These  three  terms  can  be  evaluated  based  on 
the  formulas  listed  in  Table  2.3.  A  small  correction  factor,  is  needed  to  ac¬ 
count  for  possible  errors  induced  by  the  corresponding  state  model.  The  following 
expression  based  on  the  modified  Enskog  theory  is  proposed  by  Hanley [33]. 


Zc.0  I  ^ 


(2.69) 


where  the  subscript  c  refers  to  the  critical  values.  The  correction  factor  rp,  which 


empirically  accounts  for  the  size  effect,  is  selected  to  be  1.5. 


TABLE  2.3 

Correlation  of  Viscosity  for  Methane;  Units  are  gram-cm-sec-Kelvin. 


35 


O 

o 

< 

+ 


+ 


o 

o 

3. 


I 

5£ 


u 


q; 

(h 

(U 


to 

CO 

CO 

iO 

CO 

c<» 

rH 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1 

O 

o 

O 

o 

o 

o 

o 

o 

o 

1-H 

t—H 

r-H 

t-H 

T— H 

T-— < 

t-H 

t-H 

t-H 

X 

X 

X 

X 

X 

X 

X 

X 

X 

r- 

CO 

00 

H 

o 

o 

o 

o 

c:^ 

o 

CO 

CO 

CO 

lO 

o 

t-H 

O 

CO 

o 

oo 

oo 

CO 

CM 

rH 

t-H 

o 

t-H 

Ci 

1-H 

o 

o 

00 

CM 

t-H 

r-H 

d 

1 

r— 

00 

r-H 

CO 

t-H 

Oi 

t-H 

CQ 

00 

r-H 

CM 

CO 

lO 

iO 

l>- 

'Z^' 

o 

1-H 

o 

CO 

CO 

t-H 

CM 

CO 

cs 

CO 

<N 

Oi 

cj 

CO 

CO 

r*H 

CO 

o 

]>: 

t-H 

1>^ 

CO 

T— H 

o 

CM 

1 

1 

1 

i 

C 

4- 

o 

r-t 

(N 

+ 

CO 

+ 

O 

1 

o 

r— H 

H- 

o 

T-H 

t-H 

to 

X 

X 

X 

C3 

T— 1 

00 

oo 

H 

+ 

t- 

(M 

o 

CO 

CO 

»o 

Ci 

iO 

CO 

TT  ' 

00 

CM 

ir> 

6 

o 

CO 

Cv» 

aL 

Oi 

CO 

IS 

CO 

CO 

+ 

1 

r— H 

r— 1 

I 

tj 

1 

a. 

a> 

_ 

<N 

CO 

CO 

aL 

+ 

+ 

+ 

+ 

1 

+ 

+ 

<0 

'' — ^ 

o 

o 

o 

O 

o 

o 

O 

H 

r-H 

1— H 

T— H 

t-H 

T— H 

II 

+ 

X 

X 

X 

X 

X 

X 

X 

!>- 

Oi 

CM 

o 

CS) 

CO 

00 

Oi 

CO 

CM 

Oi 

CO 

CM 

T— H 

t- 

t- 

o 

CM 

iO 

00 

CM 

CTi 

Oi 

CO 

CM 

CM 

o 

CO 

CO 

t— H 

00 

lO 

CO 

lO 

CM 

CO 

OO 

CM 

o 

CO 

00 

t- 

o 

Q 

CO 

CM 

CO 

CO 

o 

CM 

C<I 

CO 

CM 

o 

o 

00 

CO 

CD 

<I 

t-H 

1 

T— H 

t-H 

CM 

1 

r-H 

t-H 

l6 

•  <s» 

1-H 

CM 

CO 

lO 

CO 

t- 

oo 

36 


Thermal  Conductivity 

Evaluation  of  thermal  conductivity  must  be  carried  out  carefully  for  two 
reasons:  (1)  the  one-fluid  model  must  ignore  the  contribution  of  diffusion  to  con¬ 
ductivity,  and  (2)  the  effect  of  internal  degrees  of  freedom  on  thermal  conductivity 
cannot  be  correctly  taken  into  account  by  the  corresponding-state  argument.  As 
a  result,  thermal  conductivity  of  a  pure  substance  or  mixture  is  generally  divided 
into  two  contributions  [22], 

A(p,  r,  X.,  MO  =  A'^(r,  Xi,  Mi)  +  xlip,  t,  x,,  mo  (2.70) 

The  former  arises  from  the  energy  transfer  via  the  internal  degree  of  freedom,  while 
the  latter  is  due  to  the  effects  of  molecular  collision  or  translation.  For  a  mixture, 
can  be  evaluated  via  a  semi-empirical  mixing  rule. 

AU  =  E  i;  (2.71) 

i=l  j=l 

where  A^-  is  binary  thermal  conductivity  and  is  deflned  as 

(a;,)''=2[(A;)-'+(A')-‘]  (2.72) 

The  internal  contribution  A[-  is  calculated  using  the  modified  Eucken  correlation  [34] 
for  polyatomic  gases  as 

A'  =  1.32^(C»i  -  jK)  (2.73) 

where  the  properties  Mi,  and  are  the  d  ite  gas  viscosity,  molecular  weight, 
and  ideal-gas  heat  capacity  of  component  i,  respectively. 

In  Eq.  (2.70),  the  collisional  and  translational  contribution  X'ln{T,p)  can  be 
evaluated  via  the  extended  corresponding-state  method. 
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C(T,p)  =  >.’i{p«,To)TxX^ 


(2.74) 


Following  the  corresponding-state  principle,  the  scaling  factor  J-x  is  defined  herein 


(2.75) 


where  the  pseudo  molecular  weight  is 


N  N 


Y,  T,  ff 

^=1  j=l 


(2.76) 


The  correction  for  non- correspondence  is  applied  over  the  entire  collisional  or  trans¬ 


lational  contribution  to  the  thermal  conductivity  and  is  defined  as 


T  fdfrr 


1- 


2  =  1 


(2.77) 


This  expression  is  derived  from  the  Enskog  theory  so  that  the  density  dependence  of 
the  thermal  conductivity  is  simply  depended  on  the  derivative  term  of  {dfmIdT),,. 
The  value  \q{po,Tq)  is  the  equivalent  contribution  for  the  reference  fluid  at  the 
corresponding  state.  Based  on  Hanley’s  model[22],  this  reference  value  can  be 
divided  into  three  contributors. 


A'o'(po,To)  =  AW(ro)  +  \^o\To)po  +  AAo(po,To) 


(2.78) 


where  aI^^  is  the  value  at  the  dilute-gas  limit  and  can  be  evaluated  by 
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157?.^  Q 


(2.79) 


where  (jLq  is  the  ideal  gas  viscosity.  The  quantities,  AAq  ^  and  AAq  are  the  first  den¬ 
sity  and  high  dense  corrections.  The  mathematical  formulas  of  these  correlations 
for  methane  can  be  found  in  Table  2.4. 

Modern  theory  of  transport  phenomena  predicts  an  infinite  thermal  con¬ 
ductivity  at  the  pure-fluid  critical  point  and  a  large  thermal  conductivity  in  the 
vicinity  of  the  critical  point [35,36].  Although  this  near-critical  enhancement  has 
been  observed  experimentally  for  several  fluids,  the  theory  of  non-classical  critical 
behavior  predicts  that  the  difference  among  critical  criteria  for  different  compo¬ 
nents  may  diminish  the  significant  of  the  critical  anomaly [35,37, 38].  Because  the 
effect  of  the  critical  point  contribution  for  a  mixture  is  not  well  defined  and  is  likely 
to  be  small,  it  is  not  included  in  the  current  analysis  of  droplet  vaporization. 


Mass  DifFusivity  of  Gas  Phase 

Estimation  of  the  binary  mass  diffusivity  for  a  gaseous  mixture  at  high 
density  represents  a  more  challenging  task  than  evaluating  the  other  transport 
properties [39],  due  to  the  lack  of  a  formal  theory  or  even  a  theoretically  based 
correlation.  For  a  multicomponent  mixture,  the  effective  diffusion  coefficient  Vim 
for  each  species  i  can  be  related  to  the  binary  diffusion  coefficients  through  the 
following  formula: 


TABLE  2.4 

Correlation  of  Thermal  Conductivity  for  Methane;  Units  are  gram-cm-sec-Kelvin. 
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(2.80) 


where  Xi  is  the  mole  fraction  of  species  i.  The  binary  mass  diifusivity  'Dij  for 
the  gas  phase  is  obtained  using  the  corresponding-state  model  proposed  by  Taka- 
hashi[40].  The  approach  proceeds  in  two  steps.  First,  the  binary  mass  diffusivity 
of  a  dilute  mixture  is  obtained  using  the  Chapman-Enskog  theory  in  conjunction 
with  the  Lennard- Jones  intermolecular  potential-energy  function[41]. 


Vij  =  0.0018583 


-F  IjMj) 

PO-'ij^D 


0.5 


(2.81) 


and 


1.06036  0.193  1.03587  1.76474 

D  -  (j^*)o.i56i  +  exp(0.47635  T*)  ea:p(1.52996  T*)  exp(3.89411  T*y  ‘  ’ 

rr.*  kT 

T*  =  — 

6ij 

tij  = 

(^ij  =  +  (2-83) 


where  k  is  the  Boltzmann  constant,  and  e  and  a  are  the  Lennard-Jones  parameters 
which  quantify  respectively  the  energy  and  length  scale  of  gas  molecules.  These 
values  can  be  found  in  Appendix  B  of  Ref.  [34]  for  a  variety  of  species.  The 
calculated  data  for  a  dilute  mixture  is  then  corrected  with  a  pressure  correction 
which  is  expressed  in  a  generalized  chart  in  terms  of  reduced  temperature  and 
pressure,  as  shown  in  Fig.  2.6.  This  scheme  appears  to  be  the  most  complete 
to  date,  and  has  demonstrated  moderate  success  in  the  limited  number  of  tests 
conducted. 


For  the  liquid  phase,  the  mass  difFusivity  is  evaluated  using  a  model  pro¬ 


posed  by  Perkins  and  Geankoplis[42]. 


=  E  (2.84) 

where  Hm  and  fij  are  the  viscosities  for  the  mixture  and  component  j,  respectively. 
The  binary  diffusivity  D,y  is  evaluated  using  the  Hayduk  and  Minhas  model  [43]. 

Vij  =  13.3  *  10~^— pp  (2.85) 

^  li 

where 

e 


(10.2/vi)  -  0.791 
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fij  =  viscosity  of  solvent  jt,  cp. 

Vi  =  molar  volume  of  solute  i  at  its  normal  boiling  temperature,  cm^  fmole. 
T  =  temperature,  K. 

2.4  Liquid/Vapor  Phase  Equilibrium 


Prior  to  occurrence  of  the  critical  mixing  state,  thermodynamic  vapor-liquid 
phase  equilibrium  is  assumed  to  prevail  at  the  droplet  surface.  This  condition 
first  requires  that  thermal  equilibrium  be  achieved,  with  the  temperatures  of  both 
phases  being  identical. 

where  the  superscripts  v  and  £  represent  the  values  in  the  vapor  and  liquid  phases, 
respectively.  It  is  also  evident  that  no  unbalanced  mechanical  forces  exist  between 
the  two  systems,  and  as  such  the  second  condition  for  phase  equilibrium  is  that 
the  pressures  of  both  phases  remain  equal. 

p'  =  p'  (2.87) 

For  a  multicomponent  system,  the  condition  of  thermodynamic  equilibrium  is  en¬ 
forced  by  requiring  that  the  chemical  potential  of  each  component  be  the  same  in 
both  phases  at  a  given  pressure  and  temperature. 


(2.86) 


(2.88) 


If  this  condition  does  not  hold,  there  would  be  a  tendency  for  mass  to  pass  from 
one  phase  to  the  other [30].  Since  the  partial  molar  Gibbs  function  is  equivalent  to 
the  chemical  potential,  Eq.  (2.88)  can  be  rewritten  as 
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gj  =  gl  i  =  l,2,...,iV 


(2.89) 


Following  the  definition  of  fugacity  of  a  component  in  a  mixture 


7^rd(ln/j)j’  =  dgi 
Ylm(^fi/X^p^  =  1 

Equation  (2.89)  can  be  expressed  in  terms  of  fugacity  as 


(2.90) 

(2.91) 


ln/?  =  ln/f  .■=1,2,... ,iV 


(2.92) 


The  fugacity  fi  is  related  to  the  volumetric  properties  of  the  mixture  through  the 


following  equation. 


(2.93) 


where  p  is  the  density  of  the  system,  and  n,-  the  number  of  mole  of  component  i. 
Note  that  /,  depends  explicitly  on  pressure,  temperature,  and  mixture  composition. 
As  a  result,  an  adequate  equation  of  state  is  required  to  complete  the  formulation 
for  the  phcise-equilibrium  analysis. 

Using  Soave-Redlich-Kwong  equation  of  state,  Eq.  (2.93)  can  be  integrated 
to  yield  an  explicit  expression  for  the  fugacity  of  each  species  i. 


HfilpXi)  =  -I(Z  -  l)  -  ln(z  -  r)  -  - 


p2 

A  j^i _  ^ 

B  aa  h 


-i  '"('+?) 


where  A  and  B  are  defined  in  Eq.  (2.50).  Equation  (2.92)  can  be  rewritten  as 


k/f  =  ln/f  .'=1,2,. ..,iV 


(2.95) 


where 
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In/f  =  lnXf +lnp  + 


In  //  =  In  Xf  +  In  j9  + 


In  order  to  close  the  problem,  two  more  equations  are  given 

N 

-  1 
i=i 

'txj  =  1  (2.98) 

j=i 

For  a  binary  system  {N  =  2)  at  a  given  pressure  and  temperature,  there  are  four 
unknowns  involved  in  solving  the  phase  equilibrium  conditions.  These  unknowns 
include  the  mole  fraction  of  each  component  in  the  gas  phase  (Xf  and  Xf)  and 
those  of  liquid  phase  (Xf  and  Xl).  Therefore,  four  independent  equations  are 
required  to  solve  the  unknowns. 

To  solve  these  mutually  coupled  equations,  an  efficient  numerical  algorithm 
using  the  Newton-Raphson  iteration  technique  is  adopted.  Each  unknown  is  up¬ 
dated  by  its  respective  correction  during  the  iteration.  At  (n  -(-  1)*^  level,  those 
unknowns  are  expressed  as 


A 

B 


2  XJaijaij  , 


N 


A 

B 


p2  ^2  X^OLijd 


3=\ 


aa 


(2.97) 


(2.99) 


Substitute  Eqs.  (2.99)-(2.102)  into  (2.98)  and  (2.95)-(2.97),  and  the  corrections  for 
the  (n  +  1)**  iteration  level  can  be  solved  using  the  following  vector  form 


where 
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dA 


JL  JL  f)X; 

+  JlYl  -  kij)^otiaiajaj 

i=lj=l 

N  N 

—  2  6ii;Xj(^l  kij)yJaiaiQjaj^ 

1=1  j=\ 

p  daa 


dXk  niT^  dx. 


dB 


p  db 


dXk  KuTdXk 


dZ 


i-Z  +  B)+-^{Z  +  2BZ  +  A) 


dXk 


dXk 


dXk  SZ^-2Z  +  A-B-  B^  (2.103) 

To  facilitate  the  numerical  analysis,  initial  values  of  X^,  Xl,  X\  and  are 
required.  At  a  given  XI,  the  mole  fraction  of  the  other  component  is  determined. 


XI  =  \-  XI 


(2.104) 


Application  of  the  Raoult’s  law  leads  to 


jf'  =  ^i’4 

p\ 


X‘  =  (2.105) 

P2 

where  pf  is  the  saturation  pressure  of  component  i. 

2.4.1  Enthalpy  of  Vaporization 

Enthalpy  of  vaporization  is  defined  as  the  heat  absorbed  by  a  component 
to  transform  from  its  liquid  phase  to  vapor  at  equilibrium  temperature  and  corre¬ 
sponding  vapor  pressure.  It  can  be  evaluated  by 
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AK,i  =  h\-h\ 


(2.106) 


and 


h^  =  h^i-nT^—[HfiipXi)\ 


(2.107) 

where  superscript  0  denotes  quantity  in  the  ideal  state.  The  enthalpy  of  vaporiza¬ 
tion  for  component  i  is 


=  TIT^ 


(2.108) 


dT  dT 

The  explicit  form  of  Eq.  (2.108)  based  on  SRK  equation  of  state  is  given  in  Ap¬ 
pendix  B. 


2.4.2  Results  and  Discussion 

In  1875,  Gibbs  proposed  a  phase  rule[30]  which  takes  the  form 

V  =  -P  +  C4-2  (2.109) 

where  V  and  C  are  the  numbers  of  phases  and  components  in  the  system  of  interest, 
respectively,  and  the  quantity  V  is  the  variance  which  designates  the  number  of 
intensive  properties  that  must  be  specified  to  fix  the  state  of  the  system.  For 
example,  a  two-phase  mixture  of  n-pentane  [CzH^)  and  nitrogen  (iV2),  the  number 
of  phases  V  equals  2,  the  number  of  component  C  is  2,  and  therefore  the  value  for 
the  variance  V  becomes  2.  This  means  that  once  any  two  of  the  intensive  properties 
are  given(i.e.,  pressure,  temperature,  and  mole  fraction  of  a  given  component  in 
any  phase,  etc.),  the  state  of  the  system  is  fixed,  and  all  other  properties  can  be 
determined  using  thermodynamic  relationships. 
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Figure  2.7.  Phase  Equilibrium  Compositions  for  O2IH2  System  at  Various  Pres¬ 
sures. 


Figure  2.7  shows  the  equilibrium  compositions  for  H2IO2  system  in  the 
pressure  range  of  5-250  atm,  where  pr  stands  for  the  reduced  pressure  of  oxygen. 
In  the  subcritical  regime  {pr  <  1),  the  liquid  phase  consists  almost  entirely  of 
oxygen.  The  amount  of  hydrogen  dissolved  in  the  liquid  oxygen  is  quite  limited, 
decreasing  progressively  with  increasing  temperature  and  dropping  to  zero  at  the 
boiling  point.  However,  at  supercritical  pressures  (pr  >  1),  the  hydrogen  gas 
solubility  becomes  substantial  and  increases  with  temperature.  For  example,  the 
liquid  phase  contains  8.2-45  %  hydrogen  (by  mole  fraction)  for  pr  =  5,  depending  on 
temperature.  In  spite  of  its  high  molar  concentration,  the  mass  fraction  of  dissolved 
hydrogen  is  much  smaller  because  of  the  disparity  in  molecular  weight  between 
hydrogen  and  oxygen.  The  mass  fraction  of  hydrogen  in  the  liquid  phase  is  only  5 
%  at  the  critical  mixing  point  for  pr  =  5,  as  opposed  to  the  mole  fraction  of  45  %. 
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Figure  2.8.  Pressure- Temperature  Diagram  for  Phase  Behavior  of  O2I System 
in  Equilibrium. 

The  phase  equilibrium  results  also  indicate  that  the  critical  mixing  temperature 
decreases  with  pressure.  Consequently,  a  droplet  can  reach  its  critical  condition 
more  easily  at  a  higher  pressure. 

The  overall  phase-behavior  in  equilibrium  is  best  summarized  by  the  pressure- 
temperature  diagram  presented  in  Fig.  2.8,  showing  how  the  phase  transition  oc¬ 
curs  under  different  thermodynamic  conditions.  The  boiling  line  is  made  up  of 
boiling  points  for  subcritical  pressures.  As  the  temperature  increases,  an  equi¬ 
librium  vapor-liquid  mixture  may  transit  to  super-heated  vapor  across  this  line. 
The  critical  mixing  line  registers  the  variation  of  the  critical  mixing  temperature 
with  pressure,  and  intersects  the  boiling  line  at  the  critical  point  of  pure  oxygen, 
the  highest  temperature  at  which  the  vapor  and  liquid  phases  of  an  O2/H2  binary 
system  can  coexist  in  equilibrium. 
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Figure  2.9.  Effect  of  Pressure  on  Enthalpy  of  Vaporization  of  O2  in  an  Equilibrium 
Mixture  of  O2  and  /?2- 

Figure  2.9  shows  the  variations  of  the  enthalpy  of  vaporization  of  oxygen  in 
an  equilibrium  mixture  of  hydrogen  and  oxygen.  In  contrast  to  the  conventional 
approach  for  low-pressure  applications  where  the  ambient  ga.s  solubility  is  ignored 
and  the  heat  of  vaporization  depends  on  temperature  only,  the  enthalpy  of  va¬ 
porization  is  a  strong  function  of  both  temperature  and  pressure,  decreasing  with 
increasing  pressure  and  dropping  rapidly  to  zero  at  the  critical  mixing  point. 

2.5  Chemical  Reactions 

Finite-rate  chemical  reaction  model  is  employed  in  the  present  study.  When 
the  droplet  is  introduced  into  a  high-pressure  environment,  the  gasification  rate 
increases  due  to  diminishing  latent  heat  of  vaporization.  The  transport  processes 
in  both  the  liquid  and  gas  become  compatible  with  each  other.  As  a  result,  the 
mixing  of  fuels  and  oxidizers  may  occur  in  a  time  scale  commensurate  with  the 
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chemical-reaction  time.  The  use  of  a  diffusion-controlled  complete  reaction  (flame- 
sheet)  model  may  lead  to  erroneous  results  for  predicting  the  droplet  burning 
mechanisms. 

For  combustion  of  gaseous  hydrocarbon  fuels  with  oxygen,  a  single-step 
global  reaction  mechanism  suggested  by  Westbrook  and  Dryer[44]  is  adopted.  This 
model  is  a  simplification  of  reality,  but  provides  reasonable  predictions  of  chemical 
reaction  rates  over  a  wide  range  of  equivalent  ratio.  The  reaction  step  takes  the 
form 

CnHm  +  (n  +  j)02  ^  nC02  +  (2.110) 

The  rates  of  production  of  species  are 


<^CnHm  = 

(2.111) 

u;o2  = 

(2.112) 

WC’02  = 

/  Mc02  \  ■ 

(2.113) 

^H20  = 

(m\  (  MH2O  \  . 

(2.114) 

where  M  and  [C\  are  the  molecular  weight  and  species  concentration.  The  acti¬ 
vation  energy  Ea  and  factors  a  and  5,  are  dependent  on  the  type  of  fuel  species 
considered  and  their  values  are  listed  in  Table  2.5. 


TABLE  2.5 

Reaction  Rate  Parameters  Proposed  by  Westbrook  and  Dryer  in  1981;  Units  are  cm-sec-mole-Kcal- Kelvin 
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CHAPTER  3 


NUMERICAL  SCHEMES 


The  present  analysis  of  droplet  vaporization  and  combustion  poses  several 
difficulties  to  contemporary  numerical  algorithms.  First,  the  model  of  interest  is 
a  subcritical  liquid  droplet  being  introduced  into  a  supercritical  gaseous  environ¬ 
ment.  Due  to  the  large  temperature  difference  between  the  droplet  and  ambient 
gases,  the  density  varies  significantly  from  dilute  gas  to  compressed  liquid,  and 
as  such  precludes  the  conventional  approaches  developed  for  incompressible  flows. 
On  the  other  hand,  most  existing  numerical  methods  developed  for  compressible 
flow  usually  exhibit  their  best  performance  at  moderate  and  high  Mach  numbers. 
When  the  Mach  number  becomes  very  small  as  in  the  present  case,  the  sensible 
part  of  the  fluid  energy  manifests  itself  in  the  form  of  pressure,  with  the  energy 
equation  decoupled  from  the  rest  of  the  governing  equations.  The  associated  dif¬ 
ference  among  the  eigenvalues  adversely  affects  the  convergence  characteristics. 
Second,  as  velocity  approaches  zero,  the  pressure  term  in  the  momentum  equa¬ 
tion  causes  large  round-off  errors  that  prevent  the  solution  from  reaching  sufficient 
accuracy.  This  error  arises  in  the  calculation  of  the  pressure  and  convective  mo¬ 
mentum  flux  terms  in  the  momentum  equation  when  the  computer  is  forced  to 
truncate  valuable  digits.  Since  these  two  terms  differ  greatly  in  magnitude  at 
low  Mach  number,  significant  errors  are  introduced  that  consequently  contaminate 
the  numerical  solution.  Chemical  reactions  exhibit  another  category  of  numerical 
difficulties  because  of  the  wide  ranges  of  time  and  length  scales  involved  in  the 
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calculations.  In  regions  where  reaction  rates  are  high,  species  concentrations  may 
vary  rapidly  over  a  much  shorter  period  of  time  than  those  of  the  fluid  dynamic 
processes.  This  phenomenon  hinders  the  numerical  convergence  efficiency,  and 
requires  a  prohibitively  large  number  of  iterations  to  obtain  converged  solution. 

The  above  difficulties  are  common  to  most  studies  of  droplet  vaporization 
and  combustion,  but  become  more  pronounced  at  elevated  pressure  due  to  the  high 
gasification  rate  of  the  droplet  and  rapid  variations  of  fluid  properties  with  tem¬ 
perature.  To  circumvent  these  problems,  a  dual  time-stepping  integration  method 
with  a  low  Mach-number  treatment  is  chosen  for  the  current  study[45,46].  The 
scheme  is  constructed  in  two  steps.  First,  since  the  roundoff  error  is  caused  by  the 
singular  pressure  gradient  term  in  the  momentum  equation  (the  pressure  term  is  of 
order  of  while  the  convective  term  is  of  order  unity  in  the  non-dimensional 

momentum  equations)  [46],  the  pressure  is  decomposed  as  the  sum  of  a  constant 
reference  pressure  and  a  gauge  pressure. 

p{x,y,t)  =  pQ^j>g{x,y,t)  (3.1) 

The  second  step  incorporates  pseudo-time  iteration  by  adding  a  well- conditioned 
artificial  time  derivative  term  to  the  conservation  laws.  The  equations  are  well 
posed  by  preconditioning  the  eigenvalues,  effectively  circumventing  the  singularity 
problem  at  low  Mach  numbers.  An  implicit  iterative  procedure  is  then  used  in 
pseudo  time.  This  measure  is  equivalent  to  preconditioning  the  time  derivative 
terms  of  the  species  concentration  equations  so  that  all  the  chemical  and  convec¬ 
tive  processes  proceed  at  approximately  the  same  numerical  rate.  As  a  result, 
a  converged  solution  in  pseudo  time  corresponds  to  a  time-accurate  solution  in 
physical  time. 
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3.1  One-Dimensional  Treatment 

In  formulating  the  numerical  integration  procedure,  Eq.  (2.9)  is  modified 
by  introducing  pseudo-time  terms,  composed  of  a  set  of  primitive  variables  (p^, 
Ur,  h,  and  Yi),  to  the  non-conservative  form  of  the  governing  equations.  In  order 
to  achieve  better  conservative  characteristics,  the  conservative  form  of  governing 
equations  is  restored  by  adding  the  modified  continuity  equation  to  the  momentum, 
energy,  and  species  equations.  The  resulting  conservation  equations  for  a  system 
with  iV  —  1  species  equations  in  1-D  spherical  coordinate  are  written  as 

/  ZdV  +  ^  [  QdV  +  f  {E-  E^)dA  -ug  f  QdA=  f  SdV  (3.2) 
dr  JaV  dtJAV  JA  Ja  Jav 

where 


(  Pg  y 

0 

0 

0  ... 

0  \ 

Uf 

Url^ 

P 

0 

0  ... 

0 

h 

and  r  = 

hl^-l 

pUr 

p 

0  ... 

0 

Yi 

Yil^ 

0 

0 

p  ... 

0 

K  Yn-1  ) 

Yn-i/(3 

0 

0 

0  ... 

p  / 

where  the  scaling  factor  ^  is  used  to  ensure  that  all  the  eigenvalues  of  the  system 
have  the  same  order  of  magnitude,  a  condition  normally  required  for  eflficient  nu¬ 
merical  convergence.  The  vectors  Q,  E,  Ey,  Eg,  and  S  have  the  same  definition 
as  those  in  Eq.  (2.9),  except  the  absolute  pressure  in  the  momentum  equation  is 
replaced  by  the  gauge  pressure  pg. 
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The  terms  at  the  (p+l)*^  time  level  need  to  be  linearized  for  the  construction 
of  an  implicit  time  marching  scheme.  The  linearized  representation  of  convective 
and  diifusion  fluxes  and  source  term  vectors  are  as  follows. 


QP+i  =  QP  +  TAZ 

EP+i  =  EP  +  AAZ 
dZ^^^ 

jjj+i  =  r'+Y 

or 

SP+i  =  5^’  +  DAZ  (3.3) 

where  T,  A,  R,  and  D  are  the  Jacobian  matrices  defined  as 


T 

A 

R 

D 


dZ 

dZ 

dE^ 


dZ 


Complete  expressions  of  these  matrices  can  be  found  in  Appendix  C. 

For  second-order  accuracy  in  physical  time,  Eq.  (3.2)  can  be  discretize  in  a 
fully  implicit  form  as 


rpP 

3TP 

.At 

2At 

)\aZ  +  477  f  \iA- UgT)AZ -R^^4^]dA  =  -{RHSy  (3.4) 
A  V  Jaa  L  or  i 


where 
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{RHSy  = 


3QP  -  4Q”  +  Q 
2At 


n—1 


-s  + 


^JJ 


E  —  UgQ)  —  Ev 


dA  +  AD  (3.5) 


AD  is  the  artificial  dissipation  which  is  a  blending  of  second-diiference  and  fourth- 
difference  dissipation  terms.  The  second- difference  terms  are  used  to  prevent  os¬ 
cillations  at  discontinuity  in  the  fiowfield,  while  the  fourth- difference  terms  are 
important  for  numerical  stability  and  convergence[47]. 

3.2  Two-Dimensional  Treatment 


To  facilitate  numerical  manipulation,  the  conservation  laws  of  physics,  Eqs.  (2.18)- 
(2.22),  can  be  transform  into  a  generalized  coordinate  system  in  the  following 
vector  form. 

dQ  d{E  -  E,)  d(P  -  F.)  , 

9?+  d(  +— - ^ 

where  the  quantities  ^  and  rj  are  the  spatial  coordinates  in  generalized  coordinates, 
and  vectors  Fy  and  S  are  defined  as 


0 


(3.7) 


E  =  + 


(3.8) 


K  =  j{izE,  +  (,F,) 


(3.9) 


f  =  j{VzE  +  7i,F), 


(3.10) 
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Fv  =  j{t)xEv  +TJyFv)  (3.11) 

S  =  \s  (3.12) 

u 


with  J  the  coordinate  transformation  Jacobian,  and  Vxi  Vy  metric 

terms.  The  conserved  variable  vector  Q  is  defined  as 


/  P  \ 


Q  = 


pu 

pv 

pet 


V  pY^  ) 

The  convective-flux  vectors,  E  and  jF,  in  the  axial  and  radial  directions  take  the 
forms 


pu  \ 

(  pv  \ 

to 

+ 

puv 

puv 

F  = 

pv^  d-  p 

{pet  +  p)u 

{pet  -f  p)v 

puYi  y 

\  pvYi  ) 

The  corresponding  diffusion-flux  vectors,  E^  and  Fv,  are 
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0  ^ 

0  \ 

^xx 

Trx 

Txr 

F,= 

Trr 

'^'^xx  "1“  ^'^XT  “1“  (9e)a; 

UTxt  +  VTrr  +  (?e)r 

\  {Qm,i)x  / 

\  iQTn,i)r  / 

The  source  vector  S  includes  parts  associated  with  the  axisymmetric  geometry  and 
remaining  viscous  terms  and  can  be  written  as 


/ 


5  = 


Iff  JL  n 


di) 


Pa  ~  +  Vx^)  -  |y(6|f  +  Vr^) 

dri  ) 


2  2 

dfiuv  I  _  duuv  I  ^  diiv  .  _  dav  \ 


V 


0 


} 


To  obtain  time-accurate  solutions  for  time-evolving  problems,  a  dual  time¬ 
stepping  integration  method  is  chosen.  Thus,  the  conservative  equations  in  gener¬ 
alized  curvilinear  coordinates  ,  tj)  are  modified  as 

dZ  BQ  d{E-E^)  d{F  -  F^) 


ar+ 


+ 


Bj] 


=  S 


(3.13) 


where  r  represents  the  pseudo  time.  The  primitive  variable  vector  Z  and  the 
preconditioning  matrix  F  are  expressed  as 
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(  P9  \ 


Z  =  r 


u 


h 


Yy 


\  Yn-,  j 


and  r  = 


(  1//?  0  0  0  0  ...  0  \ 

u! ji  p  0  0  0  ...  0 

vj  jd  0  p  0  0  ...  0 

hj ^  \  pu  pv  p  Q  . . .  Q 

Yil^  0  0  0  /o  ...  0 

[Yn-i/0  0  0  0  0  ...  pj 


where  the  vectors  Q,  E.,  E^,  F,  and  S  remain  the  same  definitions  as  those 
for  Eq.  (3.6).  The  quantities  J,  rjx,  //r,  ix,  and  are  derived  directly  from  the 
coordinate  transformation. 

The  equation  for  the  (p+1)^*  iteration  in  pseudo  time  and  the  (n+1)**  level 
in  physical  time  is  written  as 


r - 1 - 1 - —  + 


At 


At 


+ 


a(f- fir-  =5^^.  (3.14) 

07} 


Since  the  time  step  is  set  to  be  constant  and  the  resolution  accuracy  is  second- 
order,  ai,  02,  and  03  are  chosen  as 


a\  =  2’  «2  =  -2,  «3  =  2  (3-15) 

The  terms  at  the  (p  +  1)*^  time  level  in  Eq.  (3.14)  need  to  be  linearized  for  the 
construction  of  an  implicit  time  marching  scheme.  Using  a  first-order  Tayler  series 
expansion,  the  vectors  Q,  E,  Ey,  F,  F„,  and  S  may  be  linearized  with  respect  to 


Z  as  follows 
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gp+i  =  QP  +  TPAZ 
^P+i  ^  EP  +  APAZ 

El*'  =  b;  +  r^ja(|5) 

^P+1  ^  pP  +  BPAZ 

Fi*'  =  f;  +  ri„a{^) 

5P+1  =  SP  +  B^AZ  (3.16) 

The  eigenvalues  (in  ^  direction)  in  pseudo  time  can  be  obtained  from  the  matrix 

r~^A 

A.  =  :/, C/,  1  t/(l  +  |)  ±  ^(/2(l  -  +  4/?(af  +  al)  (3,17) 

where  c  is  the  speed  of  sound  with  ai  =  02  =  ^-nd  U  is  the  contravari- 

ant  velocity  component  defined  as  =  aiw  +  a2U.  To  obtain  well- conditioned 
eigenvalues,  the  scaling  factor  /?  is  taken  to  be 
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-  + AD  (3.20) 


where  A  and  B  are  the  convective  flux  Jacobians  in  ^  and  t)  directions,  respectively, 
D  the  source  term,  and  and  R,,,,  the  viscous  term  Jacobians.  Explicit  expres¬ 
sions  of  these  matrices  can  be  found  in  Appendix  C.  The  last  term  in  Eq.  (3.20) 
represents  the  fourth-order  dissipation  defined  as 


(3.21) 


where  cad  is  the  artificial  dissipation  rate. 

One  significant  advantage  of  the  dual  time-stepping  method  is  that  the 
convergence  of  the  iterative  process  is  determined  by  the  eigenvalue  characteristics 
in  the  pseudo-time  domain,  and  not  by  the  original  eigenvalues  which  become 
stiff  at  low  Mach  numbers.  This  feature  allows  flexibility  in  the  selection  of  time 
step  sizes  in  both  time  frames.  The  physical  time  step  is  determined  based  on 
the  characteristic  evolution  of  the  unsteady  flow  under  consideration,  while  the 
pseudo-time  step  depends  on  the  numerical  stability  of  the  algorithm  and  can  be 
adjusted  to  obtain  the  optimum  convergence. 


3.2.1  Modified  Strongly  Implicit  Solver 

A  new  iterative  method  was  proposed  by  Stone  et  al.  [48].  The  algorithm 
involving  a  5-point  configuration  is  briefly  listed  below  for  reference. 


+  aijAZi+ij  +  G?,AZ,,,  +  G'ijAZi,,-,  +  (3.22) 


where  Gij  is  an  nxn  block  matrix  and  qij  is  a  vector  with  length  n.  The  super¬ 
scripts  w,  e,  s,  n,  and  p  denote  respectively  the  cells  on  western,  eastern,  southern. 
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I 

I 


and  northern  sides  with  respect  to  the  p  cell.  The  above  equation  can  be  rewritten 
in  a  vector  form  as 

GAZ  =  q  (3.23) 

where 


/ 


Gn 


\ 


G  = 


GW  QP  ge 


\  G^  G^  GP  / 

This  matrix  can  be  approximated  into  two  triangle  matrice(L  and  U),  their  struc¬ 
ture  are 


V 
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\ 


U  = 


e 


where 


V  ••• ) 

^i,j  —  Glj 

=  Gfj 

Cij  =  G\j  —  —  bijdi-ij 

di,j  =  (cij)  ^G\j 

—  i^ij)  ^Gfj  (3-24) 


The  solution  of  AZ  can  be  obtained  by  forward  and  backward  substitutions  as 
follows. 


LAZ*  =  q  (3.25) 

UAZ  =  AZ*  (3.26) 


where  Z*  is  the  provisional  vector  and  is  calculated  from  Eq.  (3.25)  first  to  serve 
as  right  hand  side  of  Eq.  (3.26). 


CHAPTER  4 


DROPLET  COMBUSTION  IN  QUIESCENT 
ENVIRONMENTS 


The  combustion  of  an  isolated  n-pentane  fuel  droplet  in  air,  as  shown  in  Fig. 
4.1,  over  a  wide  range  of  pressure  has  been  studied  systematically.  The  droplet 
is  introduced  instantaneously  to  a  quiescent  air  environment  with  temperature  of 
1500  K.  The  initial  droplet  temperature  and  diameter  are  300  K  and  100  fxm,  re¬ 
spectively,  and  the  ambient  pressures  are  either  in  the  subcritical  or  supercritical 
regime  of  the  fuel  species.  The  critical  properties  of  the  major  constituent  species 
in  the  gas  and  liquid  phases  are  listed  in  Table  4.1.  Figure  4.2  shows  the  time  histo¬ 
ries  of  the  droplet  surface  temperature  at  various  pressures.  The  impulsive  rise  of 
the  surface  temperature  in  the  early  stage  of  the  droplet  lifetime  is  caused  by  large 
temperature  gradient  near  the  droplet  surface  (i.e.,  a  step  function  was  used  as 
the  initial  temperature  distribution).  Because  the  droplet  temperature  is  initially 
lower  than  that  of  the  ambient  gases,  heat  is  transferred  from  the  hot  air  toward 
the  droplet.  Part  of  the  heat  is  used  for  the  phase  transition  from  the  liquid  to 
the  vapor  phase,  and  the  remainder  goes  to  heating  up  the  droplet  interior.  Once 
ignition  is  achieved  in  the  gas  phase,  the  energy  feedback  from  the  flame  results  in 
a  rapid  increase  in  the  droplet  surface  temperature  which  can  be  observed  by  the 
change  of  slope  in  the  time  histories  of  surface  temperature  around  the  ignition 
time.  At  low  pressures  (p  <  20  atm),  the  surface  temperature  varies  slowly  follow¬ 
ing  onset  of  flame  development.  It  then  levels  olF  at  a  temperature  slightly  lower 
than  the  droplet  boiling  temperature,  a  condition  usually  referred  to  as  the  pseudo 
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Figure  4.1.  Schematic  of  a  Burning  Fuel  Droplet. 

wet-bulb  state[49].  Under  this  condition,  nearly  all  of  the  energy  transferred  from 
the  ambient  gas  is  used  for  the  vaporization,  with  only  a  small  fraction  is  used 
for  droplet  heatup.  As  the  pressure  increases  to  33.3  atm,  the  transport  processes 
between  gas  and  liquid  are  commensurate  with  each  other,  the  entire  gasification 
process  becomes  transient  with  a  continuous  increase  of  liquid  temperature,  result¬ 
ing  in  the  disappearance  of  the  pseudo  wet-bulh  state.  At  p  =  80  atm,  the  droplet 
surface  even  reaches  the  critical  mixing  point  (at  t  =  3.61  ms)  before  the  end  of 
its  lifetime.  When  this  occurs,  the  fluid  density  and  its  gradient  vary  smoothly 
throughout  the  domain,  the  entire  flowfield  becomes  a  continuous  medium  with  no 
distinct  liquid/vapor  interface  as  can  be  observed  for  the  subcritical  cases.  Since 
there  is  no  phase-equilibrium  criterion  posed  on  the  droplet  surface,  the  surface 
which  attains  the  critical  mixing  temperature  usually  regresses  at  a  speed  differ¬ 
ent  from  that  with  the  critical  mixing  composition  due  to  the  non-unity  of  the 


Figure  4.2.  Time  Variations  of  Droplet  Surface  Temperature  at  Various  Pressures, 
n-Pentane/Air  System. 

Lewis  number.  For  the  convenience  of  discussion,  the  droplet  gasification  rate  is 
defined  herein  in  accordance  with  the  motion  of  the  surface  which  attains  the  crit¬ 
ical  mixing  temperature.  Within  the  present  treatment,  since  only  the  dissolution 
of  nitrogen  in  the  liquid  fuel  is  taken  into  account  in  the  phase-equilibrium  anal¬ 
ysis,  the  critical  mixing  temperature  depends  solely  on  the  pressures  for  such  a 
binary  system,  as  indicated  by  the  horizontal  straight  lines  in  the  plot  for  p  =  80, 
110,  and  140  atm.  The  critical  properties  and  milestone  times  associated  with  the 
fuel-droplet  combustion  are  given  in  Table  4.2.  As  the  ambient  pressure  increases, 
the  high  concentrations  of  oxygen  in  the  gas  phase  and  the  fuel  vapor  issued  from 
the  droplet  surface  result  in  a  high  chemical  reaction  rate,  consequently  causing  a 
progressive  decrease  in  ignition  time  (defined  as  the  time  when  the  peak  gas-phase 
temperature  exceeds  its  initial  value  by  100  K,  see  Table  4.2).  In  addition,  the 
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TABLE  4.1 

Critical  Properties  of  Major  Constituent  Species  in  the  n-Pentane/Air  System. 


n-Pentane 

Nitrogen 

Oxygen 

Carbon  Dioxide 

Water 

T,  (K) 

469.6 

126.2 

154.5 

304.2 

647.3 

Pc  {atm) 

33.3 

33.5 

49.8 

72.8 

217.6 

Vc  {cm^fg  —  mole) 

304 

89.5 

73.4 

94 

56 

Zc 

0.262 

0.290 

0.288 

0.274 

0.229 

surface  temperature  jump  during  the  flame-development  stage  appears  an  increase 
with  increasing  pressure.  Since  the  critical  mixing  temperature  decreases  with 
pressure,  the  droplet  reaches  its  critical  condition  more  easily  at  higher  pressures, 
almost  immediately  following  establishment  of  the  diffusion  flame  in  the  gas  phase 
for  p  >  110  atm.  The  results  displayed  here  are  not  very  sensitive  to  the  ignition 
transient  for  ambient  pressure  greater  than  20  atm.  If  the  preexponential  constant 
in  the  Arrhenius  law,  Eq.  (2.114)  artificially  increases  by  a  factor  of  10,  the  burnout 
times  changes  by  less  than  5  %  for  p  >  20  atm.  The  lack  of  sensitivity  in  droplet 
lifetime  calculations  to  the  ignition  process  and  chemical  kinetics  at  high  pressure 
may  be  attributed  to  the  fact  that  the  ignition  transient  period  is  taken  out  only 
a  very  small  fraction  of  the  total  burnout  time. 

Figure  4.3  shows  the  temporal  variations  of  the  latent  heat  of  vaporization 
at  various  pressures.  Since  the  latent  heat  of  vaporization  depends  on  both  tem¬ 
perature  and  pressure,  it  decreases  with  increasing  temperature  as  well  as  pressure, 
and  reduces  acceleratively  to  zero  as  the  critical  mixing  point  is  approached[49]. 
The  sudden  decrease  in  each  of  the  latent-heat  profiles  after  ignition  arises  from 
the  rapid  increase  in  the  droplet  surface  temperature  during  the  flame  develop¬ 
ment  stage.  At  p  =  80  atm,  another  rapid  decrease  in  the  latent  heat  is  observed 


TABLE  4.2 

Milestone  Times  and  Critical  Properties  for  Combustion  of  n-Pentane  Fuel  Droplet  in  Air  at  Various  Pressures; 
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Figure  4.3.  Time  Variations  of  Latent  Heat  of  Vaporization  at  Various  Pressures, 
n-Pentane/Air  System. 

when  the  droplet  surface  approaches  its  critical  mixing  state.  The  situation  be¬ 
comes  more  conspicuous  at  higher  pressures  (e.g.,  p  =  110,  140  atm),  because  the 
ignition  takes  place  slightly  prior  to  occurrence  of  the  critical  mixing  point  at  the 
droplet  surface  (see  Table  4.2  for  numerical  values).  The  overlapping  effect  causes 
the  latent  heat  to  drop  almost  instantaneously  to  zero  after  onset  of  the  flame. 

Figure  4.4  depicts  the  time  variations  of  the  gaseous  n-pentane  mole  fraction 
at  the  droplet  surface.  In  the  initial  stage  during  which  the  droplet  surface  tem¬ 
perature  is  relatively  low,  the  n-pentane  mole  fraction  decreases  with  increasing 
pressure.  Figure  4.4(a)  shows  that  as  the  droplet  heats  up,  the  rapid  increase  in  the 
surface  temperature  at  high  pressures  (e.g.,  p  =  50  atm)  allows  the  fuel- vapor  mole 
fraction  to  grow  at  a  faster  rate  and  consequently  exceeds  that  at  lower  pressures 
(e.g.,  p  =  5  atm).  When  the  droplet  surface  approaches  the  critical  mixing  state 
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for  >  80  atm,  the  n-pentane  concentration  increases  exceedingly  fast  and  finally 
stabilizes  at  a  level  set  by  the  critical  mixing  condition,  as  shown  in  Fig.  4.4(b). 

The  time  histories  of  liquid  n-pentane  fuel  and  dissolved  nitrogen  at  the 
droplet  surface  are  shown  in  Fig.  4.5.  In  the  subcritical  regime,  with  the  ambient 
pressure  less  than  the  critical  pressure  of  pure  n-pentane,  the  amount  of  nitrogen 
dissolved  in  the  liquid  phase  is  very  limited.  It  decreases  with  increasing  temper¬ 
ature  and  finally  vanishes  at  the  boiling  temperature  of  the  fuel.  However,  the 
situation  in  a  supercritical  environment  exhibits  the  opposite  trend.  The  ambient 
gas  solubility  in  the  droplet  is  quite  significant,  and  increases  with  temperature. 
For  instance,  the  liquid  phase  contains  more  than  35  %  nitrogen  at  the  surface  for 
p  =  110  atm  when  the  critical  mixing  condition  is  reached.  The  ensuing  change  in 
the  thermophysical  properties  of  the  droplet  may  cause  an  intrinsic  modification 
of  the  fuel  gasification  behavior. 

Figure  4.6  presents  the  histories  of  the  square  of  droplet  diameter  at  various 
pressures.  The  droplet  surface  regression  rate  increases  persistently  with  pressure, 
mainly  due  to  the  reduced  latent  heat  of  vaporization  at  high  pressures.  Although  a 
decrease  of  mass  diffusivity  takes  place  at  high  pressures  and  gives  rise  to  an  adverse 
influence  on  the  gasification  process,  this  effect  is  usually  overridden  by  the  reduced 
latent  heat  of  vaporization.  Consequently,  a  net  increase  of  the  surface  regression 
rate  is  observed.  At  high  pressures,  the  relatively  fast  increase  in  temperature 
leads  to  a  rapid  decrease  in  liquid  density  during  the  initial  heatup  period.  The 
resulting  volume  expansion  may  offset  the  regression  arising  from  the  vaporization 
process,  and  cause  the  droplet  to  exceed  its  initial  droplet  diameter  which  varies 
slowly  during  the  initial  heatup  period,  with  very  little  gas-phase  reaction  involved 
in  the  gasification  mechanism.  It  then  regresses  at  a  much  faster  rate  subsequent 


Figure  4.4.  Time  Variations  of  n-Pentane  Mole  Fraction  on  Droplet  Surface  at 
Various  Pressures,  n-Pentane/Air  System:  (a)  p  <  50  atm.  (b)  p  >  50  atm. 
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Figure  4.5.  Time  Variations  of  Liquid  Mole  Fraction  on  Droplet  Surface  at  Various 
Pressures,  n-Pentane/Air  System:  (a)  n-Pentane.  (b)  Nitrogen. 
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Figure  4.6.  Time  Variations  of  the  Square  of  Droplet  Diameter  at  Various  Pres¬ 
sures,  n- Pentane/ Air  System. 

to  onset  of  the  flame  in  the  gas  phase,  mainly  because  of  the  enhanced  heat  transfer 
to  the  liquid  phase  and  the  increased  mass  transfer  arising  from  the  higher  fuel- 
vapor  concentration  at  the  surface.  At  low  pressures  (p  <  20  atm),  when  the 
droplet  reaches  its  pseudo  wet-bulb  state,  the  surface  area  varies  almost  linearly 
with  time,  and  the  classical  d^-law  provides  a  faithful  description  of  the  regression 
behavior.  However,  the  situation  at  high  pressures  becomes  quite  different.  The 
continuous  variations  of  the  droplet  interior  and  surface  properties  prohibit  the 
use  of  the  d^-law  in  predicting  the  surface  regression  rate.  Furthermore,  as  the 
droplet  surface  reaches  its  critical  mixing  point,  the  latent  heat  of  vaporization 
drops  rapidly  to  zero  and  the  regression  rate  accelerates,  as  indicated  by  the  kinks 
in  the  proflles  for  p  >  80  atm.  The  rapid  increase  of  the  regression  rate  caused 
by  the  change  of  gasification  mechanism  after  reaching  the  critical  mixing  point  is 
a  predominant  phenomenon  in  high-pressure  vaporization  processes  and  must  be 
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Figure  4.7.  Time  Variations  of  Fuel  Droplet  Mass  at  Various  Pressures,  n- 
Pentane/  Air  System. 


considered  in  a  supercritical  spray  analysis.  When  the  ambient  pressure  increases 
from  110  to  140  atm,  the  reduced  mass  diffusion  results  in  a  slowdown  in  the 
increase  of  regression  rate  with  pressure. 

The  time  variations  of  the  liquid  fuel  mass  are  displayed  in  Fig.  4.7.  Although 
Fig.  4.6  shows  the  time  history  of  the  surface  regression  rate,  little  information 
about  the  fuel  mass  gasification  rate  is  given  unless  the  variable-property  effect  is 
properly  taken  into  account.  This  point  can  be  easily  elaborated  by  considering 
the  global  mass  balance  of  the  liquid  fuel. 
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where  mp  is  the  total  mass  of  the  liquid  fuel.  The  gasification  rate  contains  contri¬ 
butions  from  both  surface  regression  and  fuel-density  variation,  as  indicated  by  the 
first  and  second  terms  on  the  right-hand  side  of  the  above  equation,  respectively. 

Figures  4.8(a)  and  (b)  show  the  spatial  distributions  of  temperature  and 
density  in  the  droplet  interior,  respectively.  The  droplet  surface  temperature  in¬ 
creases  rapidly  in  the  initial  stage,  then  remains  fixed  after  it  reaches  the  critical 
mixing  state.  As  a  result,  all  of  the  thermal  energy  transferred  from  the  ambient 
gas  is  used  to  heat  up  the  droplet,  which,  in  turn,  causes  the  critical  surface  to 
propagate  toward  the  interior  of  the  liquid  phase.  The  temperature  also  indicates 
a  continuous  reduction  of  liquid  density,  as  shown  in  Fig.  4.8(b).  In  particular,  a 
steep  variation  of  density  occurs  near  the  droplet  critical  surface,  since  the  liquid 
density  decreases  exceedingly  fast  with  temperature  in  the  vicinity  of  the  critical 
mixing  point.  Two  points  must  be  noted  during  this  process.  First,  the  time- 
varying  reduction  of  liquid  density  causes  the  droplet  to  swell;  the  ensuing  volume 
dilatation  is  significant  and  must  be  accounted  for  in  determining  the  droplet  re¬ 
gression  behavior.  Second,  the  non-uniformity  of  density  distribution  results  in  a 
significant  outward  convection  that  may  effectively  retard  energy  transfer  to  the 
droplet  interior. 

The  time  histories  of  flame  development  can  be  examined  by  plotting  the 
distributions  of  gas-phase  temperature  and  species  concentrations  at  various  times. 
Figures  4.9(a)  and  (b)  show  a  typical  set  of  temperature  profiles  for  p  =  80  atm. 
The  corresponding  distributions  of  n-pentane  species  concentration  are  given  in 
Figs.  4.10(a)  and  (b).  The  droplet  ignition  is  achieved  around  t  =  1.28  ms.  The 


DENSITY,  kg/m^  TEMPERATURE 


550 

^  500 

450 

400 

350 

300 

ao  0.2  0.4  0.6  0.8  1.0 

REDUCED  RADIAL  DISTANCE.  r/R(t) 

800 

700 

600 

500 

400 

300 

200 

0.0  0.2  0.4  0.6  0.8  1.0 

REDUCED  RADIAL  DISTANCE.  r/R(t) 


Figure  4.8.  Distributions  of  Liquid-Phase  Properties  at  Various  Times,  n- 
Pentane/  Air  System  at  p  =  80  atm:  (a)  Temperature,  (b)  Density. 
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Figure  4.9.  Distributions  of  Gas-Phase  Temperature  at  Various  Times,  n- 
Pentane/  Air  System  at  p  =  80  atm:  (a)  prior  to  End  of  Droplet  Gasification 
Lifetime,  (b)  after  End  of  Droplet  Gasification  Lifetime. 
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Figure  4.10.  Distributions  of  Gas-Phase  Species  Concentrations  at  Various  Times, 
n-Pentane/Air  System  at  p  =  80  atm:  (a)  prior  to  End  of  Droplet  Gasification 
Lifetime,  (b)  after  End  of  Droplet  Gasification  Lifetime. 
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resulting  flame  then  spreads  outward  progressively  during  the  droplet  gasification 
lifetime  {t  <  7.94  ms).  After  completion  of  the  fuel  gasification  process,  the  residual 
fuel  vapor  in  the  gas  phase  continues  to  react  with  air  and  causes  the  flame  to 
shrink.  The  gas-phase  reaction  ends  at  t  =  11.90  ms,  at  which  point  the  fuel  vapor 
is  completely  consumed.  The  overall  flame  development  is  further  illustrated  by 
the  photographs  in  Fig.  4.11.  The  top  figure  shows  the  time  history  of  the  droplet 
surface  temperature.  The  instantaneous  temperature  contours,  including  both  the 
droplet  interior  and  ambient  gases,  are  presented  at  six  different  times.  The  color 
scale  at  the  left  indicates  temperature,  increasing  from  the  blue  (T  =  300  K) 
toward  the  red  {T  —  2900  K).  The  dark  blue  region  at  the  center  corresponds  to 
the  liquid  droplet.  The  motion  of  the  flame  front  is  clearly  observed,  as  shown  by 
the  red  zone. 

The  flame  development  process  is  further  studied  by  examining  the  loci  of 
the  peak  gas-phase  temperatures  at  various  ambient  pressures.  Figures  4.12(a)  and 
(b)  show  the  flame  movement  (as  depicted  by  the  location  of  the  peak  temperature) 
in  terms  of  the  reduced  distance.  Owing  to  the  relatively  higher  concentrations  of 
fuel  vapor  and  oxygen  and  slower  diffusion  of  fuel  vapor  away  from  the  surface  at 
higher  pressures,  the  gas-phase  ignition  occurs  at  a  radial  position  which  contin¬ 
ues  to  move  closer  to  the  droplet  surface  with  increasing  pressure.  Immediately 
after  ignition,  the  flame  front  moves  inward  due  to  the  shift  in  the  location  of  the 
stoichiometric  composition.  It  then  spreads  outward  because  of  the  increase  in 
the  gasification  rate,  until  the  liquid  droplet  disappears.  Finally,  the  termination 
of  the  fuel  supply  from  the  liquid  droplet  causes  the  flame  to  shrink  and  eventu¬ 
ally  disappear  at  the  burnout  point.  During  the  flame  history,  the  peak  gas-phase 
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temperature  increases  continuously,  except  for  the  slight  set  back  following  the  on¬ 
set  of  the  flame.  The  peak  temperature  increases  with  pressure  in  the  subcritical 
region (p  <  50  atm),  due  to  the  increased  chemical  reaction  rates  at  higher  pres¬ 
sures.  On  the  other  hand,  for  p  >  50  atm,  the  peak  flame  temperature  decreases 
as  pressure  increases.  This  is  mainly  caused  by  the  decrease  in  fuel-vapor  mass 
diflFusivity  at  high  pressures,  which  subsequently  results  in  an  increase  in  the  Lewis 
number  (note  that  pressure  exerts  a  much  stronger  adverse  influence  on  mass  dif¬ 
fusion  than  on  thermal  diffusion).  As  a  consequence  of  the  reduced  rate  of  mass 
diffusion  of  the  gaseous  fuel  to  the  flame  zone,  the  rate  of  thermal  energy  produced 
by  chemical  reactions  also  decreases,  rendering  a  lower  flame  temperature  as  the 
rate  of  heat  transferred  away  from  the  flame  varies  only  slightly  with  pressure. 

Figure  4.13(a)  presents  the  effect  of  pressure  on  various  milestone  times 
associated  with  droplet  gasification  and  burning  process.  Here  gasification  life¬ 
time  is  the  time  required  for  complete  gasification,  droplet  burning  lifetime  is  the 
gasification  lifetime  minus  ignition  time,  single-phase  combustion  lifetime  is  the 
time  duration  from  complete  gasification  to  burnout  of  all  fuel  vapor,  and  combus¬ 
tion  lifetime  is  the  sum  of  single-phase  combustion  lifetime  and  droplet  burning 
lifetime.  The  gasification  lifetime  decreases  continuously  with  pressure,  whereas 
the  single-phase  combustion  lifetime  increases  progressively  with  pressure  due  to 
its  adverse  effect  on  mass  diffusion.  More  importantly,  the  pressure  dependence 
of  combustion  lifetime  exhibits  irregular  behavior.  This  phenomenon  may  be  at¬ 
tributed  to  the  overlapping  effect  of  reduced  enthalpy  of  vaporization  and  mass 
diffusion  with  increasing  pressure. 

Since  the  time  scales  for  diffusion  processes  are  inversely  proportional  to 
the  droplet  diameter  squared,  it  is  important  to  examine  the  effect  of  droplet 
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size  on  the  burning  characteristics.  In  this  regard,  calculations  were  conducted 
for  large  droplets  having  an  initial  diameter  of  1000  /zm,  which  is  comparable  to 
the  sizes  considered  in  most  experimental  studies  of  supercritical  droplet  combus- 
tion[ll,50].  Furthermore,  the  ignition  transient  occupies  only  a  very  small  fraction 
of  the  entire  droplet  lifetime  for  large  droplets.  The  uncertainties  associated  with 
the  ignition  procedure  in  determining  the  characteristics  of  droplet  gasification 
can  be  minimized.  Consequently,  a  more  meaningful  comparison  with  experimen¬ 
tal  data  can  be  made.  Figure  4.13(b)  shows  the  milestone  times  for  Dq  =  1000  fim. 
The  combustion  behavior  of  a  large  droplet  reveals  several  characteristics  distinct 
from  those  of  a  small  droplet.  First,  ignition  for  large  droplets  occurs  in  the  very 
early  stage  of  the  entire  droplet  lifetime.  The  influence  of  gasification  prior  to 
ignition  on  the  overall  burning  mechanisms  appears  to  be  quite  limited.  Second, 
the  combustion  lifetime  decreases  with  increasing  pressure,  reaching  a  minimum 
near  the  critical  pressure  of  the  liquid  fuel.  As  the  pressure  further  increases,  the 
combustion  time  increases  due  to  reduced  mass  diffusivity  at  high  pressures.  The 
gasification  lifetime  decreases  continuously  with  pressure,  whereas  the  single-phase 
combustion  lifetime  increases  progressively  with  pressure.  This  phenomenon  is  in 
qualitative  agreement  with  the  experimental  observations  of  Faeth  et  a/.[ll]  and 
Sato  et  a/.  [50],  and  suggests  a  change  of  the  burning  characteristics  across  the  crit¬ 
ical  pressure.  At  low  pressures,  the  gasification  of  liquid  fuel  primarily  controls  the 
combustion  process,  while  in  a  supercritical  environment,  the  transient  gas-phase 
diffusion  plays  a  more  important  role. 

While  the  present  work  provides  new  insights  to  and  enhances  current  un¬ 
derstanding  of  supercritical  combustion  of  fuel  droplets,  a  quantitative  comparison 
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Figure  4.13.  Effect  of  Pressure  on  Milestone  Times  Associated  with  Droplet 
Gasification  and  Burning  Processes,  n-Pentane/Air  System:  (a)  Dq  =  100  ^m. 
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of  calculated  results  with  experimental  data  is  rather  difficult  to  make,  mainly  be¬ 
cause  of  the  neglect  of  radiation  heat  transfer  in  the  analysis.  The  lack  of  reliable 
radiant  transport  properties  at  high  pressures  also  poses  a  problem.  Thermal  ra¬ 
diation  from  the  flame  to  the  droplet  can  be  substantial  for  large  droplets  at  high 
pressures.  Thus,  without  a  full  account  of  radiation  heat  transfer  in  the  analysis, 
any  attempt  of  quantitative  comparison  will  be  futile. 


CHAPTER  5 


PRESSURE-COUPLED  VAPORIZATION  RESPONSES 


Combustion  instability  has  haunted  development  of  many  high  performance 
combustors  for  almost  forty  years.  This  problem  severely  impairs  engine  opera¬ 
tion  and  often  leads  to  catastrophic  consequence[2].  For  the  past  three  decades,  a 
considerable  amount  of  elfort  has  been  devoted  to  seek  a  deeper  understanding  of 
this  undesirable  phenomenon.  It  is  generally  accepted  that  unsteady  motions  will 
grow  if  heat  and  mass  are  added  in  phase  with  unsteady  waves.  Conversely,  it  will 
damp  if  the  mass  and  energy  additions  are  out  of  phase  with  pressure  oscillation. 
For  most  liquid-fueled  systems,  propellants  are  delivered  into  a  combustion  cham¬ 
ber  as  a  spray  of  droplets.  Since  droplet  vaporization  represents  a  rate  controlling 
process  in  the  combustors,  the  dynamic  behavior  of  spray  combustion  is  indeed 
a  statistic  consequence  of  the  vaporization  characteristic  of  an  individual  droplet. 
A  knowledge  of  the  vaporization  response  of  an  isolated  droplet  to  ambient  flow 
oscillations  becomes  essential  in  treating  combustion  instability [51]. 

Investigation  of  dynamic  responses  of  droplet  vaporization  and  combustion 
to  ambient  flow  fluctuations  started  in  the  60s.  Strahle[52]  examined  droplet¬ 
burning  responses  in  a  forced  convective  field  in  which  a  small-amplitude  sound 
wave  was  introduced  into  the  free  stream.  Since  the  exact  solution  can  only  be 
found  near  the  forward  stagnation  point,  his  primary  interests  were  focused  on 
the  local  vaporization  responses  under  the  action  of  forced  oscillations.  He  exam¬ 
ined  the  boundary  layer  surrounding  the  droplet,  and  showed  that  the  amplitude 
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of  vaporization  response  increases  with  increasing  frequency.  However,  his  model 
erroneously  predicts  that  the  vaporization  response  to  pressure  oscillation  is  un¬ 
bounded,  approaching  infinity  as  the  square  root  of  frequency  goes  to  infinity.  Hei- 
dmann  and  Wieber[53]  studied  the  vaporization  processes  of  an  n-heptane  droplet 
in  combustors  with  superimposed  traveling  transverse  oscillations.  The  response 
functions  were  evaluated  over  fairly  wide  range  of  flow  conditions,  and  the  data 
was  correlated  with  a  dimensionless  parameter  which  is  the  frequency  normalized 
by  one  half  of  the  droplet  lifetime.  In  their  study,  the  peak  response  occurs  when 
the  droplet  lifetime  approximately  equal  to  the  oscillation  period.  Negative  gains 
(phase  lag  is  more  than  90°)  may  happen  when  the  vaporization  time  is  greater 
than  about  three  times  the  oscillation  period.  Allison  and  Fa.eth[54]  studied  the 
response  function  of  a  burning  liquid  monopropellant  surface  to  imposed  pressure 
oscillation,  and  found  that  the  combustion  response  tends  to  peak  in  two  frequency 
ranges;  one  corresponding  to  interaction  with  liquid-phase  transient  effects,  which 
could  be  verified  both  theoretically  and  experimentally;  the  other  corresponding  to 
gas-phase  transient  effects,  which  could  only  be  examined  theoretically.  Recently, 
Tong  and  Sirignano[55-57]  implemented  transient  heating  and  mass  diffusion  mod¬ 
els  for  the  liquid  phase  and  proposed  that  the  unsteady  droplet  vaporization  is  a 
potential  mechanism  for  driving  combustion  instabilities.  Although  the  preceding 
studies  have  provided  significant  information  about  the  underlying  mechanisms  in¬ 
volved  in  droplet  vaporization  response,  a  number  of  fundamental  problems  remain 
unresolved  for  high-pressure  conditions. 

In  the  following,  a  series  of  numerical  simulations  of  droplet  vaporization 
in  a  forced-oscillatory  environment  are  performed.  Because  the  model  is  based  on 
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Figure  5.1.  Schematic  of  An  Isolated  Droplet  Vaporizing  in  Forced  Oscillatory 
Environment. 

first  principles,  a  systematic  investigation  into  couplings  between  pressure  oscilla¬ 
tion  and  transient  mechanisms  involved  in  droplet  vaporization  becomes  feasible. 
Figure  5.1  shows  the  situation  examined,  a  vaporizing  droplet  in  a  quiescent  envi¬ 
ronment  superimposed  with  a  periodic  pressure  oscillation.  This  approach  allows 
researchers  to  focus  on  the  various  mechanisms  that  influence  the  dynamic  response 
of  droplet  vaporization  to  ambient  pressure  oscillations.. 

5.1  Discussion  of  Results 

Vaporization  responses  of  an  isolated  droplet  to  ambient  pressure  oscilla¬ 
tions  at  high  pressures  have  been  conducted  systematically.  The  analysis  treats  an 
n-pentane  droplet  in  nitrogen  gas.  The  initial  droplet  temperature  is  300  K  and 
the  ambient  nitrogen  temperature  is  1500  K.  The  mean  pressure  is  either  in  the 
subcritical  (10  and  20  atm)  or  supercritical  (50  and  80  atm)  regime  of  the  fuel.  In 
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the  early  stage  of  the  droplet  lifetime,  a  large  temperature  gradient  exists  across 
the  liquid/ vapor  interface,  rendering  a  strong  heat  flux  conducted  into  the  droplet 
interior.  Part  of  the  energy  is  used  to  gasify  the  liquid,  while  the  remainder  goes 
to  heat  up  the  droplet  interior.  As  a  consequence,  droplet  surface  temperature 
increases  and  may  lead  to  two  different  scenarios  (i.e.,  subcritical  and  supercritical 
vaporization).  First,  if  the  ambient  pressure  is  much  less  than  the  critical  value 
of  the  fuel,  the  surface  temperature  increases  and  approaches  the  pseudo  wet-bulb 
state[49].  The  droplet  surface  provides  a  well-defined  boundary  which  separates 
the  liquid  from  the  ambient  gases.  The  vaporization  rate  can  be  easily  quantified 
by  evaluating  the  mass  flux  through  the  surface.  Second,  if  the  ambient  tempera¬ 
ture  and  pressure  are  both  in  the  supercritical  regime  of  the  fuel,  the  characteristic 
times  of  the  transport  processes  in  both  liquid  and  gas  phases  approach  the  same 
order  of  magnitude  as  a  result  of  increasing  surface  temperature.  At  some  point 
during  the  droplet  lifetime,  the  distinction  between  these  two  phases  disappears, 
rendering  a  continuous  medium  in  the  entire  domain.  Under  this  condition,  the 
transport  processes  in  the  liquid  phase  become  as  important  as  those  in  the  gas 
phase,  the  vaporization  response  is  expected  to  exhibit  phenomena  different  from 
the  low  pressure  cases.  As  mentioned  in  Chapter  4,  the  droplet  regression  is  best 
characterized  by  the  motion  of  the  surface  which  attains  the  critical  mixing  temper¬ 
ature.  The  vaporization  rate  can  be  obtained  by  tracking  the  time  variation  of  the 
droplet  residual  mass  which  is  the  total  mass  confined  by  the  critical-temperature 
surface.  It  is  worth  noting  that  the  critical  mixing  properties  vary  with  pressure 
fluctuation.  To  avoid  unnecessary  confusion,  the  surface  for  supercritical  gasifica¬ 
tion  is  selected  based  on  the  critical  value  at  the  mean  pressure. 


91 


Figure  5.2.  Temporal  Evolution  of  Droplet  Surface  Temperature  at  Various  Mean 
Pressures,  n-Pentane/Nitrogen  System;  Dq  =  100  //m,  Too  =  1500  K,  To  =  300  K. 

Figure  5.2  shows  the  temporal  evolution  of  droplet  surface  temperature  at 
various  mean  pressures.  The  rapid  heatup  immediately  after  t  =  0  results  from 
the  large  temperature  gradient  near  the  surface  (i.e.,  a  step  function  was  used 
as  the  initial  temperature  distribution).  At  low  to  moderate  pressures  (p  <  20 
atm),  the  droplet  surface  temperature  increases  rapidly  during  the  heat  up  period, 
then  levels  olf  slight  lower  than  the  boiling  temperature  of  n-pentane.  As  the 
pressure  increases  to  50  atm,  the  surface  temperature  increases  continuously,  and 
no  pseudo  wet-bulb  state  can  be  found  during  the  droplet  lifetime.  For  p  =  80  atm, 
the  droplet  surface  even  reaches  its  critical  mixing  state  at  f  =  8.3  ms.  When 
this  occurs,  the  droplet  surface  essentially  becomes  a  critical  surface  where  the 
liquid  and  gas  phases  coexist  at  the  critical  mixing  condition.  Figure  5.3  displays 
the  time  histories  of  droplet  diameter  squared  at  various  pressures.  The  droplet 
surface  regression  rate  increases  persistently  with  pressure. 
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Figure  5.3.  Time  Histories  of  the  Droplet  Diameter  Square  at  Various  Pressures, 
n-Pentane/Nitrogen  System;  Dq  =  100  //m,  T<x>  =  1500  K,  To  =  300  K. 


To  quantify  the  dynamic  behavior  of  droplet  vaporization,  a  pressure- coupled 
response  function  is  defined  as  follows. 


Rri  = 


(w  —  w)/w 


(5.1) 


^  \ip-p)/p\ 

where  w  is  the  instantaneous  vaporization  rate  in  an  oscillating  condition,  and 
w  is  that  in  a  quiescent  environment.  Figure  5.4  presents  the  time  variations 
of  droplet  vaporization  rate  and  pressure  oscillations  for  p  =  50  atm,  marked 
respectively  by  the  solid  and  dashed  lines.  The  vaporization  rate  does  not  respond 
instantaneously  to  the  pressure  variation,  showing  a  time  lag  between  these  two 
periodic  oscillations.  The  amplitude  of  the  induced  vaporization  fluctuation  varies 
throughout  the  droplet  lifetime. 

To  clarify  the  physical  mechanisms  involved  in  the  transient  response  of 
droplet  vaporization  to  ambient  pressure  oscillation,  it  is  useful  to  examine  the 
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Figure  5.4,  Time  Variations  of  Pressure  and  Droplet  Vaporization  Rate  Fluctua¬ 
tions. 

variations  of  the  various  time  scales  throughout  the  droplet  lifetime.  Five  char¬ 
acteristic  times  are  generally  of  interest;  droplet  residual  time  ri,  droplet  lifetime 
r2,  liquid  thermal  inertia  time  T3,  liquid  thermal  diffusion  time  r4,  and  gas- phase 
thermal  diffusion  time  rs,  as  listed  in  Table  5.1.  The  droplet  lifetime  T2  should 
not  be  confused  with  the  droplet  residual  time  ri  which  is  a  measure  of  the  re¬ 
maining  lifetime  of  a  droplet  at  a  given  instant  of  time.  Figure  5.5  shows  the  time 
variations  of  these  characteristic  times  for  poo  =  10  atm.  The  droplet  residual 
time  Ti  decreases  persistently  with  time.  The  other  characteristic  times  {t2,  T3, 
T4,  and  Ts)  first  increase  in  the  early  stage  and  then  drop  to  zero  at  the  end  of 
the  droplet  lifetime.  During  the  heatup  period,  the  predominant  factor  is  the  in¬ 
crease  of  surface  temperature  induced  by  the  energy  transferred  from  the  ambient 


TABLE  5.1 

Characteristic  Times  Involved  in  Droplet  Vaporization. 


94 


Tl 

Droplet  residual  time 

- 1 

T2 

Droplet  lifetime 

T3 

Liquid  thermal  inertia 

A. 

T4 

Liquid  thermal  diffusion 

PeCp/rl 

Ts 

Gets  phase  diffusion 

psCp^s'^g 

__A _ 

gases.  The  increasing  temperature  modifies  thermophysical  properties  (i.e.,  den¬ 
sity,  specific  heat,  thermal  conductivity)  on  the  droplet  surface,  showing  increases 
of  characteristic  times.  When  the  surface  condition  approaches  the  thermodynamic 
wet-bulb  stage,  the  decreases  in  characteristic  times  are  mainly  contributed  by  the 
decreasing  droplet  radius. 

Figure  5.6  represents  the  effect  of  frequency  on  the  vaporization  responses. 
The  initial  droplet  temperature  and  diameter  are  300  K  and  100  //m,  respectively, 
and  the  ambient  mean  pressure  is  50  atm.  Three  different  frequencies(i.e.,  500, 
2000,  and  6000  Hz)  are  considered.  Since  the  physical  mechanisms  associated  with 
the  vaporization  process  are  intrinsically  transient  in  time,  the  pressure  oscilla¬ 
tion  can  easily  couple  with  these  processes  and  exert  significant  influences  on  the 
vaporization  rate.  Results  indicate  that  the  magnitude  of  vaporization  response 
increases  with  frequency.  In  the  high  frequency  regime,  the  time  period  for  pres¬ 
sure  oscillation  is  much  smaller  than  the  liquid  thermal  inertia  time  in  the  droplet 
interior.  The  variation  of  the  surface  temperature  is  too  slow  to  keep  in  pace  with 
the  ambient  pressure  variation.  Since  the  phase-equilibrium  condition  is  posed  at 
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Figure  5.5.  Time  Variations  of  Characteristic  Times  of  Transient  Mechanisms 
Involved  in  Droplet  Vaporization,  n-Pentane/Nitrogen  System. 


the  droplet  surface,  the  variation  in  surface  composition  due  to  pressure  fluctua¬ 
tion  becomes  a  dominant  factor,  showing  an  amplitude  enhancement.  At  low  to 
moderate  frequency,  the  heat  transfer  from  the  ambient  is  used  either  to  gcisify 
the  liquid  fuel  or  to  heat  up  the  droplet  interior,  the  droplet  surface  temperature 
fluctuates  as  a  result  of  unsteady  heat  transfer  into  the  droplet.  However,  the  ther¬ 
mal  diffusion  and  relaxation  processes  in  the  droplet  interior  possesses  transient 
features  characterized  respectively  by  their  characteristic  times.  A  time  delay  for 
the  surface  temperature  to  response  to  the  pressure  fluctuation  is  presented.  As 
a  result,  the  surface  properties  (including  vapor  pressure,  species  concentration, 
and  latent  heat  of  vaporization)  which  depend  strongly  on  the  surface  tempera¬ 
ture  vary  with  pressure  fluctuation.  This  variation  changes  the  driving  force  for 


96 


p  =  50  atm;  I  pi  =  5  %  p  f  =  500  Hz 


Time,  ms 

Figure  5.6.  Effects  of  Frequency  of  Pressure  Fluctuation  on  Time  Histories  of 
Vaporization  Rate  Fluctuation,  n-Pentane/Nitrogen  System. 

gasification  process,  and  as  such  influences  the  dynamic  behavior  of  droplet  va¬ 
porization  responses.  The  time  histories  show  that  the  magnitude  of  the  response 
function  increases  first,  mainly  because  of  the  initial  heat-up  period  in  the  eaxly 
stage  of  vaporization.  The  magnitude  of  vaporization  response  reaches  a  maximum 
and  then  decreases.  It  is  attributed  to  the  decreases  of  the  characteristic  time  re¬ 
sulting  from  decreasing  of  droplet  radius,  and  as  such  reduces  the  magnitude  of 
vaporization  responses. 

Figure  5.7  depicts  the  effects  of  mean  pressure  on  droplet  vaporization  re¬ 
sponses.  The  ambient  mean  pressure  exerts  significant  influence  on  the  magnitude 
of  droplet  vaporization  response.  At  a  given  frequency  of  pressure  disturbance, 
the  amplitude  of  response  function  increases  with  increasing  pressure.  Although 
all  of  the  pressure-dependent  mechanisms  (e.g.,  liquid/vapor  phase  equilibrium, 
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Figure  5.7.  Effect  of  Mean  Pressure  on  Vaporization  Rate  Fluctuation,  n-Pentane/ 
Nitrogen  System;  /  =  3000  Hz. 

heat  transfer,  mass  and  thermal  diffusions)  involved  in  the  vaporization  process 
are  contributed  to  this  phenomenon,  one  major  contributing  factor  is  the  effect 
of  pressure  on  the  latent  heat  of  vaporization.  It  is  because  the  energy  required 
for  phase  transition  decreases  with  increasing  pressure.  As  a  result,  a  large  am¬ 
plification  gain  is  observed  at  high  pressure.  When  the  droplet  surface  reaches 
the  critical  mixing  state,  the  latent  heat  of  vaporization  drops  to  zero.  All  of  the 
heat  transfer  into  the  droplet  is  used  to  increase  the  liquid  temperature  and  to 
transform  the  liquid  phase  smoothly  into  the  vapor  phase.  As  a  consequence,  the 
vaporization  response  changes  significantly  once  the  critical  mixing  conditions  is 
reached. 

Figure  5.8  shows  the  time  history  of  the  vaporization  response  function 
for  p  =  80  atm.  The  vaporization  rate  fluctuates,  with  an  increasing  amplitude 
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Figure  5.8.  Time  Variation  of  Vaporization  Rate  Fluctuation,  n-Pentane/Nitrogen 
System;  p  =  80  atm,  /  =  3000  Hz. 

in  the  early  stage  of  the  droplet  lifetime.  Interestingly,  when  the  droplet  surface 
reaches  the  critical  point,  this  fluctuation  increases  abruptly  due  to  rapid  variations 
of  thermophysical  properties  of  the  droplet.  To  explain  this  phenomenon,  the 
temperature  dependence  of  latent  heat  of  vaporization  at  three  different  pressure 
(p  =  72,80,88  atm)  is  shown  in  Fig.  5.9.  If  the  ambient  pressure  executes  an 
isothermal  oscillation  with  a  10  %  decreases  in  amplitude  (from  point  A  to  point 
B),  the  latent  heat  of  vaporization  increases  by  50  KJ/Kg.  The  ensuring  influence 
on  the  vaporization  rate  is  obviously  profound.  On  the  other  hand,  if  the  ambient 
pressure  increases  by  10  %  in  amplitude  (from  point  A  to  point  C),  the  droplet 
surface  remains  critical  with  a  vanished  latent  heat  of  vaporization.  However,  the 
continuous  changes  in  the  critical  mixing  properties  with  pressure  modify  the  heat 
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Figure  5.9.  Temperature  Dependence  of  Latent  Heat  of  Vaporization  at  Three 
Different  Pressures,  n-Pentane/Nitrogen  System. 

transfer  and  mass  diffusion  processes  near  the  droplet  surface,  and  exert  strong 
influences  on  the  vaporization  behavior. 

Figure  5.10  shows  the  magnitudes  of  instantaneous  droplet  vaporization 
response  |i2p|  as  a  function  of  normalized  liquid  thermal  inertia  time  at  various 
ambient  pressures(i.e.,  10,  20,  and  50  atm).  Curves  a  -  b  plotted  in  solid  lines 
are  the  instantaneous  magnitudes  of  the  response  function  over  a  broad  range 
of  frequency  (/  =  250,500, 1000,2000, •••  ,9000  Hz,  respectively)  with  the  ini¬ 
tial  introductory  droplet  diameter  of  100//m.  Curves  c  -  d  drawn  in  dashed  line 
are  for  droplets  with  an  initial  diameter  of  300/Ltm  and  the  frequency  range  from 
10, 50, 100, 250, 500, 1000, 2000,  •  •  • ,  to  7000  Hz,  respectively.  For  all  the  cases  con¬ 
sidered  here,  the  amplitude  of  the  pressure  oscillation  remains  at  5  %  the  mean 
value.  Results  show  that  the  droplet  heat  up  process  which  occurs  in  the  early 
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Figure  5.10.  Magnitude  of  Droplet  Vaporization  Response  at  Various  Mean  Pres¬ 
sures,  n-Pentane/Nitrogen  System. 


stage  of  vaporization  has  a  strong  influence  on  the  magnitude  of  the  droplet  va¬ 
porization  response.  When  the  droplet  surface  approaches  the  quasi-steady  state, 
all  the  lines  seems  to  merge  together  for  both  large  (Do  =  300  /im)  and  small 
(Do  =  100  iim)  droplets. 

Figure  5.11  presents  the  phases  of  the  instantaneous  droplet  vaporization 
responses  0p.  The  notations  is  the  same  as  those  in  Fig.  5.10.  During  the  ini¬ 
tial  heat-up  period,  droplet  liquid  thermal  inertia  and  transient  diffusion  strongly 
influence  the  vaporization  rate,  showing  an  increases  of  phase  angle  as  surface 
temperature  increases.  The  phase  lag  increases  from  0°  toward  180°  as  the  fre¬ 
quency  increases,  due  to  the  large  deviation  of  the  characteristic  times  between 
pressure  oscillation  and  transient  vaporization  in  the  high  frequency  range.  One 
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Figure  5.11.  Phase  of  Droplet  Vaporization  Response  at  Various  Mean  Pressures, 
n-Pentane/Nitrogen  System. 


interesting  result  is  that  the  data  obtained  for  the  large  (Do  =  300  fim)  and  small 
(Do  =  100  fim)  droplets  essentially  follow  the  same  trend  closely. 

Figures  5.12(a)  and  (b)  show  the  effects  of  pressure  on  the  magnitude  and 
phase  angle  of  droplet  vaporization  response.  The  data  points  are  sampled  sequen¬ 
tially  throughout  the  droplet  lifetime  {t/ty  =  0.1,  0.3,  0.5,  0.7,  and  0.9).  The  filled 
symbols  depict  the  data  for  droplets  with  an  initial  diameter  of  100  //m,  while  the 
hollow  symbols  depict  the  cases  for  droplets  with  an  initial  diameter  of  300  //m. 
The  effect  of  initial  heatup  on  magnitude  and  phase  angle  are  compensated  using 
the  following  corrections 
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Figure  5.12.  Correlations  of  Droplet  Vaporization  Responses,  n-Pentane/Nitrogen 
System:  (a)  Magnitude,  (b)  Phase  Angle. 
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/T  X  —  T  \  0-65 

0;  =  9,-0.9(e,  +  .r)(£^_^*) 

where  To  is  the  droplet  averaged  temperature  and  T^e/  is  defined  as 


Tref  = 


< 


pseudo  wet— bulb  temperature 
critical  mixing  temperature 


if  P  <  Pc 
otherwise 


(5.2) 


where  the  pseudo  wet-bulb  and  critical  mixing  temperatures  can  be  found  in  Ta¬ 
ble.  4.2.  Results  for  the  initial  droplet  diameters  of  100  and  300  pm,  shown  re¬ 
spectively  by  the  filled  and  hollow  symbols,  scatter  along  the  same  curve  at  a  given 
pressure,  regardless  of  the  droplet  size.  The  amplitude  of  the  response  function 
increases  with  increasing  pressure  owing  to  the  susceptibility  of  enthalpy  of  vapor¬ 
ization  to  ambient  flow  oscillations  at  high  pressures.  On  the  other  hand,  the  effect 
of  mean  pressure  on  the  phase  angle  appears  quite  limited.  The  phase  decreases 
from  zero  in  the  low-frequency  limit  to  —180°  at  high  frequency,  a  phenomenon 
which  can  be  easily  explained  by  comparing  various  time  scales  associated  with 
fluid  transport  and  ambient  disturbance. 

In  order  to  decide  whether  the  droplet  vaporization  process  is  a  driving  or 
damping  mechanism  in  a  force  oscillatory  environment,  the  real  and  imaginary 
parts  of  vaporization  response  function  are  evaluated. 


Re[RpY  =  \Rp\*cos{e;) 

Im[Rp]*  =  |Rprsm(0;)  (5.3) 


The  results  are  shown  in  Figs.  5.13(a)  and  (b).  Based  on  the  dynamic  theory,  a 
positive  real  part  represents  the  driving  of  the  instability.  On  the  other  hand,  the 
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Figure  5.13.  Correlations  of  Droplet  Vaporization  Responses  Function,  n 
Pentane/  Nitrogen  System:  (a)  Real  Part,  (b)  Imaginary  Part. 
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negative  value  depicts  a  damping  mechanism.  The  peaks  occur  when  the  normal¬ 
ized  liquid  thermal  inertia  time  reaches  a  value  around  20.  The  frequency  response 
will  not  show  these  phenomena  unless  the  period  of  the  forcing  oscillation  becomes 
commensurate  with  some  characteristic  times  of  the  vaporization  mechanism.  This 
implies  that  the  strongest  coupling  occurs  between  pressure  fluctuation  and  vapor¬ 
ization  process.  The  high-frequency  cut-off  appears  when  the  value  reaches  200, 
where  the  vaporization  rate  response  switches  from  positive  to  negative,  changing 
the  driving  mechanism  of  combustion  instability  to  damping  effect. 

5.2  Approximate  Analyses 

In  parallel  to  the  numerical  calculation,  a  linear  approximate  analysis  based 
on  simplified  assumptions  is  developed.  The  purpose  is  to  assist  the  understanding 
of  the  effects  of  various  dynamic  mechanisms  on  droplet  vaporization  responses.  In 
addition,  the  non-dimensional  parameters  identified  in  the  analysis  are  related  to 
propellant  properties  that  may  be  used  to  characterize  and  examine  the  dynamic 
behavior  of  the  vaporization  process  for  any  other  propellants.  The  analysis  can 
also  be  used  to  validate  numerical  calculations,  especially  for  low-pressure  con¬ 
ditions.  We  consider  a  droplet  with  a  mass  of  m  and  a  mean  lifetime  of  fi,  in 
a  forced-oscillatory  environment.  Due  to  the  periodic  variation  of  the  ambient 
pressure,  the  properties  of  the  flowfield  (e.g.,  velocity,  pressure,  vaporization  rate, 
surface  temperature,  surface  vapor  pressure,  etc.)  also  feature  the  oscillatory  char¬ 
acteristics.  Assuming  the  deviations  of  these  properties  from  their  mean  values  are 
small  enough  and  can  be  approximated  linearly,  they  are  decomposed  into  a  time 
averaged  term  plus  a  small  fluctuating  value.  For  any  property  P,  the  definition  of 
a  small  perturbation  takes  the  form 
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1  _  _ 

P  =  (l  +  P')—  /  Pdt  =  P  +  PP'  (5.4) 

L^X  J t 

with  P'  the  fluctuation  percentage.  In  order  to  avoid  the  influence  of  disturbance 
into  the  mean  value,  P,  the  time  period  is  chosen  to  be  relatively  large  com¬ 
pared  with  the  period  of  fluctuation.  The  equations  required  in  the  analysis  are 
summarized  in  Table  5.2,  with  the  detailed  derivation  given  in  Appendix  C.  It  is 
worth  noting  that  the  transient  thermal-diffusion  process  in  the  droplet  interior  is 
modelled  using  two  different  approaches:  the  lumped  capacity  (infinite  liquid  ther¬ 
mal  conductivity,  X(  — *  oo)  and  the  finite  thermal  diffusion  (finite  liquid  thermal 
conductivity)  models. 

lumped  capacity  model; 

The  essence  of  the  lumped  capacity  method  is  based  on  the  assumption 
that  the  temperature  of  the  droplet  interior  is  spatially  uniform  at  any  instant 
during  the  transient  process,  implying  that  the  liquid  thermal  conductivity  within 
the  droplet  are  infinite.  The  transient  temperature  response  is  determined  by 
formulating  an  overall  energy  balance  for  the  droplet. 

=  Q,^,  (5.5) 

where  Cp/  and  Qs-^t  are  the  liquid  specific  heat  and  conduction  heat  flux  conducted 
into  the  droplet  interior,  respectively.  With  the  use  of  perturbation  technique,  the 
linearized  equation  Eq.  (5.5)  is  written  as 

=  <5„(<3L<  (5.6) 

The  droplet  vaporization  response  function  can  be  obtained  as 


TABLE  5.2 

Simplified  Analysis  of  Droplet  Vaporization  Response. 


107 


the  Pitzer's  acentric  factor. 
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-3^{1  -  iTyw)  {bdi +  '^r)  + 

Poo  [2$jj  -f  STyW'^j  +  b^]  +  i  [2^/  -  Z'^RT^w  -  ZbjdT^w 


where 


^  (toc  -Ts^ 

^  Cp^eTsTyW 

^  Ahr; 

j  _  /  dAhf!\  Poo 

'  “  KI^Jts'^v 

^  fdAh„\  Ts 
^  \  dTs  JvooAK 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


V  dTs  JpooAK  ^  ' 

This  transfer  function  represents  the  dynamic  relationship  between  the  vaporiza¬ 
tion  rate  and  ambient  pressure  oscillations.  When  a  sinusoidal  pressure  fluctuation, 
p'oo{t)  =  \p'\sin{uit),  is  superimposed  to  the  ambient  condition,  the  corresponding 
vaporization  rate  variation  can  be  written  as  w'(t)  =  |w'|sm(u;t  -|-  ©p),  where  |w'| 
and  0p  are  the  amplitude  and  phase  lag,  respectively,  and  take  the  form 


(5.12) 


\Rp\  =  m 


1  -1-  (r„u;)2 

[(Wl  +  ^ij)'  +  ^| 

\ 

r 

R  +  ZtvW'^  I  +  b^ 

2'^j-Z^rTvW 

—  Zb^T.o'W 

0p  =  ian  \-T,w)  +  tan 


—  tan'' 


'  2'$r  ZtvW^i  -k  6/3 

,2$/  —  Z^RTyW  —  3fe/Sr, 


— )  (5.13) 

TyW/ 
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thermal  diffusion  model: 

The  transient  one-dimensional  heat  conduction  equation  in  a  spherical  co¬ 
ordinate  system  is  solved  to  obtain  the  droplet  surface  temperature  response  to  a 
periodic  energy  flux. 

a  dt  dr  \  dr ) 

The  energy-flux  variation  at  the  droplet  surface  takes  the  form 

Qs^iQ's-.e  =  ii^rs\tTsF{w*)Tl  (5.15) 

where 


F(w*)  =  (iw*)^^^coih(iw*)^^^  —  1  (5.16) 

w*  =  B?wla  (5.17) 

The  droplet  vaporization  response  function  takes  the  form 

/  —S^{l  —  iT„w){bAhvyfdi+^R)  +  i^l 

i?,  =  —  =  - 1 - ^ J - y 

Poo  2'^R  +  StvW'^i  -f  b^AhvVf  -1-  i  2^/  —  St^w^r  —  dTvwbj3Ahv'w 

(5.18) 

where 

=  47rRXiTsRe[F{w*)]  +  AhJ-^^)+d2AKw  (5.19) 

=  47rR\iTsIm[F{w*)]  (5.20) 


The  magnitude  and  phase  angle  of  the  response  function  Rp  are 
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9^2 

1  {Tywf 

{bAhywdi  +  ’J'p)2  -i- 

2^p  +  +  b^Ahy-w 

2 

-F 

2^j  —  3t^u;^P  —  3Tywb0AhyW 

2 

(5.21) 


©,  =  (an-‘(-r.»)  + 

-  ^TyW'^R  -  STyWb^ AKw\ 

I  2^i? +  3t^u;^/ +  6/3A/i„w  / 

To  evaluate  Eqs.  (5.5)-(5.22),  a  knowledge  of  the  surface  properties  (Ts, 
A/it,,  ps,  and  A^)  is  required.  These  unknowns  can  be  solved  using 


„  _  Cp^sjToo  —  Tg)  _  Yf,oo  ^F,s 

-  Yf,s  -  1 


(5.23) 


in  conjunction  with 


Yp  =  —  - 

Ps  nsMF,s 


ppMp 

PsMf,s 


(5.24) 


where  the  partial  pressure  of  fuel  can  be  related  to  the  surface  temperature  through 
the  Clausius-Clapeyon  equation. 


dlnpjp  Ahv 

dT  ^  Kf! 

Three  unknowns(Yf’^s,  Ts,  and  pp)  are  involved  in  Eqs.  (5.23)-(5.25),  and  the  final 
solutions  are  obtained  through  an  iterative  procedure[58]. 

Since  the  equations  governing  the  droplet  vaporization  are  linearized  us¬ 
ing  the  small  perturbation  method,  a  block  diagram  derived  from  the  simplified 
analysis  is  shown  in  Fig.  5.14.  This  plot  illustrates  the  connections  among  the 
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Figure  5.14.  Block  Diagram  of  Approximate  Analysis. 


characteristic  processes  involved  in  droplet  vaporization.  Although  an  open-loop 
vaporization  response  function  to  ambient  pressure  fluctuation  is  under  consid¬ 
eration,  the  analysis  shows  that  a  closed-loop  mechanism  controls  the  dynamic 
behavior.  When  a  forced  pressure  oscillation  is  introduced  to  the  ambiance  of  a 
vaporizing  droplet,  the  amount  of  heat  conducted  to  the  droplet  interior  is  influ¬ 
enced  by  two  factors.  First,  the  heat  transfer  rate  varies  with  the  ambient  pressure, 
showing  an  enhancement  in  the  energy  transfer  from  the  gas  phase  once  the  pres¬ 
sure  raises  above  its  mean  value.  Second,  the  equilibrium  species  concentrations 
at  the  surface  is  sensitive  to  both  pressure  and  temperature.  This  can  be  further 
illustrated  by  the  thermodynamic  phase  equilibrium  analysis,  showing  that  the 
presence  of  a  second  species  at  the  droplet  surface  may  alter  the  latent  heat  of 
vaporization,  and  as  such  induces  a  variation  in  vaporization  rate[46]. 
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In  Fig.  5.12,  the  effect  of  pressure  on  the  magnitude  of  droplet  vaporization 
response  can  be  correlated  using 

\Rp\*  =  (— )°’Vpli>=10atm  (5.26) 

^Pref  f 

where  p^ef  is  10  atm.  Figures  5.15(a)  and  (b)  show  the  comparison  of  droplet  va¬ 
porization  response  between  correlations  from  numerical  calculations  (Eq.  (5.26)) 
and  approximate  analytical  data(Eqs.  (5.7)  and  (5.20))  at  p  =  1  atm.  The  dashed 
lines  depict  the  lumped  capacity  model,  and  dotted  line  the  finite  diffusion  model. 
The  results  obtained  from  simplified  analyses  deviate  substantially  from  the  nu¬ 
merical  data,  mainly  due  to  the  spatial  non-uniformity  of  thermal  properties  in  the 
droplet  interior.  To  compensate  this  effect,  a  correction  parameter  e  is  introduced 
into  the  lumped  capacity  model. 

8T 

=  Q,-t  (5.27) 

The  linearized  perturbation  equation  is  written  as 

6mCj,,efs^  =  Qs^iQ's-^t  (5-28) 

The  droplet  vaporization  response  function  can  be  obtained  as 


{bdi  +  +  i^j 

2Wjj  +  ZTyW^i  -1- 

+  i 

2^7  —  —  Zb^TyW 

where 


(5.29) 


'^R  = 

(too  -  r,  +  “*0 

(5.30) 

^7  = 

(6.31) 

Ahv 

The  amplitude  and  phase  lag,  respectively,  and  take  the  form 


Phase  angle,  0  Magnitude  of  R 


finite  thermal  diffusion  model 
lumped  capacity  model 
modified  lumped  capacity  model 


Figure  5.15.  Comparison  of  Droplet  Vaporization  Response  Function  with  Ap¬ 
proximate  Analyses,  n-Pentane/Nitrogen  System:  (a)  Magnitude,  (b)  Phase  An 
gle. 
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(5.33) 


k2'^j  —  S'^RTyW  —  3b0Ty 

For  the  n-pentane/iV2  system,  e  is  chosen  to  be  0.01.  Results  based  on  Eqs.  (5.32) 
and  (5.33)  are  plotted  in  solid  lines,  and  match  the  numerical  calculations  quite 
well. 


CHAPTER  6 


DROPLET  VAPORIZATION  IN  FORCED-CONVECTIVE 

ENVIRONMENTS 


The  gasification  characteristics  of  liquid  oxygen  (LOX)  droplet  in  a  super¬ 
critical  hydrogen  flow  is  studied  in  this  chapter.  The  initial  droplet  tempera¬ 
ture  and  diameter  are  100  K  and  100  fim^  respectively,  and  the  ambient  hydro¬ 
gen  temperature  is  1000  K.  In  order  to  examine  effects  of  the  convective  ve¬ 
locity  and  ambient  pressure  on  the  transient  behavior  of  the  droplet,  including 
interphase  mass,  momentum,  and  energy  transport,  a  parametric  study  is  con¬ 
ducted.  The  flow  conditions  cover  a  wide  range  of  pressure(p=100  —  400  atm)  and 
velocity(C/oo  =  0.1  — 15  mjs).  The  black  dots  in  Fig.  6.1  marks  the  scenarios  under 
consideration,  with  the  corresponding  Reynolds  number  shown  in  the  y-axis. 

Immediately  after  the  droplet  confronts  the  supercritical  hydrogen  flow, 
part  of  the  heat  transferred  from  the  gas  phase  goes  into  the  phase  transition  from 
liquid  to  vapor,  while  the  remainder  heats  up  of  the  droplet  interior.  In  contrast 
to  the  subcritical  droplet  vaporization,  no  thermodynamic  web-bulb  state  exists 
at  the  surface,  and  the  entire  gasification  process  becomes  transient  with  a  contin¬ 
uous  increase  of  the  surface  temperature.  In  addition,  because  the  critical  mixing 
temperature  of  the  oxygen /hydrogen  system  is  slightly  higher  than  the  droplet 
initial  temperature  and  decreases  with  increasing  pressure  (e.g.,  Tc  =  116  K  at 
p=400  atm),  the  thermodynamic  criticality  is  reached  at  the  droplet  surface  al¬ 
most  instantaneous  upon  its  introduction  to  the  flow.  Once  this  occurs,  the  latent 


116 


Figure  6.1.  Droplet  Reynolds  Number  as  a  Function  of  Freestream  Velocity  and 
Ambient  Pressure;  Do  =  100/xm,  T^o  =  lOOOiC,  To  =  lOOA. 


heat  of  vaporization  and  surface  tension  of  the  droplet  vanish,  rendering  essen¬ 
tially  a  continuous  medium  with  no  abrupt  phase  change  as  exists  for  a  subcritical 
droplet.  Although  the  interior  region  of  the  droplet  remains  at  a  liquid  state  with 
a  subcritical  temperature  distribution,  the  fluid  properties  and  their  gradients  vary 
continuously  between  the  liquid  core  and  the  ambient  gas.  Since  there  is  no  dis¬ 
tinct  boundary  existing  between  the  liquid  and  vapor  pha.ses,  for  convenience  of 
analysis,  the  surface  of  the  droplet  is  deliberately  defined  as  the  surface  at  which 
the  temperature  attains  the  value  of  the  critical  mixing  temperature. 

Several  snapshots  of  the  flowfields  under  various  conditions  are  presented 
in  Figs.  6.2-6.4.  Figure  6.2  shows  six  frames  of  isothermal  and  isoconcentra¬ 
tion  of  oxygen  at  a  convective  velocity  of  2.5  m/s  and  an  ambient  pressure  of  100 
atm.  The  free-stream  Reynolds  number.  Re,  is  30.56  based  on  the  initial  droplet 
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diameter.  The  evolution  of  the  temperature  field  exhibits  features  distinct  from 
that  of  the  concentration  field  due  to  the  disparate  time  scales  between  mass  and 
thermal  diffusion  processes  (Lewis  number  7^  1).  The  gasified  oxygen  is  convected 
downstream  from  its  edge.  Meanwhile,  the  liquid  core  with  large  momentum  in¬ 
ertia,  moves  slower  compared  with  the  gasified  oxygen.  The  oxygen  concentration 
contours  reveal  crescent  shapes  with  the  edge  bent  downstream.  At  t=1.08  ms, 
the  subcritical  liquid  core  disappears,  leaving  behind  a  puff  of  dense  oxygen  fluid 
which  is  convected  further  downstream  with  increasing  velocity  until  it  reaches  the 
momentum  equilibrium  with  the  ambient  hydrogen  flow. 

Figure  6.3  depicts  a  scenario  identical  to  that  of  Fig.  6.2,  but  with  an 
increased  convective  velocity  of  15  m/s.  The  droplet  evolution  is  substantially 
different,  exhibiting  several  distinct  modes  such  as  deformation,  viscous  stripping, 
and  secondary  breakup.  When  the  droplet  is  introduced  into  the  strong  convective 
flow  (Re=  183.36),  fast  vorticity  generation  resulting  from  the  velocity  difference 
between  the  liquid  core  and  ambient  flow  promotes  the  formation  of  an  attached 
eddy  behind  the  droplet.  This  recirculating  eddy  induces  a  slightly  higher  pres¬ 
sure  region  in  the  attached  wake,  and  exerts  dynamic  loading  over  the  rear  of  the 
droplet.  Meanwhile,  the  free  stream  also  exert  forces  to  the  front  surface.  As  a 
result,  the  droplet  deforms  and  extends  in  the  direction  normal  to  the  approaching 
flow.  The  flattened  edge  then  turns  to  the  streamwise  direction,  stretches  down¬ 
stream,  and  forms  a  skirt.  The  length  of  the  skirt  increases  with  time,  and  is 
controlled  by  the  balance  of  viscous  shear  force  between  the  dense  oxygen  and  the 
ambient  flow.  At  t=92  fis,  an  annulus  of  gaseous  oxygen  is  detached  from  the 
growing  skirt  due  to  the  local  flow  motion  and  volume  dilatation  resulting  from 
penetration  of  the  thermal  wave  into  the  skirt.  The  stripped  oxygen  from  the  edge 
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of  the  droplet  is  entrained  by  the  recirculating  flow,  forming  an  oxygen  tongue. 
The  dense  oxygen  influenced  by  the  eddy  recirculation  may  bend  toward  the  rear 
of  the  droplet.  This  process  carries  the  momentum  from  the  eddy  to  the  rear  of 
the  droplet  and  deforms  it  into  a  spherical  cape.  The  associated  cusp-like  rim  may 
catalyze  the  vorticity  generation  and  further  enhance  the  strength  of  the  recircu¬ 
lating  eddy.  Since  the  viscosity  of  the  droplet  interior  is  greatly  reduced  by  the 
energy  transfer  from  the  hot  stream,  the  liquid  oxygen  eventually  breaks  up  and 
forms  a  secondary  ring  around  its  core.  Another  signiflcant  phenomenon  is  the 
distribution  of  oxygen  concentration.  The  gasified  oxygen  stripped  from  the  edge 
of  the  droplet  leaves  the  recirculating  region  and  is  swept  downstream,  forming  a 
cylindrical  tail. 

Figure  6.4  presents  a  sequence  of  droplet  evolution  under  the  same  flow 
conditions  as  in  Fig.  6.2,  but  at  an  ambient  pressure  of  400  atm.  The  effect  of  pres¬ 
sure  on  droplet  dynamics  can  be  determined  by  a  direct  comparison  between  these 
two  figures.  Two  distinct  phenomena  are  observed.  First,  droplet  deformation  is 
enhanced  by  the  increase  in  ambient  pressure.  Here,  the  corresponding  Reynolds 
number  is  3.75  times  that  of  the  100  atm  case.  The  hydrogen  stream  carries  higher 
dynamic  momentum  and  consequently  exerts  stronger  force  on  the  droplet,  and 
as  such  promotes  its  deformation.  Second,  the  lifetime  of  the  droplet  decreases 
with  increasing  ambient  pressure  mainly  because  that  the  higher  Reynolds  num¬ 
ber  at  higher  pressure  increases  the  convective  heat  transfer  rate  and  as  a  result 
accelerates  evaporization. 

Figures  6.5  and  6.6  show  the  evolution  of  droplet  critical  surfaces.  The 
solid  lines  represent  the  instantaneous  isotherms  of  the  critical  mixing  tempera¬ 
ture,  while  the  dashed  lines  mark  the  surfaces  of  critical  mixing  composition.  For  a 


Figure  6.5.  Effects  of  Ambient  Velocity  on  Evolution  of  Droplet  Critical  Surfaces,  LOXf  H2  System  at  p  =  100  atm. 


Ambient  Pressure  on  Evolution  of  Droplet  Critical  Surfaces 
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binary  system  of  oxygen  and  hydrogen,  the  critical  mixing  temperature  and  com¬ 
position  can  be  obtained  through  a  phase  equilibrium  analysis  and  are  given  in  the 
Table  6.1.  Because  of  the  non-unity  Lewis  number  (the  ambient  Lewis  number  is 


TABLE  6.1 

Critical  Mixing  Properties  of  O2I H2  System. 


pressure 

critical  mixing 

critical  mixing 

(atm) 

temperature  (k) 

composition 

100 

142.79 

0.735 

200 

127.21 

0.561 

400 

116.0 

0.496 

about  2.0  at  100  atm,  and  increases  toward  the  droplet  surface),  the  critical  mixing 
isotherm  propagates  at  a  speed  different  from  the  surface  of  critical  mixing  com¬ 
position.  The  critical  isotherm  disappears  before  the  ambient  hydrogen  penetrates 
into  the  droplet  center.  Figure  6.5  depicts  the  effects  of  convective  velocity  on  the 
evolution  of  droplet  critical  surface.  The  dynamic  deformation  of  the  droplet  is 
substantially  enhanced  by  increasing  the  momentum  carried  by  the  ambient  flow. 
Figure  6.6  illustrates  the  effects  of  pressure  on  the  evolution  of  the  critical  surfaces. 
Entrainment  of  gasified  oxygen  into  the  attached  recirculation  eddy  seems  to  be 
augmented  with  increasing  pressure.  At  high  pressures,  the  motion  of  the  skirt 
even  reveals  a  periodic  pattern.  The  gaseous  oxygen  entrapped  by  the  recirculat¬ 
ing  flow  tends  to  move  with  the  strong  eddy  and  may  run  forward  into  the  rear  of 
the  droplet.  This  process  drives  the  skirt  to  expand  in  the  cross  stream  direction. 
Yet,  the  convective  flow  suppresses  the  expansion  and  forces  the  skirt  to  bend  and 


stretch  downstream.  As  a  result,  more  oxygen  is  trapped  into  the  vortex  eddy  and 
therefore,  starts  another  locomotion  again. 

Figures  6.7(a)-(d)  show  the  streamline  patterns  at  different  flow  conditions. 
Four  distinct  modes  are  observed,  namely  spherical,  deformation,  stripping,  and 
breakup  modes.  Unlike  low-pressure  cases,  no  discernible  flow  recirculation  takes 
place  in  the  droplet  interior,  regardless  of  the  Reynolds  number.  This  may  be 
attributed  to  the  diminished  surface  tension  at  the  droplet  surface  at  supercritical 
conditions.  In  addition,  the  blowing  effect  induced  by  droplet  vaporization  in¬ 
creases  the  boundary  layer  thickness  and  obstructs  the  formation  of  a  shear  vortex 
in  the  droplet  interior.  The  rapid  change  of  the  droplet  shape  precludes  a  stable 
shear  stress  field  in  the  liquid  core  which  in  turn  contributes  to  the  absence  of 
liquid  recirculation. 

For  all  the  cases  considered  here,  the  convective  stream  is  from  the  left  to 
right,  and  a  stagnation  point  where  the  velocity  induced  by  gasification  is  cancelled 
by  the  approaching  flow  can  be  located  on  the  left  side  of  the  droplet.  Stagna¬ 
tion  streamlines  originated  from  the  stagnation  point  divide  the  flowfield  into  two 
regions:  internal  and  external  flows.  The  spherical  mode  shown  in  Fig.  6.7(a)  is 
typically  observed  when  the  droplet  is  introduced  into  an  environment  with  a  very 
low  Reynolds  number  (Re=2.5).  Although  flow  separation  is  encouraged  by  va¬ 
porization  blowing,  no  recirculation  eddy  is  found  in  the  wake  behind  the  droplet. 
This  is  because  the  vorticity  generated  from  the  shear  flow  is  too  weak  to  form  any 
confined  eddy.  Further  increase  of  the  convective  Reynolds  number  may  cause  the 
formation  of  recirculation  eddy,  as  shown  in  Fig.  6.7(b).  This  plot  is  taken  from  the 
LOX  droplet  vaporization  in  a  convective  hydrogen  flow  with  an  initial  Reynolds 
number  of  18.34  (p=100  atm,  Uoo  =  1-5  m/s).  At  t=0.61  ms,  the  droplet  deforms 


et  Vaporization  in  Supercritical  Hydrogen  Flow  at  p  =  100  atm:  (a)  Spherical 
Deformation  Mode;  f/oo  =  1-5  m/s,  t=0.61  ms.  (c)  Stripping  Mode;  Uoo  =  5  mj 


into  an  olive  shape  with  a  spherical-like  recirculation  ring  attached  behind  it.  It 
is  worth  noting  that  the  droplet  deformation  and  surface  blowing  effect  tend  to 
elevate  the  threshold  Reynolds  number  above  which  the  recirculation  eddy  forms. 
This  result  is  in  close  agreement  with  numerical  solutions  by  other  researchers 
(Dennis  et  al.  1971;  LeClair  1970;  LeClair  et  al.  1970;  Masliyah  1970;  Woo  1971). 
Figure  6.7(c)  depicts  the  flow  structure  with  viscous  stripping,  showing  an  oblate 
droplet  with  a  deformed  vortex  ring.  Lane  (1951)  has  shown  experimentally  that, 
and  Hinze  (1955)  has  confirmed  theoretically  when  a  droplet  is  suddenly  exposed 
to  a  high-velocity  gas  stream,  the  droplet  aspect  ratio  increases  with  time.  The 
flattened  edge  of  the  droplet  enhances  the  strength  of  the  recirculating  eddies  and 
as  such  increases  the  viscous  shear  stress  dramatically.  As  a  consequence,  a  thin 
sheet  of  dense  oxygen  may  be  stripped  from  the  edge  of  the  droplet  and  swept  to¬ 
ward  the  outer  boundary  of  the  recirculation  eddy.  The  secondary  breakup  mode 
is  clearly  illustrated  in  Fig.  6.7(d).  At  Re  =  183.36,  hydrogen  penetrates  through 
the  liquid  phase,  breaking  the  droplet  into  two  parts:  the  core  disk  and  the  sur¬ 
rounding  ring.  The  flow  structure  also  indicates  that  the  volume  of  the  attached 
vortex  ring  increases  significantly. 

Figure  6.8(a)  presents  the  effects  of  free-stream  velocity  on  the  instanta¬ 
neous  variation  of  droplet  mass  at  p=100  atm.  The  droplet  residual  mass  is  de¬ 
fined  herein  as  the  mass  confined  by  the  isothermal  surface  at  the  critical  mixing 
temperature.  In  a  slow  convective  stream  {Uoo  =  0.2  m/s),  heat  conduction  plays 
an  important  role  during  the  vaporization  process.  The  rapid  droplet  gasifica¬ 
tion  in  the  early  stage  of  its  lifetime  results  from  the  large  temperature  gradient 
near  the  liquid/vapor  interface.  After  the  initial  heatup  period,  the  blowing  from 
the  droplet  surface  hinders  the  thermal  wave  propagation  and  results  in  a  slower 


(a) 


Figure  6.8.  Time  Variations  of  Droplet  Residual  Mass  at  Various  Convective 
Velocities,  LOXIH2  System:  (a)  Effect  of  Convective  Velocity  at  p  =  100  atm. 
(b)  EflFect  of  Ambient  Pressure  at  Uca  =  5-0  m/s. 
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Figure  6.9.  Dependences  of  Droplet  Lifetime  on  Freestream  Velocity  and  Ambient 
Pressure,  LOXIH2  System. 


gasification  rate.  In  a  high  Reynolds-number  flow,  increased  free-stream  velocity 
enhances  the  convective  heat  transfer,  and  as  an  outcome,  the  droplet  vaporizes 
faster.  Additionally,  droplet  deformation  resulting  from  stronger  dynamic  loading 
from  the  approaching  flow  also  encourages  the  vaporization  processes.  The  effects 
of  pressure  on  the  droplet  vaporization  behavior  are  illustrated  in  Fig.  6.8(b).  As 
the  pressure  increases,  the  convective  heat  transfer  which  depends  strongly  on  the 
Reynolds  number  increases  since  the  Reynolds  number  also  increases  with  pressure. 
In  addition,  the  critical  mixing  temperature  used  to  define  the  droplet  surface  de¬ 
creases  with  increasing  pressure.  As  a  result,  the  droplet  gasification  rate  increases 
with  the  pressure. 

Figure  6.9  presents  the  droplet  lifetime  as  a  function  of  ambient  pressure 
and  velocity.  When  a  droplet  is  vaporizing  in  a  quiescent  environment  (Re=0), 
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Figure  6.10.  Correlation  of  Droplet  Lifetime  as  a  Function  of  Free-Stream 
Reynolds  Number  and  Pressure. 

i 

I 

the  primary  mode  of  energy  transfer  is  heat  conduction.  The  gasification  rate  is 
enhanced  by  convection,  when  a  relative  motion  is  present  between  the  droplet  and 
surrounding  gases.  A  correlation  is  developed  for  the  convective  correction  to  the 
droplet  lifetime  as  a  function  of  ambient  pressure  and  droplet  Reynolds  number. 
The  correlation  is  given  as 


Tf,Re==0  l  +  0.165634Reii(Pr,O2)-°-®® 

where  Re  is  the  Reynolds  number  based  on  the  initial  droplet  diameter  and  Pt,02 
the  reduced  pressure  in  reference  to  the  critical  pressure  of  oxygen.  The  comparison 
of  the  correlation  with  the  numerical  data  is  shown  in  Fig.  6.10.  Equation  (6.1) 
bears  a  resemblance  to  the  popular  Ranz  and  Marshall  correlation  [59]  for  droplet 
heat  transfer  correction  due  to  convective  effect,  which  takes  the  form 
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1  +  0.3/?ei/2pri/3 


(6.2) 


However,  the  Ranz  and  Marshall  correlation  is  only  applicable  to  low  pressure 
flows,  and  shows  a  weaker  Reynolds-number  dependency. 

Figures  6.11(a)  and  (b)  show  the  location  of  droplet  center  of  gravity  at 
different  free-stream  velocities  and  pressures.  Because  there  is  no  distinct  va¬ 
por/liquid  interface  for  a  supercritical  droplet,  droplet  motion  is  best  characterized 
by  the  movement  of  the  droplet  center  of  gravity  which  can  be  located  by  dividing 
the  first  moment  of  the  droplet  mass  by  its  total  mass  confined  by  the  critical 
surface.  Results  show  that  droplet  travel  distance  depends  on  both  the  drag  force 
and  the  total  lifetime  over  which  the  force  acts  on  the  droplet.  For  instance,  at 
p=100  atm,  even  though  the  convective  Reynolds  number  at  Uoo  =  5  m/s  is  larger 
than  that  at  Uoo  =  1.5  m/s,  the  droplet  move  further  for  the  low  velocity  case, 
simply  due  to  its  longer  lifetime. 

The  acceleration  of  the  droplet  center  of  gravity  can  be  determined  by  tak¬ 
ing  the  time  derivative  of  droplet  velocity,  which  is  obtained  by  dividing  the  total 
momentum  of  liquid  droplet  by  its  mass.  The  results  show  that  the  droplet  ac¬ 
celerates  faster  with  increasing  Reynolds  number.  At  a  given  ambient  velocity, 
momentum  transfer  between  the  ambient  flow  and  droplet  is  also  encouraged  by 
increasing  the  ambient  pressure.  By  applying  a  least-square  curve  fitting  to  these 
data,  the  droplet  acceleration  is  expressed  by  a  first  order  polynomial  of  time. 
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Figure  6.11.  Time  Variations  of  Location  of  Droplet  Center  of  Gravity,  LOXj H'z 
System:  (a)  Effects  of  Convective  Velocity  at  p  =  100  atm.  (b)  Effects  of  Ambient 
Pressure  at  Uoo  =  2.5  mis. 
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a  —  aa  -{■  Pat  (^•3) 

where  eta  and  Pa  are  constants  given  in  Figs.  6.12(a)  and  (b),  being  functions  of 
initial  Reynolds  number. 

The  motion  of  a  droplet  is  essentially  a  result  of  forces  acting  on  its  surface, 
including  form  drag,  aerodynamic  normal  force,  and  viscous  shear  force.  Since  the 
droplet  motion  is  tracked  by  its  center  of  gravity,  the  drag  force  is  best  evaluated  by 
the  net  momentum  transfer  to  the  droplet.  Figure  6.13  presents  the  time  histories 
of  the  total  drag  force,  aerodynamic  loading,  and  drag  coefficient  at  p  =  100  atm 
and  Uoo  —  5.0.  The  drag  coefficient  is  defined  as  the  ratio  of  the  drag  force  to  the 
aerodynamic  normal  force. 
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where  {Uoo  —  Udrop)  is  the  relative  convective  velocity  and  Rc  is  the  radius  of  the 
mass-equivalent  sphere  with  a  uniform  density  distribution  of  its  initial  value.  The 
drag  force  acting  on  the  droplet  first  increases  with  time  due  to  the  large  form 
drag  arising  from  the  high  aspect  ratio.  It  then  drops  to  zero  at  the  end  of  droplet 
lifetime. 

Chen  and  Yuen[60]  found  that  the  drag  coefficients  of  evaporating  droplets 
are  small  compared  to  those  of  a  non-vaporizing  solid  sphere  at  the  same  Reynolds 
numbers.  Several  researchers  [61,62]  numerically  analyzed  evaporating  droplet  mo¬ 
tion  by  solving  the  Navier- Stokes  equations  and  proposed  the  following  correlation 


P>  m/s  ct.  m/s 


Reynolds  Number,  Remj 


Figure  6.12.  Dependences  of  Acceleration  Coefficients  on  Free-stream  Reynolds 
Number  and  Ambient  Pressure,  LOXIH2  System:  (a)  Droplet  Acceleration  at 
i  =  0  sec.  (b)  Droplet  Acceleration  Rate. 
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Time,  ms 

Figure  6.13.  Time  Variations  of  Droplet  Total  Drag  Force,  Aerodynamic  Loading, 
and  Drag  Coefficient,  LOX/H2  System;  p=100  atm,  Foo  =  5  m/s. 


=  ilfw 

where  denotes  the  drag  coefficient  for  non- vaporizing  hard  sphere,  and  6  a  con¬ 
stant  which  has  a  value  of  0.2  for  Renksitzbulut’s  model[61]  and  0.32  for  Chiang’s 
model[62].  A  transfer  number  B  is  adopted  to  account  for  the  effect  of  blowing  on 
momentum  transfer  to  the  droplet.  For  droplet  vaporization  at  low  to  moderate 
pressures  {pr  <  0.5),  Spalding  transfer  number  is  widely  used  to  characterize  the 
vaporization  rate. 


Cp{T^  -  Ts) 


with  Cp,  Too,  and  Ts  the  constant-pressure  specific  heat,  ambient  temperature,  and 
droplet  surface  temperature,  respectively.  The  variable  Ahv  represents  the  latent 
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Reynolds  Number 


Figure  6.14.  Effects  of  Ambient  Reynolds  Number,  Pressure,  and  Vaporization 
Number  on  Droplet  Drag  Coefficient. 

heat  of  vaporization  which  becomes  zero  at  critical  point,  rendering  an  infinite 
value  for  the  transfer  number. 

Data  points  used  to  obtain  the  correlation  are  sampled  sequentially  through¬ 
out  the  droplet  lifetime  for  all  the  scenarios  under  consideration.  To  quantify  the 
droplet  vaporization  rate,  a  transfer  number  Bx,  which  is  suited  for  supercritical 
droplet  vaporization,  is  proposed  as 

where  the  quantities  Tg  and  denote  the  critical  mixing  temperature  and  the 
averaged  temperature  of  droplet,  respectively. 

The  drag  coefficient  can  be  correlated  as  a  function  of  Reynolds  number, 
ambient  pressure,  and  vaporization  transfer  number,  as  shown  in  Fig.  6.14. 
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Figure  6.15.  Correlation  for  Droplet  Drag  Coefficient. 


with  a  and  b  selected  to  be  0.05  and  1.592  (pr,02)  the  reference  drag 

coefficient,  for  a  hard  sphere[63]  and  takes  the  form 


The  data  points  cluster  along  the  drag  curve  in  the  low  Reynolds-number  region. 
However,  deviation  from  the  reference  drag  curve  defined  in  Eq.  (6.9)  occurs  at  high 
Reynolds  numbers  (Re  >  10).  This  is  mainly  attributed  to  the  increase  in  form 
drag  which  results  from  increasing  droplet  aspect  ratio  due  to  strong  convective 
stream.  Although  the  shape  aspect  ratio  serves  as  a  good  parameter  to  correlate 
the  data,  it  is  not  practical  to  use  in  spray  combustion  analyses.  A  correction  factor 
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CHAPTER  7 


CONCLUING  REMARKS 

A  systematic  investigation  has  been  conducted  to  study  supercritical  droplet 
gasification  and  combustion  in  both  quiescent  and  forced-convective  environments. 
The  formulation  is  based  on  the  time-dependent  conservation  equations  of  mass, 
momentum,  energy,  and  species  concentration  in  axisymmetric  coordinates  for 
both  the  droplet  interior  and  ambient  gases.  Full  account  is  taken  of  thermody¬ 
namic  non-ideality  and  transport  anomaly  during  the  transition  of  droplet  surface 
from  the  subcritical  to  supercritical  state.  In  addition,  a  unified  property  evalu¬ 
ation  scheme  based  on  the  Benedict- Webb-Rubin  (BWR)  equation  of  state  and 
the  extended  corresponding-state  principle  is  established  to  predict  fluid  thermo¬ 
physical  properties  over  the  entire  thermodynamic  regime.  The  governing  equa¬ 
tions  and  associated  boundary  conditions  are  solved  numerically  using  an  implicit 
finite- volume  scheme  with  a  dual  time-stepping  integration  technique. 

A  series  of  calculations  have  been  carried  out  to  study  the  behavior  of  both 
hydrocarbon  and  cryogenic  liquid  droplets  over  a  wide  range  of  thermodynamic 
state,  including  both  sub-  and  super-critical  conditions.  The  vaporization  and 
combustion  of  an  isolated  liquid  droplet  in  a  quiescent  environment  are  first  in¬ 
vestigated  in  Chapter  4.  This  configuration  allows  attention  to  be  focused  on 
the  effects  of  ambient  thermodynamic  states  on  the  interphase  transport  between 
droplet  and  the  surroundings.  Results  show  that  the  ambient  gas  pressure  ex¬ 
erts  significant  control  of  droplet  gasification  and  burning  processes  through  its 
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influence  on  fluid  transport,  gas-liquid  interfacial  thermodynamics,  and  chemi¬ 
cal  reactions.  The  droplet  gasification  rate  increases  progressively  with  pressure. 
However,  the  data  for  the  overall  burnout  time  exhibits  a  considerable  change  of 
combustion  mechanism  at  the  critical  pressure,  mainly  as  a  result  of  reduced  mass 
diffusivity  and  latent  heat  of  vaporization  with  increased  pressure.  The  influence 
of  droplet  size  on  burnout  characteristics  is  also  noted. 

The  dynamic  responses  of  droplet  vaporization  to  the  ambient  flow  oscil¬ 
lations  are  presented  in  Chapter  5.  The  analysis  extends  the  work  of  droplet 
vaporization  in  quiescent  environments  and  imposes  a  periodic  pressure  oscillation 
to  ambient  flow.  The  purpose  is  to  assess  the  effect  of  flow  oscillations  on  va¬ 
porization  processes  as  a  function  of  frequency  and  amplitude  of  the  fluctuation, 
as  well  as  the  type  of  oscillation.  Representative  environmental  conditions  are 
simulated.  Effects  of  ambient  pressure  and  droplet  heatup  on  the  magnitude  and 
phase  angle  of  the  vaporization  response  are  quantified  by  means  of  numerical  cal¬ 
culations.  The  amplitude  of  response  function  increases  with  increasing  pressure 
owing  to  the  susceptibility  of  enthalpy  of  vaporization  to  ambient  flow  oscillations 
at  high  pressures.  On  the  other  hand,  the  effect  of  mean  pressure  on  the  phase 
angle  appears  quite  limited.  The  phase  decreases  from  zero  in  the  low-frequency 
limit  to  —180°  at  high  frequency,  a  phenomenon  which  can  be  easily  explained  by 
comparing  various  time  scales  associated  with  fluid  transport  and  ambient  distur¬ 
bance.  When  the  droplet  surface  approaches  the  thermodynamic  critical  mixing 
state,  an  abnormal  amplification  of  vaporization  rate  caused  by  rapid  variations  of 
thermophysical  properties  is  clearly  observed. 

Vaporization  of  cryogenic  fluid  droplets  in  supercritical  streams  is  also  in¬ 
vestigated  systematically.  A  series  of  calculations  have  been  carried  out  to  study 
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the  gasification  of  an  isolated  oxygen  (LOX)  droplet  in  a  hydrogen  stream  over 
a  broad  range  of  pressure  (100-400  atm)  and  Reynolds  number  (20-300).  Results 
indicate  that  forced  convection  exerts  profound  influence  on  droplet  dynamics  and 
gasification  characteristics.  The  droplet  lifetime  decreases  with  increasing  ambient 
Reynolds  number,  mainly  due  to  increased  momentum  transfer  in  the  boundary 
layer  which  intensifies  the  convective  heat  transfer  into  the  droplet  and  conse¬ 
quently  expedites  the  gasification  process.  Furthermore,  as  a  result  of  stronger 
dynamic  loading  on  the  droplet,  high  Reynolds-number  flows  (Re  >  200)  tend  to 
facilitate  droplet  breakup  and  shattering  as  a  result  of  stronger  dynamic  loading  on 
the  droplet.  In  all  the  cases  considered  here,  no  discernible  flow  recirculation  takes 
place  in  the  droplet  interior,  regardless  of  the  Reynolds  number.  This  may  be  at¬ 
tributed  to  the  diminished  surface  tension  at  supercritical  conditions.  The  blowing 
effect  induced  by  droplet  vaporization  increases  the  boundary  layer  thickness  and 
obstructs  the  formation  of  a  shear  vortex  in  the  droplet  interior.  The  short  local 
residence  time  of  the  droplet  surface  layer  arising  from  the  large  droplet  gasification 
rate  at  high  pressure  also  contributes  to  this  phenomenon. 

Effects  of  ambient  conditions  (temperature,  pressure,  and  convective  veloc¬ 
ity)  on  droplet  gasification  lifetime  and  drag  coefficient  are  investigated  paramet¬ 
rically.  Data  obtained  are  correlated  as  functions  of  ambient  pressure,  Reynolds 
number,  and  vaporization  transfer  number. 


Tf,Re=0  1-f  0.165634i?ei-i(Pr,O2)-°-®® 

where  Re  and  Pt,02  refer  to  the  Reynolds  number  based  on  the  initial  droplet 
diameter  and  the  reduced  pressure  by  the  critical  pressure  of  oxygen. 


142 


Cd  = 


C%Re 


0.3 


(1  +  aBx) 


1.592 


(Pr,02) 


-0.7 


(7.2) 


with 


- 

Ud  - 


24 

R^[ 


1  ^  2 

1  +  -Re^ 

D 


(7.3) 


Although  the  current  research  has  provided  significant  information  to  under¬ 
standing  the  underlying  mechanisms  involved  in  supercritical  droplet  vaporization 
and  combustion,  a  systematic  verification  of  these  models  is  still  lacking.  Future 
work  in  supercritical  droplet  modelling  should  emphasize  model  validation  against 
experimental  data. 
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APPENDIX  A 


THERMODYNAMIC  PROPERTIES 

Prediction  of  thermodynamic  properties  plays  a  very  important  role  in  ana¬ 
lyzing  high-pressure  droplet  behavior.  From  a  microscopic  viewpoint,  intermolec- 
ular  mean  free  path  tends  to  decrease  with  increasing  pressure,  and  the  molecular 
size  effect  as  well  as  intermolecular  force  are  enhanced  and  cause  the  gas  and  liq¬ 
uid  behavior  deviates  away  from  the  ideal  fluid.  Property  evaluation  in  the  near 
critical  regime  also  brings  up  another  challenge.  When  approaching  the  thermody¬ 
namic  critical  point,  fluid  exhibits  anomalous  behavior  and  are  extremely  sensitive 
to  pressure  and  temperature.  For  a  multicomponent  system,  change  of  properties 
due  to  mixing  is  even  more  complicated  which  ideal  gas  theory  is  not  sufficient  to 
predict  thermodynamic  non-idealities  at  high  pressure. 

In  order  to  establish  a  unified  evaluation  scheme  for  real-gas  properties, 
a  pseudo-pure  substance  model  of  real  fluid  mixture  is  adopted.  This  approach 
treats  a  mixture  as  a  pseudo-pure  substance  having  its  own  set  of  properties.  The 
difference  between  these  properties  and  those  at  a  reference  state  at  the  same  tem¬ 
perature  are  represented  by  thermodynamic  departure  functions  which  are  func¬ 
tions  of  p-V-T  properties  of  the  fluid.  A  basic  prerequisite  of  this  treatment  is 
the  concept  and  use  of  appropriate  reference  state.  Generally,  a  thermodynamic 
state  which  reveals  ideal  gas  behavior  is  preferred.  For  an  ideal  gas,  the  internal 
energy  is  a  function  of  temperature  only.  This  means  that  an  ideal  gas  at  a  given 
temperature  has  certain  definite  thermodynamic  properties  (i.e.,  specific  internal 


energy,  specific  enthalpy,  fugacity,  etc.),  regardless  of  the  pressure.  Under  these 
conditions,  the  departure  function  becomes  a  high-density  correction  and  can  be 
evaluated  accurately  with  the  help  of  a  precise  equation  of  state.  For  an  ideal  mix¬ 
ture,  change  of  properties  on  mixing  is  also  negligibly  small.  Only  using  a  simple 
mixing  rule  is  required  to  accurately  predict  the  mixture  properties. 

specific  enthalpy 

Since  the  SRK  and  BWR  equations  of  state  are  both  pressure  explicit  in 
nature,  i.e.,  pressure  is  expressible  as  a  function  of  temperature  and  density,  the 
evaluation  of  thermodynamic  properties  starts  with  the  definition  of  the  Helmholtz 
function,  a 

a  =  e  —  Ts 

da  =  de  —  Tds  —  sdT  (A.l) 

Since  e  —  h  —  pj Eq.  (A.l)  yields 

dh^da-^Tds  +  sdT +  (A.2) 

This  can  be  further  illustrated  by  integrating  along  an  isotherm  from  a  reference 
state  (where  the  ideal  gas  behavior  can  be  assumed)  to  the  state  of  interest.  With 
the  ideal  gas  equation  of  state  {p  =  p'RT),  Eq.  (A.2)  becomes 

h  =  h°+\{a-  a°)  +  T{s  -  s°)  +  nT{Z  -  1)  (A.3) 

where  Z  is  the  compressibility  and  the  superscript  0  represents  the  values  at  the 
ideal  state.  The  second  term  of  Eq.  (A.3)  is  the  density  correction  for  the  Helmholtz 
free  energy.  Since 
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dh  =  Tds  + 


dp 


(A.4) 


Equation  (A. 2)  can  be  written  as 

da=(^^dp-sdT  (A.5) 

Along  an  isotherm,  Eq.  (A.5)  is  integrated  from  an  ideal  state  with  a  density  of 
to  a  thermodynamic  state  of  interest  with  a  density  of  p  to  form 

where  is  chosen  with  which  the  gas  reveals  the  ideal  gas  behavior  at  1  atm. 
Following  Maxwell’s  equation,  the  entropy  can  be  expressed  as  the  first  derivative 
of  the  Helmholtz  free  energy  with  respect  to  temperature  at  a  constant  density 


=-( 


da\ 

df), 


(A.7) 


Substituting  (A. 6)  into  (A.7)  to  obtain  the  departure  function  of  entropy  as 


5  =  - 


=  S 


fp  1  ^dp\ 

\dT)f  . 

JpO  p^KdTJ 

/;?(*).* 


(A.8) 


/p^  p  \UX  /  p 

Enthalpy  of  a  mixture  can  be  obtained  by  substituting  Eqs.  (A. 6)  and  (A.8)  into 
(A. 3)  to  get 

where  is  the  ideal-gas  enthalpy  which  is  a  function  of  temperature  only.  If  the 
constant-pressure  specific  heats  of  all  components  are  given,  the  specific  enthalpy 
of  an  ideal  mixture  can  be  expressed  as 
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N 


(A.10) 


A"  +  1  Cl.dT 

z=l  ^  ■^-‘re/ 

where  h^gfi  and  Yi  are  the  enthalpy  of  formation  at  a  reference  temperature  of 
Tref  and  the  mass  fraction  of  i  component,  respectively. 


specific  internal  energy 

Based  on  the  definition  of  specific  enthalpy, 


or 


.  P 

e  +  - 
P 


dh  —  dc  di 


(A.ll) 


(A.12) 


the  specific  internal  energy  can  be  obtained  by  integrating  along  an  isotherm  from 
a  reference  state  where  the  ideal  gas  behavior  can  be  assumed  to  the  state  of 
interest. 


Substituting  Eqs.  (A. 9)  and  (A.  15)  into  (A.  13)  to  obtain 


(A.13) 


e=  + 


P  1 

pO  p2  L 


T 


I'A'i 


dp 


(A.  14) 


The  first  term  on  the  right  hand  side  is  the  specific  internal  energy  of  an  ideal  gas 
which  is  a  function  of  temperature  only,  and  can  be  expressed  as 


Er.k./.i+4  cl,dT 
.=1  l  J 


(A.15) 


specific  heat 


Two  parallel  expressions  can  be  found  for  the  change  of  entropy.  First,  if 
entropy  is  treated  as  a  function  of  temperature  and  pressure 

s  =  s{T,p) 


Therefore, 


(A.16) 

From  Maxwell’s  equations, 

p’^xdTJp  xdp/T 

(A.17) 

(A.18) 

Substituting  Eqs.  (A. 17)  and  (A. 18)  into  (A. 16)  to  get 

^  ^  dT  1  (dp\  , 

+  p[wr) 

(A.19) 

Secondly,  entropy  can  also  be  expressed  as  a  function  of  temperature  and  density 
as 


s  ==  s{T,p) 


Therefore, 


*  =  (ot. 

)  dT  + 

P 

(A.20) 

Similarly,  from  Maxwell’s  equations, 

= 

(de\ 

KdTJp 

^[dTjp 

(—)  = 
\dTJp 

.(ds 
^  \df 

))t 

Equation  (A. 20)  can  be  expressed  as 
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,  _  dT  \(dp\ 


Combine  Eq.  (A. 19)  and  (A. 21)  to  obtain 


{A.21) 


dr  1  fdp\  ,  _„dT  1  fdp\ 

fT\\dTjp,  fT\\dTjp^ 


p 


(A.22) 


Temperature  can  be  expressed  as  a  function  of  pressure  and  specific  volume 


T  =  f{p,  p) 


Therefore, 


Compare  (A.22)  with  (A. 23)  to  obtain 


(A.23) 


dp)  p 
dp)  p 


(—) 

T\  KdTjp 

pycp-c^ 

(^) 

T\  \dTjp 

p^)Cp-Cv 


When  these  equations  are  solved  for  Cp  —  Cv,  they  yield  the  same  result  as 


Cp-Cv 


py\dT 


(A.24) 


dp)  p 


dp ) p\dp) T 


(A.25) 


Therefore, 
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](^Y  / 

'  fdp\ 

/  Ur  jp/ 

[dp) 

(A.26) 


For  a  pressure-explicit  equation  of  state,  the  specific  heat  at  constant  volume  can 
be  derived  using 


a  = 


de 


dTj, 


'1 


_  /de 


dp 


(A.27) 


/p^  P  \u±  -  /  p 

Therefore,  the  constant  pressure  specific  heat  can  be  obtained  by  substituting 
Eq.  A.27)  into  (A.26) 

P 


\f^)  (ep\ 


r, 

XOp/T 


(A.28) 


fugacity 

Fugacity  is  particularly  important  when  considering  mixtures  at  equilib¬ 
rium.  Fugacity  /  is  defined  herein  as[30] 


with  the  requirement  that 


ItT d(in  f)'jr  =  dgx 


lim  (—  1  =  1 

p-^0  \p, 


(A.29) 


In  a  mixture,  the  fugacity  of  component  i  can  be  expressed  in  a  similar  form. 

TZTd{lnfi)T  =  {dgi)T  i  =  1--.A'  (A. 30) 

along  with 
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=  1  i  =  l...N  (A. 31) 

\XipJ 

The  fugacity  approaches  the  partial  pressure  of  component  i  in  an  ideal  gas  mixture 
as  pressure  approaches  zero.  Integrate  Eq.  (A. 30)  at  constant  temperature  from 
the  reference  state  to  state  of  interest  to  obtain 

TIT d{\n  fi)T  =  f Q  i<igi)T  (A. 32) 

therefore 


1  /  0\ 

=  ~  > 

_  ^  Y  ^ 

TZT  _\dni) T,p,nj^^,...  \dni/  T,p,nj^^,... 

where  fii  is  called  the  chemical  potential.  Therefore,  the  fugacity  can  be  evaluated 
via 


j _ ^  r  _P  j  ■ 

TZT  drii  UpO  p2 


(A.33) 
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ENTHALPY  OF  VAPORIZATION 


where 
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where 
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APPENDIX  C 


JACOBIAN  MATRICES 


C.l  One-Dimensional  Analysis 
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where 

and  M  is  the  molecular  weight  of  the  mixture. 


(C.l) 
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/  0  0  0  0 


0  \ 


0 

0 

0 


/f  0 

UrX  A 


Cp  Cp 


0 

Ri 

0  pVlr, 


0 

Rn-1 

0 


0  0  0  0 


P'^{N  —  1)171  j 


where 

Ri  —  —  hj\j)  +  pT^i7n{hi  —  hj^)  (C.2) 

L/p 

c.2  Two-Dimensional  Analysis 

The  number  of  species  components  involved  in  the  system  is  assumed  to  be 
N .  The  Jacobian  matrices  are  given  below 
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where  M  is  the  molecular 

weight  of  the  mixture  and 
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The  convective  term  Jacobians  A  in  the  ^  direction  is  written  as 
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0  0  0  0  0  ...  0  \ 
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0  R42  R43  R44  R45  •  •  ■  R46  ’ 

0  0  0  0  i?55  •  •  •  0 

0  0  0  0  0  ...  Ree 


3«i  +  0-2  )f^ 


=  TOLlOioH 


-aiOi2H 
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Rzz  = 
Ri2  = 
Rii  = 

7^44  = 

Ra5  = 
R46  ~ 
R55  = 

7?66  = 


«1  +  ^"2)^ 


Qaj  +  ai^i^u  +  ^aia2fiv  -  +  0:2)^- 


2  ,  '^  „2 


-a\a2iiu  +  ai  +  -(^2  ”  l«i  +  “2 


p^hiVi  -  /ijv2?wj  ~  ~  ^0 
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(C.6) 

(C.7) 

(C.8) 

(C.9) 

(C.IO) 

(C.ll) 

(C.12) 

(C,13) 


where  U  =  ctiu  +  C(2V,  ol\  =  ^xi(^2  —  ty  The  Jacobian  matrices  B  and  R^tj  in 
the  Tj  direction  is  obtained  by  substitute  oi  with  rjx  and  02  with  riy  in  matrices  A 
and  respectively. 


APPENDIX  D 


SIMPLIFIED  ANALYSIS  OF  DROPLET 
VAPORIZATION  RESPONSE 


The  characteristics  of  unsteady  droplet  vaporization  are  examined  here  by 
means  of  a  simplified  analysis.  A  droplet  with  a  mass  m  and  a  mean  lifetime  fy  is 
vaporizing  in  a  forced  oscillatory  environment.  The  variations  in  vaporization  rate 
depend  on  the  difference  between  the  vapor  pressure  at  the  droplet  surface  and 
the  ambient  pressure.  However,  the  vapor  pressure  is  mainly  determined  by  the 
droplet  surface  temperature.  Heat  transferred  into  the  droplet  is  induced  by  the 
temperature  difference  between  the  ambient  gases  and  liquid  droplet.  There  are 
several  assumptions  involved  in  the  current  analysis.  First,  the  flow  is  treated  in 
a  one- dimensional  spherical  coordinate  system.  Second,  the  liquid  temperature  is 
considered  to  be  uniform  over  the  entire  droplet  but  varies  with  time.  The  mass  dif- 
fusivity  is  assumed  to  be  inversely  proportional  to  ambient  pressure.  Furthermore, 
a  quasi-steady  theory  is  needed  to  obtain  the  droplet  lifetime. 

Vaporization  Rate: 

The  droplet  vaporization  rate  can  be  determined  by  the  time  variation  of 
total  droplet  mass 


w  =  — 


dm 

lit' 


(D.l) 


where  m  and  w  are  the  instantaneous  mass  and  vaporization  rate,  respectively. 
Using  the  perturbation  concept,  Eq.  (D.l)  can  be  developed  into 
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^  —  —  \  /—  — 

— — (m  m  +  m)  =  w  w  +  w 

at 


After  time  averaging  the  above  equation,  the  mean  vaporization  rate  is 


(D.2) 


w  = 


dm 

dt 


(D.3) 


Substitute  Eq.  (D.3)  into  (D.2)  to  obtain 


m  dm' 
w  dt 


m  —  — w 


(D.4) 


or 


(1  —  Tys)m^  = 


(D.5) 


where  s  is  the  Laplace  operator  and  is  the  droplet  lifetime  and  can  be  approxi¬ 
mated  by  dividing  the  instantaneous  mass  with  the  vaporization  rate 


m 


= 


w 


(D.6) 


Assuming  that  the  period  of  oscillation  is  much  larger  than  the  local  gas  residual 
time.  Using  the  quasi-steady  state  theory,  the  droplet  lifetime  can  be  evaluated 
by  the  classical  d^-law. 


m 


Tv 


w  2Ps/95  ln(l  4- jB)’ 


(D.7) 


The  quantities  r,  pi,  ps,  and  Vs  are  instantaneous  droplet  radius,  liquid  density, 
surface  density  and  mass  diffusivity,  respectively.  The  variable,  S,  is  called  the 
Spalding  transfer  number  and  is  defined  as 


165 


R  —  '^s)  _  y^F,00  ^F,S  /£)  g\ 

AK  - 1 

Since  the  surrounding  gas  is  pure  nitrogen,  Ypoo  can  be  set  to  zero.  Other  quantities 
Ypsi  Too,  Ts,  Cp,  and  A/i^,  are  the  n-pentane  concentration  at  the  droplet  surface, 
ambient  temperature,  surface  temperature,  specific  heat  on  droplet  surface,  and 
latent  heat  of  vaporization,  respectively. 

The  vaporization  rate  predicted  by  the  quasi-steady  state  theory  can  be 
expressed  as 


w  =  ^irrlpsVs  =  ^T^rsPsT>sYi4>.  (D.9) 

At  low  to  moderate  pressures,  the  diffusion  coefficients  vary  inversely  with  density. 

Assuming  Vg  on  —  and  r  oc  Eq.  (D.9)  can  be  further  developed  into 

Ps 

w  =  +  B),  (D.IO) 

where  Ci  is  a  constant.  Again,  as  before,  apply  the  small  perturbation  technique 
to  this  equation,  a  linearized  equation  is  obtained. 
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Vapor  Pressure: 

The  vapor  pressure  on  the  droplet  surface  can  be  evaluated  by  Antoine’s 
equation. 

\nps  =  ai-—^ —  (D.13) 

Ts  -  az 

where  ai,  02  and  03  are  constants  depending  on  the  mixture  of  concern[34].  The 
linearized  perturbation  equation  is  found  to  be 

p's  ~  (D.14) 

where  b  is  defined  as 


a2Ts 

{fs  -  az)^ 


(D.15) 


Energy  Balance  on  Droplet  Surface: 

The  conservation  law  of  energy  must  be  satisfied  at  the  droplet  surface. 


Qs—^£  =  Qg-*s  —  Qs-i-g  (D.16) 

Using  the  perturbation  technique  together  with  the  quasi-steady  assumption,  the 
linearized  equation  is  written  as 

Qs^eQ’s^i  =  -  Q's-,g)  (D.17) 

Heat  Transfer  from  Ambient  Gases: 

The  heat  conducted  into  the  droplet  is  treated  as 

Qg^s  -  4xr^/i(roo  -  Ts)  (D.18) 


Therefore, 
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Q',^,  =  \m’  -  (D-W) 

O  loo  Is 

Heat  Carried  with  Vaporized  Fuel: 

Energy  convected  away  from  the  droplet  can  be  obtained  by  evaluating  the 
enthalpy  increase  in  the  vaporized  fuel. 

Qs^g  =  AhyW.  (D.20) 

The  perturbed  equation  can  be  approximated  by 

gU,  =  AK  +  w' 

=  diPoo  +  d2Tl  +  w' 

where  di  and  ^2  are  scaling  constants  and  take  the  form 

i 

I 

di 

d2 

Energy  Transferred  into  Droplet: 

Lumped  Capacity  Model: 

Under  the  assumption  that  the  droplet  is  treated  as  a  lump  of  substance 
with  an  infinite  thermal  conductivity  (A^  — >  oo).  The  transient  temperature 
variation  is  determined  by  formulating  an  overall  energy  balance  for  the  droplet. 

r)T 

(D.24) 

where  Cp^e  and  qs-^e  are  the  liquid  specific  heat  and  heat  flux  conducted  into  the 
droplet  interior,  respectively.  The  linearized  perturbation  equation  is  written  as 


^dAK\ 
~  \  dT  ) 


dT 

dAhv' 

dpoo  J  Ts  AK 


PCG  Ahy 
Poo 


(D.22) 

(D.23) 


(D.21) 
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=  Qs^eQ's-.e  (D-25) 

Combine  Eqs.  (D.5),  (D.ll),  (D.14),  (D.17),  (D.19),  (D.21),  and  (D.25),  and  sub¬ 
stitute  dt/t  with  the  Laplace  parameter,  s  =  iw,  the  droplet  vaporization  response 
function  can  be  obtained  as 

/  -3/?(l  -  iTyw)  {bdi  +  '^!r)  + 

Rp^^  =  T - r-S - ^ - T  (D.26) 

Poo  2^  Jl  +  3TyW^J  +  6/?  +  Z  2^^/  “  S^RTyW  —  3bl3TyW 

where 
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—  tan' 


■i 


2'^R  +  dTyW^I  +  bl3 

2^/  —  —  Zh^TuW 


(D.32) 


Finite  Diffusion  Model: 


The  transient  diffusion  response  of  energy  transferred  from  the  ambiance 
can  be  approximated  by  solving  the  transient  one-dimensional  heat  conduction 
equation  in  a  spherical  coordinate  with  a  forced-oscillatory  surface  temperature. 


(D.33) 


a  dt  r'^  dr  \  dr 

Assuming  the  droplet  temperature  is  uniformly  distributed  in  space  but  varies  with 
time.  The  initial  temperature  is  selected  to  be  equal  to  the  surface  temperature 

t  =  0  ,  T  =  fs 

A  prescribed  temperature  oscillation  is  first  specified  at  the  droplet  surface  to 
simulate  the  dynamic  behavior  of  energy  transfer  into  the  liquid  phase, 

r^R  ,  T  fs  +  \fsT^\e^^^ 
r  =  0  ,  T  =  Ts, 

where  u>  is  the  angular  frequency.  To  solve  this  equation,  four  dimensionless  pa¬ 
rameters  are  introduced  herein, 

"  Oit  ^  T  —  Ts 


r  = 


0 


w*  = 


o 

r  w 


R'  ‘  r2’  ^  iTsT^I  ’  a 


(D.34) 


where  R  is  the  initial  droplet  radius  and  a  is  the  thermal  diffusivity  in  the  liquid 
phase.  Thus,  Eq.  (D.33)  becomes 


dQ  1  d 


,90 


dr  T]^  dr]  L  dr] . 

subject  to  the  initial  and  boundary  conditions 


(D.35) 
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T  =  0  ,  0  =  0 

Ti  =  l  ,  0  =  e*‘"*^ 

r;  =  0  ,  0  =  0.  (D.36) 


Therefore,  Eq.  (D.35)  can  be  solved  to  obtain 


0  =  - 


sink 

\/w*/2(l  +*)?? 

sink 

yjuj* I2{1  -1-  i) 

(D.37) 


This  equation  provides  the  temperature  profile  and  can  be  used  to  calculate  the 
energy  flux  variation  at  the  droplet  surface. 


Qs-^eQ's 


s-*£  - 


dT 


dT 


r=rs 


=  iirrsXiTslT^l  ^ 
=  A7rrsXef,F{w*)T', 


Tl=l 


(D.38) 


where  is  the  liquid  thermal  conductivity,  and 


F{w*)  =  {iw*)hoth{iw*)^  -  1 


Combine  Eqs.  (D.5),  (D.ll),  (D.14),  (D.17),  (D.19),  (D.21),  and  (D.38),  and  sub¬ 
stitute  dt/t  with  the  Laplace  parameter,  s  —  iw,  the  droplet  vaporization  response 
function  can  be  obtained  as 


w'  _ 

—3/9(1  —  ir-uw) 

{bAhywdi  -f  +  i'^i 

Poo 

2^r  -I-  3r„u;^/  +  bl3Ahj,w 

+  i 

2^*/  —  —  ZTvwb/SAhvyf 

(D.39) 


where 
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'^R  =  iTrR\eTsRe[F{w*)]  +  AhJ-^^)+d2AK^ 

\J~oo  ^  S  '' 

=  ATrRXiTsIm[F{w*)] 

The  magnitude  and  phase  angle  of  the  response  function  Rp  are 
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2’I'7  —  dTyW'^R  —  ZTyWb^Ahy\ 

Qj)  =  tan  ^(— rt,u;)  +  tan  - r - — 'l 

^  V  i-  ;  \bAKwdiF'>l>Rj 
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(D.40) 

(D.41) 


(D.42) 


(D.43) 


